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SUMMARY 


Direct measurements of supersonic local skin friction, using the 
floating-element technique, are presented for Mach Numbers from 
2.0 to 4.5 and Reynolds Numbers from 3 X 10°to9 & 10% Tur 
bulent flow and transition are emphasized, although some meas 
urements in the laminar regime are included. The observed effect 
of compressibility is to reduce the magnitude of turbulent skin 
friction by a factor of two at a Mach Number of 4.5 and a 
Reynolds Number of about 10° 

The boundary-layer momentum-integral equation for constant 
pressure is verified within a few per cent by two experimental 
methods. Typical static pressure measurements are presented to 
show that transition can be detected by observing disturbances 
in pressure associated with changes in displacement thickness 
of the boundary layer. 

It is found that the turbulent boundary layer cannot be defined 
experimentally for values of “,0/v; less than about 2,000, where @ 
is the momentum thickness. For larger values of 1)6/»; there is a 
unique relationship between local friction coefficient and momen 
tum-thickness Reynolds Number at a fixed Mach Number. The 
Appendix compares the present measurements at J/ = 2.5 with 


experimental data from other sources 


INTRODUCTION 


by FLOATING-ELEMENT TECHNIQUE, which attempts 
to measure directly the friction drag on a portion of 
a body in a fluid flow, was successfully used in 1929 by 
Kempf! and later by Schultz-Grunow.” An early but 
largely unproductive effort to apply the technique at 
supersonic speeds was that of Ginsburgh.* Beginning 
in 1948, H. W. Liepmann at GALCITT undertook the 
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development of a floating-element instrument for high 
speed flow, and this work has been carried forward by 
Dhawan?‘ and by Hakkinen.® The method has also been 
used, with variations, by Weiler and Hartwig® at the 
University of Texas, by Eimer’ at GALCIT, by Brad 
field, DeCoursin, at the University of 
Minnesota, and with conspicuous success by Chapman 
and Kester’ at the Ames Aeronautical Laboratory of the 
NACA, the last two of these investigations being made 
with a body of revolution rather than with a flat-plate 


and Blumer‘ 


model. 

The writer’s research at JPL,* 
the success of Liepmann and his group, is an independ 
ent attack on the flat-plate drag problem. During the 
research full advantage was taken of the size and range 


* although inspired by 


of the 20-in. supersonic wind tunnel at the Jet Propul- 
sion Laboratory and of Liepmann’s experience with the 
floating-element technique. Design of the flat-plate 
model and instrumentation was begun in early 1951, 
and the experiments discussed here were carried 
out in the spring of 1952. The 
has been briefly reported at the 1952 Annual Summer 
the the Aeronautical Sci 
ences in Angeles and at the Eighth Inter 
national Congress on Theoretical and Applied Me 
chanics in Istanbul, and some preliminary data have 
A comprehensive dis 


present research 


Meeting of Institute of 
Los 


been published in a short note." 
cussion of the research can be found in 
thesis,'! which has been published in three reports of the 
Jet Propulsion Laboratory.'*~'* Some introductory 


remarks in the thesis have also appeared as a survey 


the writer's 


paper.” 
List OF SYMBOLS 


Cc = local friction coefficient 
Cy = mean friction coefficient 


** Jet Propulsion Laboratory, California Institute of Tech 


nology. 
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Mf = free-stream Mach Number 

p = pressure 

q = free-stream dynamic pressure, p)M™,"/2 
R = ideal Reynolds Number 

R’ = Reynolds Number per unit length, 4,/» 


R, = ideal Reynolds Number corresponding to point of ob 
served maximum surface shearing stress 


Rg = momentum-thickness Reynolds Number, 1)6/»; 
7 = temperature 
u = velocity 
x = distance from plate leading edge 
y = distance from plate surface 
y = ratio of specific heats 
6 = boundary-layer thickness 
6* = boundary-layer displacement thickness 
0 = boundary-layer momentum thickness 
# = momentum-thickness increment of tripping device 
“ == viscosity 
v = kinematic viscosity, u/p 
p = density 
T, = wall shearing stress 
Subscripts 
Q = value at stagnation condition 
1 = value in free stream 


value for 14 = 0 at fixed R 


1 = 


MODEL AND INSTRUMENTATION 
(1) Flat Plate 


Fig. 1 shows the flat-plate model installed in the test 
section of the 20-in. supersonic wind tunnel at the Jet 
Propulsion Laboratory. The plate completely spanned 
the tunnel in order to isolate the working surface from 
disturbances originating elsewhere on the model. Three 
floating-element instruments were mounted in the plate, 
on the tunnel plane of symmetry, at x = 5.5, 13.0, and 
24.0 in. from the leading edge. Three small windows of 
'/>-in. plate glass were installed in each door, so that 
flow over the plate surface and leading edge could be ob- 
served using the tunnel schlieren system. 

A series of static-pressure orifices was provided on the 
plate surface about 3 in. from the tunnel centerline. 
In addition, each of the floating-element instruments in- 
corporated three orifices for static-pressure measure- 


Fic. 1. The flat-plate model in the 20-in. supersonic wind tunnel. 
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ment; one of these orifices was the gap around th 
element itself. Three copper-constantan thermo 


couples were soldered to stainless-steel buttons buried 
in the plate surface near the floating-element stations 
and the temperature was also measured in the element 
balance chambers and in the tunnel reservoir. 

All metal components exposed to the air stream on 
the working surface of the plate were fabricated of stain. 
less steel. Particular attention was paid to working 
surface finish and to the sharpness of the leading edge 
which was a wedge of 15-deg. included angle. Various 
methods of estimating the leading-edge radius indicated 
that the value finally achieved was less than one half- 
thousandth of an inch. 


(2) Floating Element 


The skin-friction instrument used at JPL was similar 
in principle to the floating-element mechanism of Liep- 
mann and Dhawan.‘ The most important difference, 
aside from size and range, was the use at JPL of a null 
technique in measuring applied forces. 

The basic components of the instrument are shown in 
Figs. 2 and 3. The inner flexure system, consisting of 
the table, two flexlinks, and the ring, allowed a deflec- 
tion of the element under an applied load. The deflec- 
tion was measured indirectly as the displacement that 
must be given to the table in order to return the element 
to a standard null position with respect to the surround- 
ing structure. Motion of the table was effected by 
means of a precision micrometer, acting through a piano- 
wire push rod clamped at both ends, and was con- 
strained by the second or outer flexure system shown in 
the figures. 

The null position of the floating element was de- 
tected by a Schaevitz variable-reluctance _ trans- 
former!® whose coil was fixed on the instrument struc- 
ture and whose core was attached to the element ring. 
Two dashpots were provided for damping in the inner 
flexure system and were adjusted to give slightly less 
than critical damping. The choice of a rectangular sur- 
face element was made in order to approximate a local 
measurement, and the corresponding rectangular hole, 
'/, by 1'/. in., was fabricated by a combination of 
milling, machine-filing, and broaching operations. 


(3) Micrometer 

The micrometer shown in Figs. 2 and 3 was designed 
around a standard commercial lead-screw of 0.025-in. 
pitch. The lead-screw was keyed to prevent rotation 
and was actuated by a rotating nut that was driven 
through a worm and worm gear. A counter and 
graduated dial on the micrometer input shaft indicated 
the position of the lead-screw in millionths of an inch. 

Each of the micrometers was calibrated by direct 
comparison with a Pratt & Whitney standard measur- 
ing machine. It was found that the maximum microm- 
eter position error did not exceed 0.5 per cent and st 
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MEASUREMENTS OF TURBULENT FRICTION ¢ 
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ELEMENT 


Fic. 2. The floating element and micrometer installed in the plate. 
could be neglected in the present measurements. In diameter, spaced '/, in. apart, and projected about 1/1, 
general, under steady conditions, the null position as in- in. beyond the surface of the plate. For increasing 
dicated by the Schaevitz gage could be repeated within © Mach Number, there was a small but noticeable increase 
one or two millionths of an inch using ordinary care to. in the wire Reynolds Number needed to trip the 
avoid backlash in the mechanism of the micrometer. boundary layer. 

Failure to obtain the sensitivity mentioned could The device of discharging air from small holes in a 
usually be traced to the presence of foreign material in surface was proposed by Fage and Sargent” as a 
the clearance gaps. The spring constant for the inner method of stimulating transition to turbulent flow in 
flexure system was about 0.4 Ib. per in., and deflections the boundary layer. The most important property of 
ranging from 0.0005 to 0.015 in. were encountered in such a tripping device, aside from the obvious ad- 
vantage of control, is the small intrinsic drag of the jets. 
In the JPL tests, air was discharged from a row of holes 
(4) Tripping Devices 0.014 in. in diameter, spaced '/, in. apart, and located 
3/,in. from the leading edge. During some preliminary 
measurements at a Mach Number of 2.0, it was found 


the wind tunnel. 


During the JPL research, three different tripping de- 
vices were used to stimulate transition to turbulent 
flow. These were a sand strip, a leading-edge fence, and 
a row of airjets. The fence, consisting of a row of wire 
posts normal to the airflow and to the plate surface, was 
used in most of the measurements of turbulent friction. 
The sand strip was tested at J = 4.5, and the effect of 
air discharge on transition was investigated at Mach 
Numbers of 2.0 and 4.5. Fig. 4 shows the appearance 


of the boundary-layer flow at J = 4.5 with natural 
ELEMENT 


transition and with the various tripping devices. | 
The fence was found to have the virtue that the 
initial momentum loss at the wires could be computed ne 
from the known drag characteristics of cylinders in 
supersonic flow. However, there is some doubt as to 
es % ¢ DASHPOT 


whether true equilibrium was reached in the turbulent 
boundary layer far downstream, where the fence usu- 
ally contributed from 10 to 20 per cent of the total 
momentum loss. The wires used were 0.014 in. in Fic. 3. Separated view of the floating element. 
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that there was a critical discharge rate, which did not 
appear to depend on tunnel pressure. For 
rates than the critical value the airjets had a negligible 


effect on the distribution of surface friction, while for 


larger rates the effect was large and did not change 
appreciably with a further increase in the strength of 
the disturbance. At ./ = 2.0, the critical decrement in 
displacement-thickness Reynolds Number was about 10 
for values of 1,6*/», between 700 and 2,000 at the sta- 
tion of the jets. At JJ = 
transition was less sensitive to discharge rate, and the 


critical decrement in 1,6*/»,; was about 25 for initial 


1.5, the boundary-layer 


values between 3,000 and 5,000. 


FRICTION MEASUREMENTS 


(1) Local Friction 


The direct friction measurements of the present re- 
search are shown in Figs. 5-14. Both the ordinate 
Tx/q and the abscissa m,x/v; have been deduced in an 
unambiguous way from measured physical quantities. 
In particular, the static pressure in the balance cham- 
bers was used in all cases to determine the local free- 
Mach Number, temperature, velocity, and 

For this calculation the fluid was assumed to 


stream 
density. 
be a perfect gas expanding isentropically from the tun- 
nel reservoir with y = 1.400, and the Sutherland equa 
tion? was used to relate temperature and viscosity. 
Two corrections, always amounting together to less 
than 2 per cent in the shear, have been applied to the 
data of the figures. First, streamwise pressure forces, 
caused by small local pressure variations in the element 
gaps, have been compensated for on the basis of a chan- 
nel calibration discussed in reference 15. Second, the 
turbulent friction data, which are relatively sensitive to 
variations in Mach Number, have been corrected to the 
nominal Mach Numbers of the figures by means of the 


approximate relationship 


(0C,/OM)p = —m?C,/M 


where 
m? = 


That is, the function C,(.\7, R) at constant R has been 
approximated, for this calculation, by the quantity 
C,(0, R) 


VT,/To CAO, R) = [1 + (y — 1) M?2/2] 


Random error in the local turbulent friction co- 
efficients is believed to be less than 2 per cent. Con- 
tributing to this error are the uncertainty in the element 


electrical zero position under dynamic conditions and 


t The parameters C, 7, and yu, in 

UT, + C)/(T + ©)| 

have been taken as 192°R., 492°R., and 3.59 & 1077 (slugs/ft.' 
(ft./sec.)(ft.), respectively. 


u/u, = (T/T,) 


AERONAUTICAL 


smaller 


SCIENCES FULLY, 18:54 

the uncertainty in both the element static calibratio; 
and the micrometer calibration. In addition, the fric 
tion measurements at large Mach Numbers required 
the accurate measurement of tunnel static pressures of 
Both the pres 
sure and the zero position are, of course, relatively mor: 


the order of 2 to 10 mm. Hg absolute. 
uncertain at low Reynolds Numbers. It was possible, 
during the profile surveys mentioned below, to repeat 
some of the shear measurements at the first and second 
stations. These data are shown as flagged points in 
Figs. 5-14, and in general indicate excellent repro- 
ducibility. The only systematic error that is believed 
to be important arises because flow over the clearance 
gaps may contribute to the force on the element; this 


error probably does not exceed 3 per cent.!* 


(2) Mean Friction 


The execution of boundary-layer research with im 
pact-probe instrumentation will not be discussed at 
greater length than is necessary to establish the ap- 
proximations that are usual in an analysis of the data 
It should be noted, however, that the experience of 
many investigators has shown that an impact probe 
in a boundary layer is an instrument subject to several 
sources of error. The presence of the probe may in 
fluence the development of the original boundary layer 
through upstream propagation of pressure disturbances, 
especially in the case of laminar flow. The existence of a 
velocity gradient and the proximity of a wall will lead 
to distortion of the flow pattern around the probe. On 
the other hand, the use of small probes, in an effort to 
minimize the problems already mentioned, may in- 
troduce new errors associated with viscous dissipation 
near the probe nose. Finally, the effect of velocity and 
density fluctuations on the pressure indicated by the 
probe may not be negligible. 

In the present experiments the quantity sought is the 
momentum thickness, defined by 

p u 


m u 
/ ( l— ) dy (1 
J 0) pi Uy it} 


Fig. 15 shows the variation with y of the integrand in 


d= 


this expression for two typical boundary-layer profiles 
ata Mach Number of 4.5. 
ured was, of course, the impact pressure; 


The quantity directly meas 
the stati 
pressure was assumed to correspond to isentropic ex 
pansion from the tunnel reservoir with y = 1.400 and 
was taken to be constant through the boundary layer 
The latter assumption was in each instance confirmed 
by measurement of surface pressure near the probe 
station. From these pressure data the local Mach Num 
ber in the boundary layer follows. 

The stagnation temperature in the boundary layer 
has been assumed constant at the value in the tunnel 
reservoir. The local static temperature, density, ve- 
locity of sound, and fluid velocity were obtained from the 
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Fic. 14. Local friction, air-jet trip, 7 





usual equations of energy and state.+ An inspection Figs. 16 and 17 show typical pressure distributions 
of the method of data analysis shows that errors in the observed during the friction measurements of Figs 
velocity u and the momentum pu will be, respectively, 9 and 10 at a nominal Mach Number of 3.7. The ex- 
in the same and opposite direction as the error in the pansion downstream of the leading edge and the com- 
local stagnation temperature and will be of half the pression in the transition region are conspicuous when 
magnitude. Although the calculation of the integrand the observed pressures are compared to the empty test- 
in Eq. (1) thus involves a cumulative error at a given section pressure distribution. 

point, it seems unlikely that the net error in @ from this (4) Transition 
source might be larger than the probable uncertainty of ; 7: ; 
1 The region of transition from laminar to turbulent 
flow may be identified by a minimum and a maximum 
1 the curve of local shearing stress plotted against 


a few per cent in the stagnation temperature which is 

typical for most of the boundary layer. 
it 

Reynolds Number. The transition is not in general a 


(3) Static Pressure 
unique function of the Reynolds Number for the ex- 


Before discussing the friction measurements in de- : Soe ; wre 
periments reported here, and this fact increases the 
difficulty in determining the shear distribution from ob- 


servations at three points. However, the strong corre- 


tail, it is convenient to describe a second technique that 
was found useful in establishing the nature of the 
boundary-layer flow over the plate. During the present ; baie 
ane ae <i ond, lation observed between shear and pressure distribu- 
research, relatively large static-pressure disturbances .. é Sea es : ge 
: tions is useful in interpolating between measured values 
of the friction coefficient at a given stagnation pressure 
Figs. 18-24 show, for the tests with natural transition 
and with the leading-edge fence, Reynolds Numbers 


were observed on the plate surface, especially in the 
region of transition from laminar to turbulent flow. 
These disturbances are believed to be connected with 
the apparent body shape which the boundary-layer dis- 


. : of the minima and maxima in the shear and pressure 
placement thickness presents to the airflow. 


distributions as obtained from Figs. 5—12 and from pres- 
t The gas constant ® has been taken as 1,716 (lb./ft.2)/(slug/- SUTe measurements like those of Fig. 16. Some judg- 
ft.*)(°R.) in the state equation p = pRT. ment has, of course, been exercised in estimating the 
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MEASUREMENTS OF 


Reynolds Numbers for the pressure extrema, since the 
data are not always adequate to the demands imposed. 
Furthermore, for technical reasons the pressure meas- 
urements were made approximately 3 in. from the tun- 
nel plane of symmetry, and it is not certain that the 
boundary-layer flow was accurately two-dimensional. 
In spite of these defects in the pressure data, it appears 
that the minimum in the static-pressure distribution 
corresponds closely in all cases to the position of maxi- 
mum rate of change of the local friction coefficient. It 
should be remembered that both the friction and pres- 
sure measurements in the region of transition represent 
the response of a highly damped instrument to a quan- 
tity that is fluctuating rapidly with time and, therefore, 
that the present experience with transition may be 
peculiar to the particular wind tunnel and instrumenta- 


tion involved. 


NATURAL TRANSITION 


FENCE TRIP 


SAND TRIP 


TURBULENT 


FRICTION ON A FLAT PLATE 139 
Several other features of Figs. 18-24 are of special 
interest. The Reynolds Number per inch, R’, is pro- 
portional to the tunnel stagnation pressure, and the 
movement of the natural transition toward smaller 
values of x but toward larger values of ux»; with in- 
creasing tunnel pressure is evident in the coordinate 
system of the figures. On the other hand, the fence 
tripping device apparently has the property at high stag- 
nation pressures of causing transition to occur at a fixed 
Reynolds Number increment from the leading edge. It 
is therefore possible that the tendency of the correspond- 
ing local-friction measurements to define a single curve 
in the turbulent flow regime is not accidental (see, for 
example, Fig. 8, where the measurements are relatively 
insensitive to the presence of local pressure gradients 
At low values of tunnel stagnation pressure, however, 


the effect of the fence is entirely different. The large 











Fic. 4 


Schlieren photographs of the boundary layer at \J = 


4.5, R’ = 


130,000 per in 



































440 JOURNAL OF THE AERONAU'’ 
10 7 | 
TESTS 62, 76 | 
08 
| 
06 - 
Ee | (1 - *) 
P, Uy U, | 
04 
02+—— 
Oo 
(@) 
y , INCHES 
Fic. 15. Typical momentum-thickness determination, 1/7 = 4.5. 


initial increment in @ caused by wire drag apparently 
does not strongly influence the slope of the velocity pro- 
file at the wall upstream of actual transition, as may be 
verified by a comparison of Figs. 11 and 12. Further- 
more, at the lowest tunnel pressures, the transition may 
actually be moved to higher Reynolds Numbers by the 
disturbance from the fence, from the evidence, for ex- 
ample, of Figs. 23 and 24. These remarks serve to 
emphasize the lack of real understanding of the transi- 
tion phenomenon and show the need for caution in at- 
tempting to predict the influence of upstream dis- 


turbances. 
INTERPRETATION 


(1) Uniqueness 


In a survey paper,’ the writer has discussed the 
question of uniqueness for the turbulent boundary layer 
with constant density, emphasizing the phenomenologi- 
cal point of view. Unfortunately, there exists at present 
no adequate basis for a generalization of this study to 
compressible flows. The present discussion will therefore 
depend entirely on experimental information in attempt- 
ing to determine the properties of the ideal boundary 
layer. 


Experience in low-speed flow suggests that the tur- 
bulent boundary layer with uniform pressure may be 
characterized as fully developed when there is a unique 
relationship between the slope and ordinate of the 
momentum-thickness distribution 6(x). That is, the 
assumption of uniqueness implies a function 


C, = 2(dRe/dR) = g(Re) (2) 
which may be integrated, on satisfying the initial con- 
dition R = 0 when Ry = 0, to yield a relationship 
R = G(R.) for the ideal boundary layer. Since the 


surface friction is supposedly always positive, Re is a 
monotonically increasing function of R, and the initial 
condition R = Rs = 0 can alwavs be satisfied. 


FICAL 
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However, the integration can be carried out, and an 
exact definition of R achieved, only if the function g of 
Eq. (2) is known for all positive values of Ro. The 
function in question being defined experimentally, it js 
clear that any corresponding definition of R is not itself 
unique and can be valid only in an asymptotic sense 
In other words, various experimental determinations of 
R become equivalent only when Ry is large compared 
with values below which the turbulent boundary layer 
is not observed and is therefore undefined. It also 
follows, when considering various measurements of tur- 
bulent friction, that discrepancies in the final corrected 
data may not mean that one or another series of 
measurements is in error but only that they have been 
differently interpreted. 

Fig. 25 shows the relationship between 7, g and 16 /y 
at various Mach Numbers, as determined experimen- 
tally at JPL, together with the corresponding curve for 
zero Mach Number from reference 15. It appears that 
the supersonic boundary layers in question are fully de- 
veloped according to the present criterion, except per 
haps for the points that are farthest to the left in the 
figure at the lower Mach Numbers. In particular, two 
of the profiles at 17 = 4.5 were obtained with the fence 
trip, two with the airjets, and one with the sand strip. 


(2) Maximum Shearing Stress 


One interpretation of the local friction data of Figs. 
5-12 is suggested by the fact that relatively detailed in- 
formation on transition is available in Figs. 18-24. For 
the measurements in question, it has been observed by 
the writer that the increment in Reynolds Number be- 
tween the point of maximum shearing stress and the 
point of measurement is a well-defined function of the 
measured local friction coefficient, whether transition is 
Fig. 26 shows this relationship for 
The abscissa 


natural or forced. 
Mach Numbers of 2.6, 3.7, and 4.5. 
(R — R,) in the figure is the quantity ,Av »,, where 
Ax is measured between the point of maximum shearing 
stress and the downstream element station for each 
operating condition. Data at a Mach Number of 2.0 
are omitted because small local pressure gradients 
affect the surface friction and simultaneously make it 
difficult to identify the pressure disturbance at transi- 
tion. 

The local friction coefficient is experimentally a 
unique function of (R — R,) at a given Mach Number 
and also by assumption a unique function of R. It 
follows that R, can only be a function of R for constant 
AM. However, since the measurements at different sta- 
tions at a fixed value of tunnel pressure should be char- 
acterized by the same value of R,, the latter must be a 
constant for each of the three curves of Fig. 26. This is 
to say that the end of the transition region, when re- 
ferred to the ideal boundary layer, occurs at a value of 
C,, or of Re, or of R, which depends only on 1/7 whether 
transition is stimulated or not. 
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Natural transition, 1fJ = 3.70. Fic. 22. Transition with fence trip, 17 = 3.70. 
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Fic. 23. Natural transition, 1J = 4.54. 


(3) Momentum Balance 


The argument of the preceding section is clearly in- 
complete. Neither the question of uniqueness nor the 
question of momentum balance can be resolved experi- 
mentally without reference to momentum-thickness 
measurements, and these have not been cited in develop- 
ing Fig. 26. 

The momentum-integral equation for two-dimen- 
sional flow without pressure gradient is 


Cy, = t/q = 2(d6/dx) (3) 


where the measured quantities 2d@/dx and r,,/q¢ repre- 
sent, respectively, the rate of momentum loss in the 
fluid and the drag force on the surface. For the experi- 
ments with the fence trip, Eq. (3) may be integrated in 


the form 


a” 
6 = 6+ = dx (4) 
0 -q 
where 6) is the contribution of the tripping device. For 


example, if the drag coefficient of a cylinder based on 
projected area is 1.2 for the range of J and R con- 
sidered here, then!* @) = 0.0034 in. 

Multiplying Eq. (4) by the Reynolds Number per 
inch, R’, there is obtained 


Re = R’% + / om dR'x 
0 2q 


This expression is to be applied to measurements of Rg 
However, the behavior of 


(5) 


with x fixed and R’ variable. 
the transition curves in Figs. 20, 22, and 24 implies that 
T,/q in the laminar, transition, and turbulent regimes is 





TEST 15 
7 
10 10 
u,x/y, 
Fic. 24. Transition with fence trip, JJ = 4.54 
a function of m,x/», only, at least for large R’. Thus 


Eq. (5) becomes 


“a Ta Ay 
Ro, — Ro, = | ( + ) ar'x 
R' x 2q x 


It is evident at once that the measurements with the 
fence trip must violate the principle of uniqueness ex- 
pressed by Eq. (2). If Eq. (6) is correct, it follows that 
R, cannot be a unique function of 7,/g when R’ is 
variable and 6 is different from zero. It also follows 
that Fig. 26 should be viewed with suspicion, since the 
evidence is against the proposition that a unique value 
of Rs» can be associated with each point in the figure. 

Integration of Eq. (6) may be carried out conven- 
iently in Fig. 26, on replacing the variable of integration 
R'x by the equivalent quantity (R — R,). Table 1 
lists some observed values of (Re, — Re,), together with 
the corresponding quantities computed from the right- 
hand side of Eq. (6) using #. = 0.0034 in. and x = 21.5 
in. The general agreement in the table confirms the 
supposition that 7,,/g is a function of (Re — R’6) rather 
than of Ry for the experiments with the fence tripping 
device. In other words, a residual disturbance from the 
fence is detectable, at least in part, as an additive in- 
crement in Re, and the hypothesis of uniqueness is 
therefore not supported by these data. 

Whether this conclusion is important in the present 
context depends on whether the effect can be distin- 
guished in Fig. 25. In particular, Eq. (6) may be con- 
sidered merely as a formula for experimental continua- 
tion of the direct measurements of Re» given in the figure, 
the continuation being consistent with momentum bal- 


ance without assuming uniqueness. In this sense, Eq. 


TABLE | 


Investigation of Momentum Balance 


(Re. — Re) (Re. — Re,) (Re. — Ro,) (Re — Ri) (Re — R;) 
VW Tests Trip Eq. (6) Eq. (5) Observed Eq. (7) Fig. 25 
2.6 27, 26 Fence 3,540 3,600 4,200,000 4,300,000 
3.7 20, 19 Fence 3,310 3,460 4,650,000 5,050,000 
4.5 23,32 Fence 2,960 ; 3,120 4,600 ,000 5,150,000 
4.5 76, 62 Air jet 2,440 2,340 3,400,000 3,850,000 
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Fic. 26. Equivalence of natural and forced transition. 


(6) allows the data of Fig. 26 to be transferred to the 
coordinate system of Fig. 25, with the result shown in 
Fig. 27. 

Finally, a joint test of the hypotheses of uniqueness 
and momentum balance can be carried out by direct 
calculation from Eq. (5). Choosing a value for R’ of 
290,000 per in., local friction distributions can be 
estimated as shown by the dashed lines in Figs. 7-12. 
The relationship between Ro and r,./g obtained by inte- 
gration of Eq. (5) is found in each case to j oin smoothly 
onto the appropriate experimental curve, as shown in 
Fig. 27. 

For the fence trip, the computed values of Rg at the 
most downstream station are found to be a few per cent 
smaller than values obtained by experimental con- 
tinuation—i.e., than values measured directly. This 
discrepancy, while not important in Fig. 27, should be 
compared to the result of a less equivocal investigation 
of momentum balance at a Mach Number of 4.5 with 
the airjet trip. The relevant profile data at two sta- 
tions at the same tunnel pressure have already been pre- 
sented in Fig. 15. Table 1 lists the experimentally 
observed increment (Ro, — Re), together with the in- 
crement calculated from Eq. (5) with @ = 0; the ob- 
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served increase in Ry between the two probe stations js 


maller than the integral oi 9 


perhaps 4 per cent 
over the same interval. The diser« pancy of 4 per cent 
by coincidence or otherwise, agrees with an estimate 
of the probable effect of the clearance gaps around the 
floating surface element. 

An inspection of the velocity profile data tabulated in 
reference 14 indicates that the profiles with the airjet 
trip are probably free of residual effects from the trip 
ping device and may in fact be taken as typical of th 
supersonic turbulent boundary layer. Returning to Fig 
27, it follows that the fence data probably do not de 
viate from the fully developed condition by more than 
the observed discrepancy between the fence and airjet 
data at M = 4.5. 
therefore certified as representing the fully developed 


The collected points in Fig. 27 are 


turbulent boundary layer, within the degree of accuracy 
implied both by the comparison just given and by the 
estimate already made of possible systematic error in 


measurement of shearing stress. 
(4) Reynolds Number 


One method of obtaining Reynolds Numbers for the 
present data requires the supersonic measurements to be 
compatible with measurements at zero Mach Number 
and to be internally consistent in terms of the momen- 
tum-integral equation, Eq. (3). Given a pair of 
boundary-layer profiles at the same Mach Number but 
different Reynolds Numbers, an average local friction 
coefficient can be defined by 


l ine 
R, a R, I. 


It is easily shown that the average coefficient expressed 
by Eq. (7) is, for practical purposes, identical with that 


Ro, = Ro, 
R, — R 


Gi. C,dR = 2 (7) 


given by 


1 f° 
C,.. = C,dR (8) 
oe ke kin 


when a power law is assumed for the variation of C, with 
R, within the interval in question—i.e., a straight line in 
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Fic. 27. The ideal boundary layer. 
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Values of (R, — R)) calculated from Eq. (7), using 


the average coefficient from Eq. (9), are listed in Tabl 
|. As has been pointed out, the second entry for .1/ 

1.5 involves two profiles at different plate stations at the 
same tunnel pressure, with the airjet trip, and the re 
maining data involve two profiles at different pressures 
at the same station, with the fence trip. 

Consulting Fig. 25, a straight line that slopes up and 
to the right at 45° is a line on which C,(.1/),C,(O) and 
Ro(.\7) Ro(O) are equal; 
to be also a line of constant Reynolds Number. Abso 

R follow on 
asymptotic sense already defned; the result" for 1/ 
25. Values of R for the ideal supersonic 


it may be assumed initially 


integrating Eq. (2) in an 


0 


lute values of 


is shown in Fig. 
boundary laver are conveniently determined by inte 
polation, and the increment (R» — R;) thus derived for 
the various paired profiles is given in Table 1. 

Phe second method for calcul iting R 
reconciled with the first by assuming that lines of co1 
stant R in Fig. 25 have a slope which is not 45> but is at 
least independent of Reynolds Number in the range 
being considered. The problem is then reduced to one 
of determining the proper slope so that the interes pts on 
the curve .// = Oare separated by the desired increment 
in R. they corre 
spond to the values of R, and of Cy, C, CAM) /C(O 
and Cr Cr, = Re( 17) R,(0) at constant R, which are 
listed in Table 2. The result of the calculation is rela- 
tively insensitive to errors in measurement of either 
shearing stress or momentum thickness. In particular, 
it is easily shown that the value of (C;,/C;,)/(Cr Cy 
for a given pair of profiles would be unaffected if the 


Such lines are shown in the figure; 


measured values of 7,,./g or Rs were corrected by a small 
fixed percentage to account for a suspected experi- 
mental error. It follows that the development leading 
to Table 2 does not depend significantly on the question 
of uniqueness raised in the preceding sections. 

It should be pointed out that difficulty in defining an 
ideal Reynolds Number is not necessarily material in 
practical calculations of flat-plate drag. For example, 
the original values of u,x/» for a typical curve in Fig. 
27 can obviously be recovered by integrating Eq. (2), 


FRICTION ON A FLAT PLAT I 4 
» Ideal Conditioas 
) 0 OO148 0 00104 ) ) 
240 0.0 0.00 60 
MyA “Ke dRg 
—=» 10 
v Jr's, (Re 
with g( Re C,(R») defined experimentally. A cor 


responding calculation can be carried out for any plate 
flow for which the transition curve in Fig. 27 1s known or 
In practice, it may be convenient to 


of Figs. 5-14, 


can be estimated. 
begin with the coordinate system C,(R 
taking the laminar flow as an initial state and assuming 
a distribution of local friction for the first half of the 
transition region. The function R,e(R) then follows on 
integrating the conventional expression for momentum 
The result of the integration can be 


Rs) of Fig. 27 


balance, Eq. (5). 
transferred to the coordinate system C, 
and can be suitably connected to the curve describing 
the known final state of the boundary layer. Eq. (10 
used to continue the calculation in 


may then be 


definitely. 


(5) The Ideal Boundary Leyer 


The uniqueness criterion expressed by Eq. (2) 1s 
clearly a vital concept in any evaluation of experi 
mental « obtai turbulent boundary layers 
For | peed flow, this concept leads immediately to 


in intuitive definition ynolds Number and 


eventually to a quantitative definition based on certain 


characteristic properties of the mean-velocity profil 


\ccording to the evidence considered in reference 15, 


the development of a turbulent boundary layer with 
constant density may be divided roughly into four re 
gions: 

|) For momentum-thickness Reynolds Numbers Rg 
greater than about 5,000, the turbulent mean-velocity 
profile follows certain well-defined laws of similarity 
The contribution of the sublayer to the streamwise flow 
of mass or momentum is negligible. 

(2) For Reg between 2,000 and 5,000, the similarity 
laws apply, but the sublayer flow should be taken into 
account. 

(3) For Rg between 500 and 2,000, the mean-velocity 
profile is probably unique, but the similarity laws are 
not strictly valid. 

(4) For Re less than about 500, fully turbulent flow is 
rarely observed. 

In the region (1), the relationship between 7,,/q and 
R, has a simple analytic form involving four dimension- 


less parameters. If this expression is arbitrarily adopted 
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Fic. 28. Collected data at M = 2.5. 
as the uniqueness relationship of Eq. (2), then the 
Reynolds Number defined by integration” will be 
valid in an asymptotic sense for large values of Ro. 

As the generalized similarity laws for the compressi- 
ble mean-velocity profile are unknown, the correspond- 
ing definition of Reynolds Number in terms of an 
asymptotic state is not yet possible. It is therefore 
necessary to classify as tentative any quantitative re- 
sults of the present report which are not independent of, 
or at least insensitive to, the definition of an ideal 
Reynolds Number. In particular, emphasis should be 
placed on Fig. 27 rather than on Fig. 26, and more im- 
portance should be attached to the values of C;/C;, in 
Table 2 than to the values of Cp/C,r;. 

One of the original objectives of the research re- 
discriminate among several 
the compressible turbulent 


ported here was to 
theoretical analyses of 
boundary layer. However, during study of the experi- 
mental data it became obvious that the most conspicu- 
ous effect of compressibility is to move the various 
stages of boundary-layer development toward larger 
values of R». For example, the effect of surface heating 
at large Mach Numbers is known to be manifested in 
thickening of the sublayer,'* '* so that the final stage of 
turbulent boundary-layer flow must be reached at pro- 
gressively larger values of Re as M increases. This 
effect is overlooked in the theoretical literature with one 
exception.!® 

Finally, it is important to note in Fig. 27 that fully 
developed turbulent flow was apparently not observed 
in the JPL experiments for Rg smaller than about 2,000. 
In the absence of other evidence, therefore, it is recom- 
mended that values of R» substantially larger than 2,000 
be stipulated for any fundamental future study of fric- 
tion, heat transfer, or other phenomena in supersonic 
turbulent boundary layers. 


APPENDIX.—COMPARISON WITH OTHER DATA 
ari = 2:5 


For two-dimensional flow over a flat plate, it is 
readily shown that the complete form of the momen- 
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tum-integral equation (3) ist 


Tw do 
] 


6 dp 
4) / 


(2 — J\/? + 6* 1] 


q “dx q dx 


It has been found experimentally'* !**! that the value 
of the profile parameter 6*/@ is very nearly 4.2 for tur- 
bulent flow at Mach Numbers near 2.5, and it follows 
that the combination (2 — M? + 6*/6) in Eg. (11) is 
negligibly small. This property of the momentum- 
integral equation is an important reason for studying 
experimental data at M = 2.5, since the boundary layer 
should be insensitive to local variations in static pres- 
sure, 

Profile measurements in turbulent boundary layers at 
Mach Numbers between 2.4 and 2.5 have been reported 
by Wilson,'? by Monaghan and Johnson,” and by 
Rubesin, Maydew, and Varga.*!_ In each of these in- 
vestigations the momentum-thickness distribution 6(x) 
on a flat plate was observed for a fixed tunnel condition, 
and 2d6/dx can be computed immediately as a function 
of Re. Values calculated by the writer are represented 
by solid symbols in Fig. 28, and the corresponding curve 
for 7 = 0 from reference 15 is also shown. 

The data derived from profile surveys appear at first 
glance to define a unique boundary layer. Moreover, 
the fact that the experimental curves for = 0 and 
M = 2.5 are parallel suggests that C,(/)/C,(0) is, in 
fact, independent of Reynolds Number. The slope and 
ordinate of the function R,(R) are therefore affected 
equally and uniformly by compressibility, and a line of 
constant R has a slope of 45° in Fig. 28. 

Values of R obtained by interpolation in Fig. 28 are 
plotted against ™,x/»,, the Reynolds Number measured 
from the leading edge, in Fig. 29. In each case the 
ordinate and abscissa differ by a constant that agrees 
closely with the value chosen in the original paper. 
This result is not surprising, since each of the original 
analyses determined an apparent origin by referring the 
supersonic measurements to flow at JJ = 0 and by 
assuming the effect of compressibility to be nearly or 
exactly independent of Reynolds Number. 

Chapman and Kester® have recently 
directly the friction drag on the cylindrical afterbody 


measured 
of a cone-cylinder model at a Mach Number of 2.5. 
Small corrections were introduced to compensate for a 
starting length upstream of the region of measurement, 
using the fact that Cr(\/)/Cr(0) was found expert- 
mentally to be independent of Reynolds Number at 
constant Mach Number. Since the mean and local fric- 
tion coefficients satisfy a relationship C, = d(CrR)/dR, 
where by definition C, = 7,,/q and Cry = D/2zrqx, it is 
convenient to represent C; in a limited region by the 
power law Cy = AR Vm. then C; = Cr[l — (1/n)]. 
Assuming that a value of 1/m = 0.17 is appropriate to 
the data for JJ = 2.5 in Fig. 8 of reference 9, the meas- 
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urements in question lie as shown by the open points in 
Fig. 28. Because Re was not measured in these experi- 
ments, the data have been plotted against CyR/2, which 
is supposedly equivalent to Rg when the latter is suitably 
defined in cylindrical coordinates. No attempt has 
been made to correct the original values of Cy for the in- 
fluence of curvature, although it is clear that any such 
correction should affect ordinate and abscissa equally 
in Fig. 28. Values of R obtained by interpolation are 
found to agree with values given in reference 9 for a 
length ‘diameter ratio of 8 and are found to be a few 
per cent smaller for length /diameter ratios of 13 and 23. 

If any conclusion can be drawn from the experience 
of the several investigators who have measured tur- 
bulent friction in supersonic flow, it is that more study 
of experimental techniques is needed. Fig. 28 in some 
cases exhibits a serious disagreement between measure- 
ments of drag and of momentum loss at a Mach Num- 
ber of 2.5. The only direct verification of momentum 
balance known to the writer is the one reported here, 
which was obtained on a smooth but unpolished surface 
at. = 4.5 with an airjet trip. It is entirely possible 
that surface roughness may be a factor in causing the 
disagreement noted in Fig. 28. However, the writer's 
personal conviction is that too much attention, rather 
than too little, is probably paid to surface finish in 
investigations concerned primarily with flow down- 
stream of transition. 

Experience at JPL does not suggest an explanation for 
the discrepancies in Fig. 28. The effects of pressure 
gradient have already been discounted, and any sys- 
tematic error in momentum-thickness measurements 
should be common to all the experiments considered 
here. 

It should be noted in connection with profile surveys 
that the momentum thickness at large Mach Numbers is 
highly sensitive to the shape of the velocity profile at 
the outer edge of the boundary layer. On assuming 
constant stagnation temperature and static pressure, so 
that a relationship between p and u is known, the maxi- 
mum value of the integrand in Eq. (1) can be com- 
puted as a function of Mach Number. The value in 


question is 0.250 at WV = 0, 0.131 at MW = 2.55, and 
0.067 at M = 4.55. However, the integrand in Eq. (1) 
behaves like [1 — (u/u,)] near the outer edge of the 


boundary layer. On referring to Fig. 15 for a typical 
profile, it is clear that a given error in velocity must lead 
to relatively larger errors in the ordinate, and there- 
fore in 6, as the Mach Number increases. Moreover, if 
the phenomenon of intermittence observed in low-speed 
flow persists at supersonic speeds, the velocity deduced 
from impact by pressure measurements may depend on 
the nature of the intermittence and on the probe re- 
sponse to a fluctuating signal. 

It is also possible that the ambient flow in a typical 
supersonic wind tunnel deviates sufficiently far from 
the two-dimensional condition to account for all or part 
of the inconsistency observed in Fig. 28. The practice 
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of making an adjustment in the top and bottom nozzle 
contours, to account for the presence of a boundary 
layer on all four tunnel walls, leads inevitably to three- 
dimensional flow in the test section. The flow is 
slightly converging in the plan view and slightly diverg- 
ing in the elevation. The two effects cancel as far as 
the pressure distribution on the tunnel axis is con- 
cerned, yielding there a uniform value of static pressure. 
In the 20-in. tunnel at JPL, pressure and flow inclina- 
tion surveys off the nozzle centerline at a Mach Number 
of 2.6 have definitely established the presence of this 
peculiarity in the supposedly uniform flow. 
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Strength Under Combined Tension and 
Bending in the Plastic Range 


J. M. FRANKLAND* ann R. E. ROACH? 
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AS8STRACT 
A rational method is presented for the strength of rectangular 
sections under combined tension and bending in the plastic range 
The analysis is based on the accepted theory of plastic bending 
as presented by Osgood? and others. Nondimensional interac 
tion curves between axial load and bending are presented for three 
types of stress-strain characteristics representative of artificially 


aged aluminum alloys. Strength is presented in terms of allowable 


plastic strain at the extreme fiber. Some comparisons are made 


with tests on 75ST-6bars in pure bending. It is also shown how the 


change in curvature under load may be obtained, and suggestions 
ure given for the use of this to calculate strength under statically 


indeterminate conditions 


INTRODUCTION 


oe METHODS OF ANALYSIS for bending in the 
plastic range have been available for some time. 
Cozzone'! has published a practical method that has 
found wide acceptance among aeronautical engineers, 
and Osgood? has written a review of the theoretical ap- 
proach to the problem. The theory is based on the 
assumption of a linear distribution of strain over the 
depth of the section, the stresses being derived from the 
strains by use of the uniaxial stress-strain curve of the 
material. Tests, such as those by Rappleyea and East- 
man,* have confirmed the theory satisfactorily. 

Several criteria for ultimate strength in bending have 
been proposed. Cozzone' recommends that the con- 
ventional tensile strength of the material be used for the 
maximum bending stress, although it is possible to use 
his charts for any assumed value of extreme fiber stress. 
Brush, Sidebottom, and Smith‘ take the ultimate bend- 
ing moment to be that at which plastic yielding has 
penetrated halfway to the centroidal axis. Barrett° 
uses the 0.5 per cent offset yield strength for the maxi- 
mum bending stress, justifying this choice on the basis 
of British proof-load requirements. 

It should be observed that the tensile strength is not, 
in principle, a limiting stress for bending, since for duc- 
tile materials the tensile strength is merely the conven- 
tional stress at which a tensile specimen begins to neck, 
and necking does not occur in bending failures of com- 
pact specimens. 


Presented at the Structures Session, Twenty-Second Annual 
Meeting, IAS, New York, January 25-29, 1954. 

* Technical Assistant to Chief of Structures. 

+ Lead Structural Development Engineer. 


The point of view expressed in this paper derives 
mainly from experience with aluminum-alloy forged 
fittings, which in heavier sections sometimes show low 
ductility and thus have only a limited capacity for 
plastic redistribution. It has consequently seemed de 
sirable to present the strength data in terms of extreme 
fiber strain rather than in terms of stress and to use 
various values of limiting strain ranging from low to 
high ductility. Data are presented here only for the 
rectangular cross section in order to direct attention to 
the effects of ductility and the interaction of axial load 
and bending moment. 


METHOD OF ANALYSIS 


The Ramberg-Osgood® expression for the stress-strain 


curve will be used. This can be written in the form 


e= s+ (3/7)s" (1) 
where 
e —= €/€; 
S = ¢/ Cy 
oy = secant yield strength in tension, the secant 
having a slope of 0.7 
ey = oy/E 
n = a constant characteristic of the material 
o0 =} true stress 
€ = natural straint 


The first term on the right-hand side of Eq. (1) repre- 
sents the elastic component of the strain, and the second 
term represents the plastic component. 

The exponent will be taken as an odd integer, since 
this appears sufficiently accurate in view of material 
variation and it permits Eq. (1) to apply in both tension 
and compression. The assumption is made that the 
compressive stress-strain curve is symmetrical to the 
tensile branch. This is an adequate representation of 
the facts for the artificially aged alloys 145ST-6 and 
75ST-6 but is somewhat in error for work-hardened 
materials and 24S5T-3. It should be noted that Eq. (1) 
ceases to apply when unloading occurs, as unloading 
changes only the elastic component of strain. It is 

t In the strain range of interest, the difference between con- 
ventional strain, A///, and natural strain, log, [1 + (A///)], is not 
large. The reader can interpret € as conventional strain without 
appreciable error if he so prefers. 
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further assumed that the material is free of appreciable 
internal stresses. 

The undeformed cross section is a rectangle of depth 
2h and breadth 26. Rectangular coordinates x and y 
are chosen with the origin at the center, so the section is 
bounded by 


x = +6, y = +h 


The section is loaded in tension by an axial force P 
applied at the centroid and in bending by a moment ./ 
about the x-axis. The sense of the bending moment is 
such that the tensile stresses due to bending appear in 
the upper half of the section (y > 0). 

The basic assumption of the theory is that the dis- 
tribution of strain is linear over the section, and hence 


e€ = & + Cin (2) 


where 7 = vih. 

Under load, the section deforms from its original rec- 
There is a change in area and depth, 
An ap- 


tangular shape. 
and the centroid is no longer at mid-height. 
proximate treatment is adequate for the range of con- 
ditions considered in this paper. 

It will be assumed that the transverse contraction is 
half of the longitudinal extension and that the dis- 
torted section is a trapezoid. Anticlastic curvature 
and anisotropy of plastic flow are thus neglected. The 
components of displacement in the cross section are, in 
the x-direction, 


u = —(1/2)xe e- (3) 
and, in the y-direction, 
ba 
v= —s€y / e dy 
= 0 
= — (1/2)h e-[con + (1/2)eyn7] (4) 


by virtue of Eq. (2). 
Let m be the value of n at the top of the deformed sec- 
tion and m2 be the value at the bottom. The same sub- 


AERONAUTICAL 


SCIENCES-—JULY, 1954 

scripts will be used to denote extreme fiber values of ¢ 

and s. Terms containing higher powers of e; than the 

first will be dropped in the following calculations: 
From Eq. (4) one finds 


m = 1 — (1/2) elo + (1/2)q] ( 
nm = —1+ (1/2)er[c — (1/2)a]f 


The strains in the extreme fibers are 


Co + Cum 
Co + C12 


Qa = 
2 


Solving the last pair of equations for c and ¢, and 
using Eqs. (5), the following values are obtained: 


ey(€; — @2)° 
16 y (1 2 


| 
& = 5 (e; + @) + 


9 


l . 
OQ = 51 — &) tT 5 ela” — &*) (7) 


2 oe 
The centroid of the deformed section drops to ¥*, 
where 


ae 


n* = —(5/12)q & (S) 


It is convenient to characterize the strain distribu- 
tion in terms of the plastic components of the strains 
e; and e, at the extreme fibers. Let m be the maximum 
plastic strain, occurring at m, and let am be the plastic 
strain at 7p. 

Then 


—\ nm 
m == 4)s 
pe 1 t (Q) 
‘ 


( 
am = (. 


For pure bending a is approximately minus one, and, 
in the presence of axial load, it increases to the upper 
limit of plus one for pure axial load with no bending. 
Thus m is a measure of the severity of straining, and a 
gives an indication of the relative amounts of bending 
and tension. 

Now the axial load is 


P = 9 Ss dA (10) 
and the bending moment is 
M = hoy f's(n — n*) dA (11) 


the integration extending over the deformed section. 
The axial load and the bending moment will be 

measured in terms of the load Py and the moment Vy, 

which acting separately would produce a maximum 


stress of the secant yield strength—viz., 
Py = 4Abhoy 
My = (4/3)bh?ey 
Using Eq. (3), 
dA = bh(2 — e ey) dn (12) 


Whenever there is bending in addition to pure axial 
load, Eq. (2) can be solved for 7 to give 
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dn = (1/c) de 


By means of the last two equations, the integrals of 
and (11) can be written in terms of 


S sede, D= Sse? de 


evaluated by use of Eq. (1) 


Eqs. (10) 
B= fsde C= 


These integrals can be 


with the following results: 


] : e mn , 
B = — ($,7 — se*) + (s; — @So) (13) 
2 n+ 1 
l . 7+ 1 ; 
= §;3 — $o%) m (s}7 — aso") + 
C | l 5 a eS 1 QS2 
mn 
if, = 2 (14) 
2n + 
Pr et2,, 
D= — (s\4 — Sot) +m — (353 — ase*) + 
| n+3 
2n+1,, stat mn 2 
m? (s;? — ass”) + - (s} — ass) (15) 
2n + 2 an + | 
One obtains 
P/Py = (1/2c¢,)[B — (1/2) eC] (16) 
M 3 (c l D) m (“ 4 *) r 17 
= = —_ €y — 2 n (id) 
My 2c? 2 C Py 


Eqs. (16) and (17) now express the axial loading and 
the bending moment in terms of quantities that are de- 
the stress-strain curve of the material, 
ey and mn, and the extreme fiber 
The parameters c) and c,; can be 
(7), a and m are defined 


termined by 
which is defined by 
stresses and strains. 
calculated from Eqs. (6) and 


by Eqs. (9), and n* can be found from Eq. (8). 


INTERACTION CURVES FOR BENDING AND TENSION 
Nondimensional interaction charts plotting P,P, 
against .l/, .\/,; are presented in Figs. 1, 2, and 3. Each 


chart gives curves showing the combinations of axial 
load and bending for which the maximum 
tensile strain is constant—that is, for constant 
These curves of constant strain 
The actual plastic 


moment 
plastic 
values of m. are also 
curves of constant maximum stress. 
deformation is, of course, obtained from the normalized 
strain m as the product m ey. 

Each chart presents the conditions for a single ma- 
terial characterized by the strain-hardening exponent 
and 21 have been chosen as 
For 


n. Values for n of 11, 15, 
representative of artificially aged aluminum alloys. 
75ST-6, 2 can be taken while 2 = 21 is suitable 
for the 24ST-SO series and 14ST-6. For calculation of 
to deformation of the 
3 throughout. 


as 15, 
the correction terms due cross 
section, ey has been taken as 6.3 X 10 

It will be apparent that the strength is greatly affected 
by small changes in ductility if the ductility is low, say 
! or less. On the other hand, the strength be- 
comes substantially independent of ductility at larger 


for mm = 


values of extreme fiber strain, say for m > 10. In prac- 
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a function of loading for » = 11 


Fic. 5 Curvature as 


tice, these larger values of m are not often usable be 
cause of excessive deflection. For instance, for m = 10 
under pure bending, the radius of curvature of a 75ST-6 


bar would be about eight times the depth. 


PURE BENDING 
The relation between nominal bending stress and ex 
treme fiber strain can be obtained for a given material 
Nominal bend 


ing stress is defined here as the quotient of bending 


by the methods of the previous section. 


moment by section modulus. The nominal bending 


stress at failure ts the modulus of rupture. 


Curves of nominal bending stress versus the perma 


nent strain in the extreme fiber are shown in Fig.4, which 
is drawn for a material characterized by 7 15 and 
having yield strengths from 55 to 75 kips per sq.in. 


Also shown are ex 
7OST-6 


rhis is representative of 75ST-6. 
perimental data obtained from square bars of 
loaded by uniform bending moment. Because of lim 


itations of the loading jig, strains greater than 7 per 
cent could not be obtained, so that the experimental 
points at this strain were obtained on unbroken speci 
mens. Points at lower strains represent fracture and 
are moduli of rupture. 


BENDING MOMENT AND AXIAL 


LoAaD 


CURVATURE DUE TO 
From the analysis given earlier, the deformations 
accompanying various combinations of loading can be 
obtained. Of greatest interest is the curvature. 
If Ry is the radius of curvature of the centroidal axis 
in the unloaded state, and R is the radius under load, 
the difference of the extreme fiber strains is given by 


2h[(1/Ro) — (1/R)] (18) 


ey(@; — @2) = 
It is convenient to use a nondimensional curvature, x, 
which is defined as 


xk = (h/ey)[(1/Ro) — (1/R)] (19) 
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From Eq. (18), 


It should be noted that « is alw 1vs positive. 


To any pair of load values, 1/7) .\/y; and P. Py, ther 
will correspond a definite value of x. Curves of cor 
stant « are given in these coordinates for 7 1, 15 
and 21 in Fi 5,6, and 7. The curves illustrate the ip 
crease in bending flexibility at constant bending mo 
ment when the axial load is increased in the plasti 


range. 

When the relation of bending moment to axial load js 
not statically determinate, it is necessary to know th 
significant deformations in order to follow the loading 
history in bending moment and axial load coordinates 
The most important case is primary axial load a 
companied by secondary bending. 

To illustrate the use 
loading of a member that is initially curved will be 


of the x-curves, the case of axial 
considered. A rectangular bar of uniform section and 
length / is loaded by a tensile force P applied at the 
centroids of the end sections. Let x be the longitudinal 
coordinate with origin at the center of the bar, the ends 
of the bar being at +(1/2)/.. The centroidal axis under 
1 


no load is assumed to be the sinusoid given by 





Phe maximum value of y is found at mid-len; l is 
equal to a. Under load, it is assumed that the cen 
troical axis takes the shape 

\ a COS 7 
\n umption of this c! has been found suitabl 
| est ird 1 Osgo | l ia! ) uimns 
with imitial curva It uid be an t 
proximation under axicl tens! ] alues | 
a hich are not large 

The critical mis at mid-leneth, and 
Ry = s Ta 
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Fic. 6. Curvature as a function of loading for n = 15. 
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Fic. 7. Curvature as a function of loading for n = 21. 
and thus 
k = (rh/eyl*)(ay — a) (21) 
Also MV = Pa, or 
M/My = (3a/h)(P/Py) (22) 
Eliminating a between Eqs. (21) and (22), one ob- 
tains 
M 3do 3 eyxl?\ P 
= = (23) 
AW; h rh? P, 


This is the equation of a straight line through the origin 
of graphs such as Figs. 5, 6, and 7. Assuming a value of 
x, one may therefore find the associated values of 1// My 
and P/ Py for a particular problem from the intersection 
of that straight line with the appropriate « curve. 

For a specific example, take 


l= 16.2h, a = (1/3)k, e = 6.3 X 10-? 


A material for which » = 21 is assumed, so Fig. 7 ap- 
plies. 


Eq. (23) reduces to 


M My; = (1 — 0.5x)P Py 


When « is equal to 2, the parenthesis on the right- 
hand side vanishes. This is the limiting value of «x at 
which the bar becomes straight and no higher values of 


«can be obtained. 
At « = 0:2, 
M/My = 0.9P/Py 
The intersection of this line with the curve for x = 0.2 in 
Fig. 7 gives 
M/My = 0.20, 


At« = 0.4, 


M, My = 0.8P/Py 


and the intersection of this with the « curve yields 


another point 
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M/My = 0.38, P/Py = 0.47 

Repetition of this for other « values yields further points 
through which a curve can be drawn. This curve, Fig. 
8, shows the values of bending moment as a function of 


axial load while the load is increasing. 


A maximum value of J1//1\/; = 0.46 is reached at 
P/Py = 0.70, after which the bending moment de- 
creases with axial load. 
from Fig. 3 are also shown as the dashed lines in Fig. 8. 


Curves of constant m taken 


These give the maximum plastic strain in the section at 
each stage of the loading. It will be seen that a normal- 
1 is needed for elimination of 
For the 


material properties assumed, this corresponds to a 


ized plastic strain of m = 
secondary bending by plastic readjustment. 


necessary ductility of (4)(0.63) = 2.5 per cent. 


Application of this method to curved membets under 
compressive axial load is possible. However, such use 
is limited by the requirement that no significant unload- 
ing should occur in the member, and this restriction 


may rule out most of the cases of interest. 


Of course some problems will not yield to such a 
simple approach. The deflection shapes represented 
above by a single parameter family of curves might 
have to be replaced by families of two or more parame- 
ters in order to obtain a satisfactory relation between 
deflection and curvature. 
proximations by trial and liquidation of error could be 
used to satisfy the conditions at two or more stations 
along the length of the member. But no matter what 
elaborations might be found necessary in a particular 


A method of successive ap- 


case, the information provided by the « curves is funda- 
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The Transfer Function of a Rocket-Type 
Guided Missile with Consideration of Its 
Structural Elasticity 


JACK L. BEHARRELL* ann HANS R. FRIEDRICH? 


Convair-San Diego Division of General Dynamics Corporation 


ABSTRACT cated weight per pound propellant for a given payload 
and range, makes it doubtful whether the simplifica 


Assuming that the pitching or yawing motion of a rockes type 
tion, mentioned above, is still accurate enough. It js 


guided missile under the influence of thrust and air forces normal 
to its longitudinal axis can be described by the motion of an thought vitally important to consider more rigorously in 
elastic free-free beam with given mass distribution and given closed loop control investigations the elasticity of the 
stiffness distribution, the transfer function of the missile is de : 

duced. It relates the deflection angle 6 of the control agent to the 
angle 0) measured by a displacement gyro at any point x, on the 
longitudinal axis of the missile. This transfer function is to be 
used in future investigations to replace the simple transfer func- tic missile. This transfer function is to relate the de- 
i a rigid I yi » block diagr: he over: issile . ° 

tion of a rigid body in the block diagram of the overall missile flection angle of the control agent to the angle measured 


missile. 
It is the object of this paper, as a basic step in this 
direction, to formulate the transfer function of an elas 


control loop. : : 
by an angular displacement gyro at any point on the 


longitudinal axis of the missile. 
INTRODUCTION 


i» MOST THEORETICAL investigations regarding guided ae 
missile control problems, it is customary to formu- siesta 

late the transfer function of the missile itself under the The pitching and yawing motion of the missile in 
assumption of a rigid body. A selection of the funda- space under the influence of thrust and air forces nor- 
mental frequency of the closed control loop by a con- mal to its longitudinal axis can be described with good 
siderable factor smaller than the lowest natural fre- approximation by the motion of an elastic free-free 
quency of the bending modes is considered sufficient beam with given mass distribution and given stiffness 
for introducing this simplification. distribution. The following partial differential equa- 
However, the present trend to increase the overall tion is the mathematical expression for this rotary 


performance of guided missiles by reducing the fabri- motion. 


asl Oryx, 2) 1. y (x,t) Oy(x,t) O°y(x,/) _ .Ov(x,t) 
Dx? Eits)— + g(x) + m(x) + F(x)— 
ie" 


Ox? [Ov(x,t)]/Of| OF or? Ox 
0} Ov(x,/) se 
D(x) = Jé(t)A(x — x, l 
O/fL Ox 
The boundary conditions to be fulfilled for a free-free beam are 
O*v(0,t)  O%y(0,F) O*y(x,,/) O* y(x,,/) 9 ‘ 
Ox? Ox* Ox? ox? 
In Eqs. (1) and (2) the following nomenclature ap- y = deflection of any point x normal to the 
plies: desired position of the longitudinal 
s = distance of any arbitrary point on the axis of the missile 
longitudinal sie of the missile meas- EI = local stiffness distribution along the 
ured irom the Up of the nose longitudinal axis of the missile 
x) = length of the missile Pye 3 ' 
; : _ " ' —. 
t < tesbebeeene time g local distribution of the structura 
damping coefficient along the longi- 
Received June 23, 1953. Revised and received February 19 tudinal axis of the missile 
1954 
m = local mass distribution along the long 


* Dynamics Engineer. 
+t Design Specialist—Dynamics Group. tudinal axis of the missile 
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F = local coefficient of the restoring air rigid missile having the same aerodynamic shape and 
forces normal to the longitudinal axis the same mass distribution 


of the missile 


— , : ' Oy(x,/ Ove (x,t 
D = local coefficient of the damping air F(x) — F(x) — { 
forces normal to longitudinal axis On Ox 
of the missile Pa) ac , ”a) [ Ove v,t . 
Té = control force normal to longitudinal D(x) afl Ox D(x a On ‘ 
axis of missile 
5 = deflection of the control agent For the rotary motion of the equivalent rigid body, 
A(x — x,;) = dirac A function such that for the following equation holds 
> Dike sa ae) a 
SH Ae oe Xx 0 yalx,t , os 6 
t= fs A(x — x,)dx = 1 
0 where 
For the sake of simplification, the following permis X = location of the center of gravity 
sible assumptions will be introduced: @ = rigid body rotation about cg 


(1) The ‘tural damping is small enoug be : . i 
1) The structural damping is small enough to be From Eqs. (4) through (6) follows 


neglected : 
Ov(x,/ : “ 


g(x) = 0 (3) F(x) 
Ox 


normal to the longitudinal axis of the missile are de- D(x) of ome | D(x d6\t (S 
t Ox dt 


(2) The restoring air forces and damping air forces 


pendent only on the rotary motion of an equivalent 


With the first assumption mentioned above and utilizing Eqs. (1) and (6) the rotary rigid body motion is de 


scribed by the following differential equation: 


d*6(t d6(t - —_ ‘ 
m(x)(x — Xe) — + D(x) + F(x)0(t) = To(t)A(x — 4 9 
dt? dt 
Multiplication of Eq. (9) with (x — x,,), application of the Laplace transformation with regard to time /, and in 


tegration over the length of the missile results in the following equation: 


> . 7 


\ x x 
vs) f m(x)(X — Xeg)*dx + sO(s / D(x) (x — Xy)dx + (ss / Fi(x)(x — Xe)dx T (x1 — Xeg)8(s (10) 
JV J0 J 
my 
[ = / m(x) (xX — Xe)? dx (11) 
J0 


represents the moment of inertia about an axis through the cg normal to the longitudinal axis of the missile. 


where 


Introducing for the sake of abbreviation 


1 f° 
ly = D(x)(x — Xeg)dx (12 
( f V)AN Veg AX 
l PM 
jin = tJ, F(x)(x — Xeg)dx (13 
Die 
Mor = 7 T(xX1 — Xeq) (14 


the following equation is obtained from Eq. (10 


O(s) = } [uco/(s? + dos + uw) ] 5(s) | (15) 


and from Eqs. (6) and (15) follows 


VRIX,S) = t [Meo ‘(s? + dos + puw)] (x - Xag)5(s) | (16) 
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From Eqs. (1), (2), (3), (7), and (8) follows for the rotary motion of the elastic missile 


ot O° y(x,f) O*y(x,t) pe : d6(t) _ 
.| E(x) : + m(x) —— = T6(t)A(x — x,) — F(x)a(t) — D(x) (17 
ax" Ox" of? ; dt 
with the boundary conditions 
0*v(0,t)  O*V(0,t) O*y(x,t) O*%y(x,f) 
= = = = . = = = = 0 (18) 


Ox? ox? Ox? Ox? 


Application of the Laplace transformation with respect to time transforms the partial differential Eqs. (17) and 
(18) into the following set of ordinary differential equations: 


a} ... d*y(x,5) pas : 
- EI (x) : + SMX) V(X,S) = T6(s)A(x — x,) — [sD(x) + F(x) ]0(s) (19 
dx? dx? 
with the boundary conditions 
d?y(O,s) = d*y(0,s) d*y(x,s) d*®y(x),5) 
: =—“—— = 0; w= = 0 (20 
dx? dx? dx? dx? 
Inserting Eq. (15) into Eq. (19) gives 
et Oe - Meo ; 
(1 *) + s?*m(x)yv = IT a(x — *) - = [sD(x) + F(x)]$ 6(s) (21) 
dx? dx? | s? + dos + pn f 


As a first step in the direction to find the general solution of this inhomogeneous differential equation, solutions 
of the following homogeneous differential equation will be obtained. 


d? ; ad? Vu - 
(e1 ——] + s*m(x)vH = O (22 
dx? as 
with the boundary conditions 
d?y4(0,8) d*y4(0,5) 0 d*yy(X),5) d*yy(X7,8) 0 9 
- = © = e S = s = (Z ) 
dx? dx? dx? dx 
Since E/(x) and m(x) are positive functions, and since Eq. (22) is self-adjoint, the system, on setting s*> = —\ 


has a complete set of real eigenvalues AK and corresponding orthogonal eigenfunctions oO -(x,AK) reference 1, pages 
I 5 5 5 K 5 
62, 65. 


These eigenvalues are related to the natural frequencies through 


Sk= WV ‘“—Ne = +tiwx (24 


The eigenfunctions «x(x, + twx) may then be identified with the bending mode shapes. 
The general solution of the inhomogeneous differential equation Eq. (21) may be written 


y(x,s) = A(s) + Ao(s) (x — xy) + YS Ax(s\ox(x, +twx) (25 
K=1 


Combining Eqs. (21), (22), (24), and (25) results in the following three equations: 


s?m(x)A(s) = TA —-x)--— se [sD(x) + F(x) |$a(s) (26 
| s> + dos + pn f 
s?m(x)(x% — Xeg)Ao(s) = J TA(x — x) —- — [sD(x) + F(x)]a(s) 27 
: | s*? + dos + po f 
>» Ax(s) (s? + wx?)m(x)ox(x, +twx) = {TA(x — x) - = [sD(x) + F(x) ] ¢6(s (2 
Fea | s*> + dos + py f 


Integrating Eq. (26) over the length of the missile results in 


A(s) = | lt ~ = (sD + F) js 29) 
tie s°* M s? + dos + wor : : . 





wher 
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Si 
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8, me 


sk pi 


or a) 


Fro1 


(17 


(18 


7) and 


(19) 


(20 


(21) 


“tions 


(23) 


= =) 


pages 


(24) 


(26 


7 


(9Q) 
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where 


x1 
M = f m(x)dx = total mass of missile 
0 


Pui 
D = y| D(x)dx = resultant damping air force normal to longitudinal axis of missile 
0 
x! 
F = [ F(x)dx = resultant restoring air force normal to longitudinal axis of missile 
J/J0 


Eq. (29) characterizes the rigid body translatory motion of the missile. 
Multiplying Eq. (27) by (x — x,,) and integrating over the length of the missile results in 


“_ ‘ 
° ie Meo 
Ay(s)s? m(x) (x — Xq)*dx = J T(x; — Xg) -— = E D(x) (x — Xeq)dx + 
0 ( s* + dos TT Mro 0 
ae -~ . } ’ as 
F(x)(x — Xe)dx 7 6(s) (30) 
0 sa { 


Ao(s) = $[ueo/(s? + dos + pw) ]8(s)! 31) 
i 


Utilizing Eqs. (11), (12), (13), and (14) from Eq. (30), there follows 


Comparing Eq. (31) with Eq. (15) reveals that Ao(s) = @(s) describes the rotary rigid body motion of the missile. 
Multiplying Eq. (28) by ¢)(x, +iw,) and integrating over the length of the missile results in 


oo 
> Ax(s)(s? + ox’) | m(x)o; (x, +tw;)ox(xX, +1wxK)dx = | To ,(2, +iw;) — 
K=1 0 


* 


x] 7x1 
Keo . . . 
i E D(x)o(x, +twj) dx + | F(x)oj(x, +1w;)dx 
s* + dos + py 0 0 


Since the eigenfunctions ¢x(x, +iwx) are orthogonal 


ee 
~——_ 
2) 
ww“ 
ww 
to 


+ 


Tox(x)) — re ‘a ra _|s | D(x)ox dx + | F(x)or ax} 


Ax(s) = 6(s) (33) 


= 
(s> + ox’) f m(x)ox*dx 
0 
» 


From Eqs. (25), (29), (31), and (33), then follows 


|, l <0 - 
W%,s) = | 7 — : (sD + r)| +-— - (¥ — Le) + 
. - s> + dos + wn s* + dos + pry 


= OK ) (Ss 34) 


_ $ 
K=1 (s? + wr?) | m(x)oK*dx 
0 


With the help of Eq. (34) the desired transfer function, relating the control agent deflection angle 6 to the angle 


ie ; Kei he an 
Tox(x), +tlwx) — -: | | D(x)oxdx + | Posoxds 
7 - = dos TT Mri 0 0 


§, measured by any position gyro, can be found. 
Assuming that the gyro is located at x = x, on the longitudinal axis of the missile, the gyro measures the local 
slope of the function y(x,/). Thus 
) - 
6,(t) = y (x,t) (35) 
Ox ra 
or after application of the Laplace transformation 
d | 
6,(s) = ye y(x,5)} 
x Ix=x 


36) 


From Eqs. (34) and (36) finally follows 
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ie —— 1ox 1 
/ [sD(x) + F(x) loxdlx a , ToK(x _ = 1WK ) pa 
J0 


4, | Meo dx |. dx 
_ oe 7 + 37 
0 s? + dos + py xl “XI 
K=1 (5s? + ax?) m(x)dK7dx K=1 (52 + wr? m(xX dK 7dx 
i 0 Jv x 


Recognizing in ; “xy 
6 gs Fox = F(x) oxdx 1() 
oa J 0 
Mex = / m(x)bK dx (38) ; ; is 
Jo the generalized restoring air force for the kth mode and 
‘ ; in 
the generalized mass for the kth mode and in 
7x 7 _— ° 
i? co 7 K(X, + 1WK l 
Dex = / D(x)oxdx (39) - PRIX: ald i 
_ the generalized control force for the kth mode, then Eq. 
the generalized air damping for the kth mode and in (37) can be rewritten as 
6, im J Meo [ = > SDex + Fer : =| > lr dox\ (49 
6 ts? + dos + pn K=1(S* + wr?) Mex dx Ku1(S°+ wr?)MoK dx f. 7 . 


For technical applications, it will be sufficient in gen- investigate the influence of the location of the angular 
eral to consider only the first or the first and second displacement gyro and /or rate gyro on the longitudinal 
mode. stability quality of the closed control loop in order to 

find such a location that the influence is mini- 
CONCLUSIONS mized. 


From Eq. (37) it can be seen that the representation 
of the transfer function of an elastic rocket type guided REFERENCES 
missile in a general control block diagram consists of a 
basic block representing the rigid body motion of the 
missile and parallel to it one additional block for each 


' Collatz, L., Eigenwert Probleme Und Ihre Numerische B 
handlung, Akademische Verlagsgesellschaft, Leipzig, pp. 59-65, 
. = ; 1945 
vibration mode considered. : : * Kamke, E., Differentialgleichungen: Loesungsmethoden Und 
With such a block diagram describing the pitching — Zyesungen. I. Akademische Verlagsgesellschaft, Leipzig, p 


or yawing motion of the missile, it is now possible to 537, 1944 
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Effect of Large Temperature Changes 
Including Viscous Heating) upon Laminar 
Boundary Layers with Variable Free- 
Stream Velocity’ 


SOLOMON LEVY? 
University of California at Berkeley 


SUMMARY 


This is a theoretical investigation concerned with the deter 
mination of skin friction and heat transfer for a laminar boundary- 
The effects of 


rictional heating and the variation of properties are considered 


layer flow with variable free-stream velocity 


for a constant wall temperature and a constant product of vis 
cosity and density 

Characteristics of the flow are computed from the laminar 
boundary-layer partial differential equations, which are first 
transformed into ordinary differential equations. Numerical 
solutions to these equations were calculated on the differential 
inalyzer for low Mach Number flows (Mach Number approaching 
0) with a Prandtl Number of 0.7 and 1.0 and for the entire Mach 
Number range with a Prandtl Number of 1. Approximate solu- 
tions are presented for high Mach Number flows (Mach Number 
larger than 5) with a Prandtl Number of 0.7. All the computa 
tions deal with a power variation of the free-stream Mach Num 
ber, MW, = P(2s)?/, and with a permeable and an impermeable 
surface either insulated or maintained at constant temperature. 
A calculation method is presented for the determination of heat 
transfer and skin friction for an isothermal surface and an arbi 
trary free-stream velocity variation, utilizing the results derived 
from the numerical analysis of the differential equations 

Comparison of the various solutions, as well as their application 
toa specific flow example, leads to the following conclusions: 

1) The temperature difference between the wall and the free 
stream exerts a substantial effect on the skin-friction coefficient 
The heat-transfer affected radically. The 
separation point moves to regions of greater deceleration as the 


coefficient is not 


body surface is cooled. 

2) For a constant product of viscosity and density the skin 
friction increases in an accelerating flow and decreases in a de- 
celerating flow when the dissipation term is included in the 
energy equation. The same trend prevails for the heat-transfer 
coefficient but on a smaller scale 

3) For accelerating flows the skin friction increases as the 
Prandtl Number is reduced from 1 to 0.7. With extremely small 
or large Mach Number values the increase in the skin friction is 
of the order of 1 or 10 per cent, respectively, and the solutions for 
t Prandtl Number of 1 may be utilized to yield an approximate 
answer if the Prandtl Number is not too different from 0.7. On 
the other hand, the changes in the heat-transfer coefficient are 
excessive for the case of dissipative flow, and multiplication of the 

Received July 21, 1953 Revised and received February 19, 
1954 

* The author wishes to acknowledge his indebtedness to Dr 
R. A. Seban for his guidance throughout the conduct of the in- 
vestigation. The author is indebted especially to J. Killeen for 
performing all the calculations on the differential analyzer 

t Presently, Engineer at Knolls Atomic Power Laboratory, 
General Electric Company, Schenectady, N.Y. 


results for a Prandtl Number of 1 by the factor (Prandtl Num 
ber) /* does not lead to the correct answer with the exception of 


low Mach Number flows (nondissipative flows 


NOMENCLATURI 


ts = specific heat at constant pressure, B.t.u./lb. °F 

D- = significant dimension of the body, ft 

F = function appearing in Eq. (42 

f = function of the independent variable 7 

g = temperature ratio, T;/7), as 

h = heat-transfer coefficient, B.t.u./hr. ft.2/°F 
= constant equivalent to | lol ile 

k = thermal conductivity, B.t.u./hr. ft.2 °F. /ft 

VW, = free-stream Mach Number 

m = exponent appearing in Eq. (50) 

n = constant appearing in Eq. (50 

p = pressure in the boundary-layer, Ibs. per sq_ft 

q = local heat rate, B.t.u./hr. ft.? 

R= constant appearing in perfect gas relation 

r = recovery factor 

5 = function of x defined by ds = wi pi u dx 

7 = temperature in the boundary-layer, °F 

Tow = temperature of insulated surface, °F 

7 = stagnation temperature, °F 

il = velocity component in the x direction, ft. per sec 

7] = nondimensional velocity, “;/uMa 

i = velocity perpendicular to body surface, ft. per sec 

i = suction velocity, ft. per sec 

x = distance along the body surface, ft 

4 = nondimensional distance, x/D 

y = distance normal to the body surface, ft 

s* = nondimensional displacement thickness 

: = nondimensional momentum thickness 

8 = velocity distribution parameter 

6* = displacement thickness, ft 

= independent variable defined by Eq. (5) 

7] = momentum thickness, ft 

6 = nondimensional momentum thickness, (@/D) V ueD/va 

m = absolute viscosity, Ib. sec. /ft.* 

v = kinematic viscosity, ft.?/sec 

Dy = nondimensional viscosity, »/ ve 

p = density of the flowing fluid, slugs per cu.ft 

o = Prandtl Number, nondimensional 

T = shear stress, Ibs. per sq.ft 

x = nondimensional ratio, we/mp) 

y = stream function 

Subscripts 

] = condition at the edge of the boundary layer 

w = condition at the surface of the body 

a = condition at arbitrarily selected point 
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INTRODUCTION 


gaps THE SKIN FRICTION and heat transfer in 
a laminar boundary-layer flow is important in de- 
signing high-speed aircraft. Since Prandtl first derived 
the boundary-layer equations in 1904, many investi- 
gators have tried to derive new solutions to these equa- 
tions. Their aim was to obtain the boundary-layer 
characteristics for an arbitrary body in high-speed 
flight. A complete analysis of the problem, including 
dissipation, variation of properties, free-stream ve- 
locity, and wall temperature, is difficult. For a con- 
stant free-stream velocity, complete solutions have 
been presented by Chapman and Rubesin'! and by 
Crocco.” Complexity of the boundary-layer equations 
is substantially increased with variable free-stream ve- 
locity, and the problem must be attacked under addi- 
tional simplifying assumptions. 

In this investigation solutions are presented for a non- 
uniform free-stream velocity, assuming a constant wall 
temperature and a constant product of density and 
viscosity. The solutions are computed from the par- 
tial differential equations, which are first transformed 
into ordinary differential equations, and these results 
are utilized in the integral methods to supply a cal- 
culation method for cases in which the boundary con- 
ditions are more general. Ordinary differential equa- 
tions are obtained for the following conditions: 

(1) Nondissipative flows (J; — 0) for any constant 
value of the Prandtl Number. The free-stream ve- 
locity is a power function of the distance measured 
along the body, and the equations specify the effects of 
the variation of the properties. 

(2) Dissipative flows over the entire Mach Number 
range for a Prandtl Number of unity and for a power 
function of the free-stream Mach Number, J, = 
P(2s)?”*, 

(3) Dissipative flows at high Mach Numbers (.\/,; > 
5) for any constant value of the Prandtl Number and 
for M, = P(2s)**. Even though solutions over this 
Mach Number range have little practical value be- 
cause of shock-wave interaction and large pressure 
gradient across the boundary layer the obtained results 
help to specify the changes in the dissipation flow char- 
acteristics as the Prandtl Number is reduced from | 
to 0.7, typical of air. 

Previously, exact numerical solutions were obtained 
by Brown and Donoughe* only for the nondissipative 
flows of item (1). The other two types of flow have 
not been investigated. The purpose of this work is 
to recompute the nondissipative solutions of Brown 
and Donoughe on a differential analyzer for the simpli- 
fied case of pu = constant; to obtain numerical solu- 
tions for the flows of items (2) and (3); to extend the 
results for all three types of flow to the more general 
case of flow about an arbitrary body; and to determine 
the effects of dissipation and variation of properties 
through a specific numerical example. 
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PART I—SOLUTIONS OF THE DIFFERENTIAL EQUATIONS 


Basic Equations 

The boundary-layer forms of the continuity, mo- 
mentum, and energy equations for two-dimensional 
steady flow of a gas with constant Prandtl Number 
and specific heat are 


) a ‘ 
(pl) (pv) = ) 
Ox _ Oy , i. 
( Ou 4 = Op 4 O ( ~) : 
u v =— 
4 Ox oy Ox Oy * oy “) 


( or, =) ef =) 
u v = 
OM Ox dy dy \* ay 


(=) 22) a 
o Oy ‘ Oy a 


where 7’, is the local stagnation temperature defined as 
T, = T + (u?/2cp) 

For 0p/Oy = 0, the equation of state, p = Rp7T, spec- 
ifies that the density p at any x position is inversely 
proportional to the temperature 7°. 

The variation of the free-stream pressure p and free- 
stream stagnation temperature 7’, is obtained from 
Eqs. (2) and (3) for large values of y, 


—Op/Ox = pitti(du,/dx) 
with 


pili (dT, /dx) = O 


(Ou/OV)y.+ «0 = (O°u/Oy"),.. = 0 


The continuity equation is satisfied if the stream 
function y is defined so that 


pu = Oy Ov, pv = —Oy/Ox (4) 


In reference 5 Illingworth transformed the momen- 
tum and energy equations without the dissipation term 
into ordinary differential equations. The same trans- 
formation is applied here to Eqs. (2) and (3), including 
the dissipation term. 

The variables x and y are replaced by the variables 


s and n, 


z U1p1 77; 
s= Mipil, dx, 7 = a? —dy (5) 
0 (2s) o I 


and the change of variables is governed by the rela- 


tions 
O _ pty O Os O 
Ov ~ (2s) On Oy Os 
Oo (° On 4 °) 
= 1 
ae Non ds Os 


The velocity and the stagnation temperature are 
assumed to be defined by 
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EFFECT OF LARGE TEMPERATURE 
u = u;f'(n) (6) 
T, = T;, g(n) (7) 


where primes indicate derivatives with respect to the 
independent variable 7. The form of Eqs. (6) and (7) 
presumes that Eqs. (2) and (3) can be transformed into 
ordinary differential equations. The end results of 
the transformation validate this assumption, but they 
do not show by any means that the obtained equations 
are unique.t 

For » = 0, Eq. (7) specifies that the wall temper- 
ature 7), is equal to 7’,g(0), and the solutions are re- 
stricted to a constant wall temperature, since the free- 
stream stagnation temperature previously was found 
to be invariant with x. 

Integration of Eq. (6) with respect to 7 gives 


y = (2s) 


The terms of Eqs. (2) and 


‘f{(m) (S) 


(3) can now be evaluated 
in terms of s and 7 from Eqs. (5), (6), (7), and (8), and 


the momentum and energy equations become 


+ Krzywoblocki® recently investigated the general problem of 


uniqueness of the boundary-layer solutions and showed that the 
Jacobian of the transformation 


x(S), WS, 7) 
J . 
Ss, 7 


vanishes identically and that the transformation is not one-to- 


one, indicating the corresponding solutions may not be unique. 
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) os 4 S8@hf.. T 
(xf) +f - f‘> — Tr = 0 (9) 
1 


On u, ds 
7 (x ) 
7 = 0 (10) 
7 ,, On 


For the purpose of simplification x is assumed equal 


oT 


l 
( x! — fo’ ( — 1 ) 
xg ) . = on T; 


On 


with x = up up. 


to 1. Chapman and Rubesin' have shown that solu- 
tions of the boundary-layer equations for flat-plate 
flow with x = | agree well with the more exact solu 


an appropriate constant is used in the skin 
They took up = C 


tions if 
friction and heat-transfer relations. 
with the constant C being evaluated from the wall tem- 
perature condition and being incorporated in the skin- 
friction and heat-transfer expressions. 


The ratio 7'/ 7, is computed from Eq. ( 
ues ; 
= g(n) + M,? [g -—f” 11) 


where 
2] M,? = i," 


2,1, 


Substituting for 7/7, the momentum and energy 
equations become 





TABLE 1 - SOLUTIONS OF THE DIFFERENTIAL EQUATIONS 














Investigator Equation Velocity & Property Prandtl Number Thermal Condi- Type of Transformation 
Solved Mach Number Variation tion at Wall Sclution 
tartree'”) Momentum u)= px 4p Constant Any constant Variable Exact Y bed < —¥ 
Mj)z 0 7 
tokert ‘®) Energy Constant 
Temperature 
(9) poet?) 
Schuh Levy) Energy Variation of 
Tifford-Chu Power Type 
Sckort=Drewite 1”) Energy ’ Insulated 
trowa-Densughe*! Momentum 923 ~~ . 
Energy Variable 0.7 Constant 
Temperature 
Present Solutions Momentum v PPpeconstant 0.7, 1.0 7 tly" +4 
— sal ppudz 
croceo (2) Uyeconstant Variable 0.725, 1.0 
My constant 
Chapman-Rube sin‘ +) ppsconstent Any constant Variation of 
Power Type 
” . “/EFPi(’ Ley 
¥ 
Stewartson “omentum w= P( 2s) Hfeconstant 1.0 Insulated Fh Ter, Sa 
| x=[ ae ds 
Present Solutions Momentum ppeconstant 1.0 Constant 
Enerry Temperature T 
Levy-Seban‘!4) =e p/p ies 0.7, 1.0 Approximate Gl b 94 
1=0 
Present Solutions Momentum Constant Om, ad. 
Enerry = P(2s * pps constant Any constant Temperature & t 
Insulated 
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,, 25 dy y¥-!1 ue For all three cases the momentum equation is Eq. 
a FT — ~ nae , 4 
u, ds : 9 = x (14). The corresponding energy equations are, respec- 
(f’? —g) =0 (12) “Vvely, 
I T\ 7-1 g'’ + ofg' = 0 (16) 
(the +(— - \) = (1 + 7 Mi')e - Das 1 
a i = g + fg = () (17) 
‘fg ] y=] c . 
: : ( — | ( ust) -f'* =0 (13) g’’ + afg’ — (1 — a) (d*/dn*)f’? = 0 (18) 
i, 1 / 2 dn? = 
’ i aah The associated boundary conditions are, for all three 
The terms 7)\/7), and : 
cases, 
2s du l ; 
(1 + L M,? 7» = UV, f = Oorf =const., f’ =0 g=1-A 
u, ds 2 fr . 
9-> @, = |, g=1 


appearing in Eqs. (12) and (13) are written : 
which correspond to 


ie _— 1] P >. or 
— = |] (1 + . = MM; y = 0, v = Oorv = w(x), u e i, (1 — A)T,, 
T,, | = yo &, u=m, T, = T 
= 5) 2 ’ ie i 
(1 rs Y un) 2s duy aoe dM, _ / Exact solutions of Eqs. (14), (16), and (17) have been 
- u, ds M, ds computed on the differential analyzer for various values 


; i vn ‘ of 8 and a large range of A values. On the other hand, 
with the latter expression obtained from the differentia- 5 


tion of [(y — 1)/2].W\? = u,?/2c,7) and from the evalu- 
ation of d7’, ds from the statement that the free-stream 


the system of Eqs. (14) and (1S) was solved by an ap- 
proximate method. Before proceeding to discussion of 


ane : . the numerical results, a review of available solutions is 
stagnation temperature is invariant with x. 





a : ; shown in Table 1. The present types of solutions are 
| The final forms of the momentum and energy equa- 5. aided for comparison purposes. Examination of 
uous are Table 1 reveals that the case of 14, ~ 0 and o = 0.7 

f''’ + ff’ + Big — f’*) =0 (14) has been investigated previously. The other two types 
of numerical solutions presented here have not been 
fd hem 0 (y — 1)/2]40? considered by other investigators. 

1+ [(y — 1)/2|M? 

a - 
f*+=@Q (15) 
dn? 


Eqs. (14) and (15) must be ordinary differential equa- 
tions if the definitions of u and 7, in terms of a single 
variable n are to be valid. This occurs for 8 = constant 
and for any of the following conditions: 

(1) \4, — O or nondissipative flow. The Prandtl 
Number can assume any constant value, and the free- 
stream velocity satisfies the relation 

uy = P(2s)** 

(2) The Prandtl Number is equal to 1 and JJ, = 
P(2s)°*. No restriction is imposed upon the magni- 
tude of the Mach Number, and the dissipation term is 
included. 

(3) M, ~ » and M, = P(2s)? *. The effects of 
dissipation are again included, but the flow is limited to 
high Mach Numbers. Any value of Prandtl Number 
can be used. The boundary-layer equations are no 
longer valid at high Mach Numbers, and the usefulness 
of the results as J, ~ © is limited. The equations 





g(O)=1.5 | | 











presented here can be utilized, however, to specify the fe) | J 

characteristics of a dissipative boundary layer with a 02 0 0.2 °° 0.6 08 19 

Prandtl Number different from 1 and will be used pri- ee : : 
Fic. 9. Velocity gradients at the wall for constant 8 flows with 


marily for comparison purposes. f(0) = 0 
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Exact Solutions—Nondissipative Flows 


In the absence of frictional heating the free-stream 


Mach Number is 0 and 7/7); = g(n). 


and 8 = (2s/u,) (du,/ds). There results 


2 (8/2—B) 
u, = P(2s)?/* = Px /?-8 


This corresponds to the wedge flows considered by 
Hartree,’ the parameter 8 being identical in the present 
equation and in Hartree’s solutions. 

The momentum equation and energy equation are 
Eqs. (14) and (16), respectively, which were derived 
previously. They are identical to the equations ob- 
tained by Illingworth.’ Exact solutions to the sys- 
tem of Eqs. (14) and (16) were computed on the differ- 
ential analyzer. Calculations were carried out for a 
sufficient number of values of 8 and g(0) to describe the 
combined effect of free-stream velocity and properties 


variation. These cases are indicated in the following 


tabulation: 
B = 1.0, 0.5, 0.2, 0, —0.14 
g(0) = 0.6, 0.8, 1.0, 1.5, 2.0 
o = 0.7 
f(0) =0 
B = 1.0, 0.5, 0.2, 0 


g(0) = 0.8, 1.0, 2.0 

a = 1.0 

f(0) =0 

B = 1.0, 0.5, 0.2, 0 
g(0) = 0.8 

a = 1.0 

f(0) = —0.5 


The large number of cases considered prohibits pres- 
entation of all the results. Typical velocity profiles 
are shown in Figs. 1, 2, 3, 5, 6, and 7, and these demon- 
strate the effects of changing 6 or 7,,/7;. The changes 
in these profiles with respect to @ are illustrated in Figs. 
4 and 8 for the case of f(0) = —0.5, which corresponds 
to a permeable wall with a suction velocity of the form 

Myttif(O) Lew 


= (3s)'" = uf(0) 


U1 Pu 


] 
=o 


(19) 
2 Uy PyX 


The gradients of the velocity and temperature func- 
tions at the wall are shown for all cases in Figs. 9, 10, 11, 
and 12. The knowledge of these gradients enables the 
calculation of the skin friction and heat rate at the 
wall, 


Tw _ Kw (= = f’’(0) | Mu (20) 
Pitty” — pots” \Ov'’y29 =V2— 8B Yurpux ~ 
Iu x hx ‘(0 UX p, 


Ze oo T, ky Ry ns KvV2 —_ B Ku 


AERONAU 


The free-stream 
variation is determined from the relations ds = jy p,u, dx 
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Fic. 10. Temperature gradients at the wall for constant 8 flows 


with o = 0.7, ¢ = 1.0, and f(0) = 0 


The displacement thickness 6* and the momentum 
thickness # can be computed if they are expressed in 
their nondimensional form. 
the definitions 


Introduction of 7 into 


yields 


* were calculated from the known 


velocity and temperature profiles and are presented in 
the second part of the paper, where they are utilized to 
solve the momentum integral equation. 
15 show these results. 

All the solutions presented here show the effects of 
the variation of properties. 


Values of z; and z 


Figs. 14 and 


Examination of Figs. 9 to 
12 reveals that the skin friction increases considerably 
from the incompressible value when the surface of the 
body is heated or cooled. Conversely, the heat-transfer 
coefficient is only slightly affected. The effects of the 
variation of properties decreases until no effect is ob- 
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EFFECT OF LARGE TEMPERATURE 
0), which verifies the 
Even 


served for flat-plate flow (8 = 
classical result for this flow for up a constant. 
though the calculations were not carried out for sepa- 
ration flow, /’’(0) = O, the results for 8 = —0.14 show 
that separation occurs for larger unfavorable pressure 
gradients when the body is cooled. 

The solutions for a Prandtl Number of 0.7 for the 
simplified case of px = constant are a repetition of the 
calculations of Brown and Donoughe.* However, the 
solutions are obtained here by a different numerical 
technique and might help to clarify some of the diffi- 
culties of solution reported by Brown and Livingood,‘ 
did not approach zero as 


who noticed that /’’ and f’’’ 
» became large for the large values of 7,,/7) considered 
in reference 3. Comparison of the present results with 
those of Brown and Donoughe shows agreement within 
5 per cent for the skin friction and heat transfer. Cor- 
respondence in the values of the nondimensional dis- 
placement and momentum thickness is not too satis- 
factory as the ratio 7\,/7, becomes large, with a 
maximum deviation of 20 per cent noted for 7,,/7, = 
20. As 7,,/7T, increases, the velocity profiles obtained 
here and in reference 3 tend to approach | from above 
asn—> ©. The corresponding values of f’’(0) and g’(0) 
are defined well, but a small change in the fourth signifi- 
cant figure of these gradients leads to a 5 per cent 
change in the velocity profiles for » > 2.0. The func- 
tions f’’ and f’’’ will approach zero at large values of 
n only if the gradients at the wall can be specified within 
five to six significant figures, and these exceed the ac- 
curacy of the differential analyzer. Asymptotic solu- 
tions were substituted instead to determine the velocity 
profiles at large values of 7. 

Levy and Seban'‘ have obtained an approximate solu- 
tion of Eqs. (14) and (16) with a considerable reduction 
in the amount of calculations. The degree of accuracy 
of these approximate results can be determined by com- 
paring them with the present exact solutions. Agree- 
ment within at least 5 per cent exists between the ap- 
proximate and the exact solutions for skin friction and 
heat transfer, the discrepancy becoming larger as 7,/ 7; 
increases; the same trend is noticed for the values of the 
displacement and momentum thicknesses, the miaxi- 
mum deviation being approximately 10 per cent. 

The results presented here are exact for Wj, = 0. 
Their application to 7, > 0 is possible if the term 


[((y —1)/2]M?  d@ 
1+ [(y — 1)/2]M)? dn?’ 


Pre (22) 


is small as compared to the other terms of Eq. (15). 
For a Mach Number of 1, Eq. (22) has a maximum 
value 0.4 (1 — o)f’’? (0) at n = 0. This term becomes 
larger as 7.,/7, or 8 increases. For stagnation flow 
with o = 0.7, it is equal to 0.27 and 0.37 for T,,/T, = 
1.5 and 2.0, respectively. Since the other terms of the 
energy equation are equal to zero at 7 = 0, the neglect 
of the term of Eq. (22) is open to question unless the 
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Prandtl Number is much closer to 1 or the Mach Num- 


ber not too large (.\J,; < 0.7). 


Exact Solutions—Dissipative Flow with a Prandtl Number 

of 1 

For a Prandtl Number of | the solutions to the bound- 
ary-layer equations become much simpler. For flat 
plate flow, the Crocco-Stoddola relation specifies that 
the stagnation temperature is proportional to the ve 
locity. For the case of variable free-stream velocity, 
Stewartson'* has shown that the dissipative flow system 
can be reduced to the incompressible wedge flows if the 
surface is insulated. For no heat transfer at the wall, 
a particular solution of Eq. (17), which satisfies all the 
boundary conditions, is g(n) = 1, and the momentum 


equation becomes 
rr eet ay — ft) = 0 (23) 


Solutions to Eq. (23) have already been obtained by 
Hartree’ and can be used to compute the skin friction 
for dissipative flow over an insulated surface. 

The correspondence noted between the dissipative 
and nondissipative solutions for an insulated surface 
also exists for a surface maintained at constant tem- 
perature. When the wall constant, 
Eqs. (14) and (17) must be solved simultaneously to 
These 


temperature is 


obtain the temperature and velocity profiles. 
equations are identical to the low Mach Number Eqs. 
(14) and (16) if the Prandtl Number is replaced by | 
in Eq. (16). The solutions of the previous section for 
nondissipative flow with a Prandtl Number of | are 
applicable therefore to the dissipative flow considered 
here if the ratio 7, 7,, is substituted for 7 7\ and the 
velocity “ is replaced by ./; in the expression ™, = 
P(2s)*’*. Therefore, Figs. 3 and 6 represent typical 
velocity and stagnation temperature profiles for various 
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Fic. 11. Velocity gradients at the wall for constant 8 flows with 
o0 = 1.0 and a permeable wall with f(0) = —0.5 
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EFFECT OF LARGE TEMPERATURE 
values of @ and 7,,/7,,, and Figs. 4 and 8 show the 
changes in these profiles with respect to 8 when the 
wall becomes permeable and has a suction velocity 
—(1/2). 


defined by f(0) = The suction velocity is ob- 


tained now from 


} 


| uy Ky 
279 "" 
= Py u, dx 
{ 


Myf (0) 


ra (24) 
(2s) 


The velocity and temperature gradients at the wall 
are identical to those presented for a Prandtl Number 
of | in the previous section. The skin friction and heat- 
transfer coefficient are obtained from these gradients 


by the expressions 


Tx My 7; U1" - 
— — f’'"(0) = 
Pll” Py)” - (2s) 
" | 
f°" (O) / My 
Se (25) 
V2 
- \ Pu | My dx 
0 
/ uy dx h / My, dx 
q J0 - JV _ 
T, —T uy, ky, Rly 
| " 
, D u, dx 
g (0) Pu I 1 
= (26) 
V2 K Ku 


The local heat-transfer coefficient is defined as 
dw qu 
since, for an insulated surface with a Prandtl Number 


ee & 


For 1/, ~ 0, Eqs. (24), (25), and (26) reduce to the 


= 7),,, the adiabatic wall temperature. 


corresponding nondissipative equations. However, this 
correspondence between the dissipative and nondissi- 
pative flows does not help to specify the effects of dis- 
sipation, since the free-stream velocity distributions are 
different for the two flows. In the case of dissipative 
flow, the free-stream velocity 1; is 


(y — 1)/2]M/? 
uy = (2, Ts li dhe :) (27) 
1+ [(y — 1)/2|M;" 


Substitution of Eq. (27) into the definition dx = ds 
up\4; determines x in terms of s 


Approximate Solutions—High-Speed Dissipative Flow 


The energy equation, Eq. (15), can be written 


d* 
A il 


g’’ + ofg’ — (1 — o) 
dn 
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ds I + [(y — 1)/2])M;? 
dx = (28) 
mip V 2c, T, [(y — 1)/2]|M;? 


for 


, 
P(2s a2 (YQ) 


Once x is known as a function of s, the variation of 
the free-stream velocity “, with x is obtained from Eq. 
(27). The resulting expression is entirely different 


Px 3/s 


In spite of this difference in 


from the simple relation u; = obtained for 
the nondissipative flows. 
the free-stream velocity variation from one set of solu 
tions to the other, a qualitative description of the dis 
sipation effects is possible by referring to the variable 
s. The dissipative flow free-stream velocity approaches 
a constant as s goes to infinity, while for the nondissi 
pative flow it continually increases with s. At s equal 
to zero, (2s/u,) (du, ds) = 8 for the two flows, and at 
the origin the two curves for the free-stream velocity are 
tangent to each other, as shown in Sketch | for 8 = 8; > 
0. Curve Aa corresponds to the dissipative flow and 
curve Bb to the nondissipative one. 

For the same value of s, the skin friction is identical 
for the nondissipative flow and the dissipative flow 
Sketch 1. On the other hand, there also 
exists a nondissipative flow with 8 not constant which 


shown in 


has a free-stream velocity distribution of the type Aa. 
At the point A such a flow has a 8 value smaller than £, 
because du, ‘ds there is less than at the point B, and the 
skin friction for this variable nondissipative flow is 
lower than the value corresponding to the point B—1.e., 
the 
flows, addition of the dissipation term to the energy 


less than dissipative value. For accelerating 
equation corresponds to an increase in the skin friction. 
The reverse is true for decelerating flows, since the posi 
tions of the curves Aa and Bb are interchanged 

The number of numerical solutions obtained on the 
differential analyzer for dissipative flow with a Prandtl 
Number of 1 is relatively small. Most flows involve a 
wider variation of wall temperatures and of suction 
and blowing conditions. Approximate values can be 
obtained for all these conditions by using the method of 
Examination of Figs. 9 and 10 reveals 
another alternative, the 
affected by the Prandtl Number, while the heat-transfer 


coefficient for the two Prandtl Numbers considered can 


reference 14. 


since skin friction is not 


be deduced from one another by the use of the factor 
a *. Utilization of the results for a Prandtl Number 
of 0.7 becomes possible if proper care is taken in their 
conversion to a Prandtl Number of unity. 


l x d* = 9 
1+ [(y—1)/2] M2 dn? 


30) 
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As the free-stream Mach Number becomes larger and approaches infinity, the last term of Eq. (30) vanishes and 
the energy equation becomes 


g’’ + ofg’ — (1 — a) (d2/dn?) f’2 = 0 U8 


For a Mach Number of 5 or 10, the neglected factor amounts to '/¢ or |», of the term (1 — a) (d?/dn?)f” retained 
in Eq. (18). Then, the solutions of Eq. (18) are expected to be valid for 4, > 5 and can be computed for an in- 
sulated wall or a constant wall temperature. 

For no heat transfer at the wall Eq. (18) can be integrated to yield the temperature distribution in terms of the 
velocity function 


. » 


/ ie ; of. "fd 7 « se "dn a? 
g=- =1-(1- oa) / e So" san (/ e Ss a hie in) dn (3] 
ls, . 0 dn? 


Let » = 0, and the adiabatic wall temperature is 


figes ? —  ¢ ln *n ee sa @ 
. 5 . So e So 7 jie in) dn 
fo a 4 dn? 


The recovery factor r = (Tan — 11)/(u:7/2c,) becomes 


| | 7 fo" fan (f ¢ fo’ fan d? ,, ) 
‘m1 — (1 — ¢) 21 + ; R si i 7 
0) 4 (vy — b/212S Jo 0 dn? dn} dn 


and it increases as the Mach Number increases. As the free-stream Mach Number approaches infinity, the re- 
covery factor reaches the maximum value 


wd "ed Pe a Ng a? 
r=1-—(1-—a) I e So" san if e So? ian = te dn| dn (33 
J/J0 0 dy? f 


For flat-plate flow, 8 = 0, Eq. (33) leads to the same solution as the complete energy equation, Eq. (15) or (30), 
and the recovery factor is equal to the square root of the Prandtl Number. For 6 different from zero, the evalu- 
ation of the recovery factor requires successive approximations. The first approximation is obtained by substi- 
tuting the incompressible value of the function f in Eq. (33) and gives r = 0.798 for ¢ = 0.7 and 8 = 1.0. This 
value may vary as more approximations are carried out, but it is close enough to the incompressible values" to 
indicate that the effect of the Mach Number upon the recovery factor is relatively small. 


When the wall is not insulated, Eqs. (14) and (18) may be solved simultaneously according to the methods of 
Ringleb’ and of reference 14. Exponential functions are assumed for both the velocity and temperature profiles. 
These in turn are determined by satisfying Eqs. (14) and (18) along 7 = 0 as far as possible. The details of the 
method are given in reference 14. The gradients at the wall f’’(0) = —a and g’(0) = —KA are specified from the 
following two equations: 


108KA (2a? — B) = 24a° — 5a*[128(1 — K) — 2Be + 1] + 26a[158C1 — AK)? — 2(1 — K) — 380(1 — K)] 
(34 
: 120 ae 3 = D4 
24A® + —— (1 — o)A*a? — 5] o — 7A(1 — o) + — (1 — o) (1 — K)BIA% + 
K K 


0 ; 120 f 
A | - (1 — o) (1 — A)26? + — (1 — o)’a' + o(1 — KNB + 1911 — o) (1 — Ks" = 


K K? 
-é E | 1 — K)p — 28 
= == ) ae. | as ( on SB — Zoe + 
K Co K a) aN 


Examination of Eq. (34) reveals that it reduces to in the heat-transfer coefficient may be expected. These 
the low Mach Number equation as presented in refer- effects are demonstrated in Fig. 13 where the skin fric- 


ence 14 when the Prandtl Number is reduced to 1. tion and heat transfer are shown for 8 = 1 and various 
The similarity between the low and high Mach Number values of A. The heat-transfer coefficient is now de- 
equations for the term a indicates that the skin friction fined as h = q,/T,, — T,,. This change permits the 
for the same constant £ flow is not affected relatively as use of Eq. (26). The skin friction is still defined by 
the Mach Number becomes large. However, the form Eq. (25). 

of Eq. (35) is different entirely from the corresponding The curves for dissipative flow with a Prandtl Num- 


low Mach Number flow equation, and large changes ber of 1 are not included in Fig. 13. Comparison of 





1es and 
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‘tained 
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of the 


1€ re- 


(30), 
valu- 
‘bsti- 
This 
1 to 
Is of 
files, 


the 
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(34) 
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Figs. 9 and 13 reveals that for 8 = 1 the skin friction 
of a dissipative flow increases at most by 10 per cent 
as the Prandtl Number is reduced from 1 to 0.7. This 
result, coupled with the practical independence of the 
skin friction upon the Prandtl Number for the non- 
dissipative case, shows that the results for a Prandtl 
Number of | can be utilized to predict the skin friction 
The 
resulting error will not exceed 10 per cent for Mach 
The same 


over a body maintained at constant temperature. 


Number values between 0 and 5 or larger. 
argument is not valid for the heat-transfer coefficient 
because the curves of Figs. 10 and 13 are too different. 
The results presented in this section have led to some 
interesting conclusions, but they should be used with 
extreme caution. The approximate method utilized 
to solve the system of Eqs. (14) and (18) is subject to 
error. This error cannot be estimated, since no exact 
solutions to the same problem have been obtained. 
Further, the pressure gradient in the direction normal 
to the body cannot be neglected, and the derived solu- 
tions are not correct. They serve to demonstrate, 
however, that, as the Prandtl Number is reduced from 
1.0 to 0.7, the skin friction does not change at large 
Mach Numbers if the pressure is taken constant across 
the boundary layer. 
THE INTEGRAL EQUATIONS 


Part II—-SOLUTIONS OF 


The exact solutions presented in the previous section 
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and numerical work is necessary to solve the laminar 
Yet the flows dealt with a 
Particu- 


boundary-layer equations. 
special type of free-stream velocity variation. 
larly, the problem of flow about an arbitrary body can- 
When the free-stream velocity 
Karman-Pohlhausen in- 


not be solved exactly. 

varies arbitrarily, the von 
tegral equations are used generally to yield an approxi 
mate answer to the problem. Several solutions based 
upon the integral methods have been proposed in the 
past and are briefly reviewed in Table 2. The solutions 
to be obtained in this section are listed for comparison 
purposes. While have 
been analyzed already,** * the present methods will 


two of the flows considered 


vield different answers from those proposed by the other 
authors for all cases specified in Table 2. 
Basic Equations 


The method of solution starts from the integral mo- 


mentum equation 


duy , d ny dU 
21) pi + u° (p19) + pitt; 6 - 
dx dx dx 


( =) an 
(. ) 
, Oy/ y=0 is 


together with the boundary condition at the wall 


kK ( =) du, = 
= } (94) 
Oy 4 Ov/ Jy=0 ni dx 











show that a considerable amount of mathematical and with the definitions 
TABLE 2 - SOLUTIONS OF TEE INTEGRAL EQUATIONS 
Investirator Equation Solved Mach Number Property Prandtl Number Thermal Condi- Additional 
Variation tion at Wall Assumptions 
Pohlhausen’?7) Momentum M)= 0 Constant Any constant Variable Polynomial Profiles 
Kroujiline (8) Energy Constant 
P ., (19) , 
Grusch itz , Momentum Variable Pfeconstant 1.0 Insulated 
Low 21 Wei 20) 
‘ 
Yorduchow( 22 Momentum Enerry Constant Ratio of hydrodynamical 
Temperature to thermal boundary 
| layer is constant 
Eckert (©) N,s 0 Constant Any constant One-Parameter Profiles 
Setan( 25) Variable 
Levy-Seban‘ 14) 0.7, 1.0 Constant 
Temperature 
Bckert-Livingood 24) Variable 0.7 
Present Solutions j v pporenetent 0.7, 21.0 
(oc 
Rott 25) ( Momentum Variable 0.7 Insulated 
Rott-Crabtres 26) 1.0 
Present Solutions Momentun. 1.0 Insulated 
Prssent Solutions Momentum Energy Constant 
Temperature 
Present Solutions v Vv Any constant Constant y 








Temperature & 
Insulated 
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is p u 
o* = | (: — ) ay (38) 
0 pi Uy 
‘i u u 
6 = | E (1 — ) ay (39) 
0 pi ly uy} 


If one chooses to determine the momentum thickness 
# in terms of x, Eq. (36) becomes a first-order differen- 
tial equation for the unknown variable @. Solutions 
to this differential equation are based upon the simpli- 
fying assumption that the velocity and temperature 
profiles about an arbitrary body belong to a one- 
parameter family of profiles 


T./T,, = g(n, 8), um = f'(n, B) (40) 


where the proper free-stream velocity “, and the proper 
equivalent value of s have been substituted in the defi- 
nition of n. The functions f’ and g are also dependent 
upon the specified thermal condition at the wall. 


Introduction of Eq. (40) into Eqs. (36) to (39) yields 


i — | 1 
Bz? = (1 i sae ui) — (41) 
My 2 dx 
1 (6° (y — 1)/2]|.M" 
Te ( Pt) = F(B2;7) + 262," lO Ls is 
dx \ wy 1+ [(y — 1)/2]M: 


(42) 


where for 


> 


—&= | (1 — f’)dn and ¢ = / (1 — g)dn 
0 0 
” 9 é ae ¢ 
2 [ (O)s; — Bz;" (2 5 )] 


Equations (41) and (42), together with the known 
distribution of /, and the prescribed thermal condition 


F(gz,7) = 


at the wall, determine the development of 6 in terms of 
x. Estimates of 6p; /u yield of 82;" from Eq. (41) and 
thus determine (d/dx) (6°p;/4,) from Eq. (42). 

The local values of the skin friction and heat-transfer 
coefficient are obtained from the computed values of 
the momentum thickness #6. From Eqs. (25) and (26) 
there results, for 2s = 67 p,71)7/2,?, 


Ti! (1/2) pots? = 2f"(0)2;/ (10 p1/ 11) (43) 
h = k, [g’(0)2,/0K (44) 


where h = q,/(7T,, — T,,). 
The numerical solutions are considerably simplified 
when Eqs. (41), (42), (43), and (44) are written in non- 


dimensional form. These are 


A — 1 l 
d22 =| (1 +7 Mi?) = (45) 


V; \ dx 


1 (@° — 1)/2]M? 
My = ( ) = F(@z,7) + 262;? ly Ths sak (46) 
dx \i 1+ [(y — 1)/2]M;? 
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Tw lu,D 2f"(0)z; 
a | = — 17 
(1 2) pytt;" Va 1,0 V; : 
hD Va g’(O)z; 
= 18 


ky uD 6K 


where the nondimensional quantities are defined as 
follows: 


Uy, = U/Ug, y= yy Va 
6. fu,D 

d= j x = x/D 
D Va 


Here D is some significant dimension of the body and 
the subscript a refers to conditions at an arbitrarily 
selected point. 

The system of Eqs. (45) to (48) determines all the 
important characteristics of the flow about an arbi- 
trary body for a specified initial condition of 6. The 
initial value of @ is obtained from a limiting process at 
the stagnation point forx— 0. Fora small region near 
the stagnation point, the free-stream velocity is propor- 
tional to the distance x and the momentum thickness 
9 is obtained for ¥ approaching zero from the expression 


Is = 07 p;7u1" 23? 
6 = 3z,/v (du,/dz)x =o 19 


The system of Eqs. (45) to (48) is reduced further 


if F(@z,;7) is a linear function of 6z,? 


F(6z;7) = —mBz;? — n 


and there results 
A ,1 + [(y —1)/2] M2 °°”? 
_ n' —" “— xX 


9) /9 


de ~ m l Y ini l 9 “s ats ~ - 
iy 1+ s M,? dx (00 
J0 é 


where the constants vary with the particular type of 
flow considered and with the prescribed thermal condi- 
tion at the wall. 

The functions {”(0), g’(0), and F(@z;2), which are es- 
sential to the solution of the problem, vary with the 
specific type of flow considered. On the other hand, 
Eq. (45) shows that the form parameter 62,” is defined 
from the same equation for all the flows. Specific flow 
cases will be considered now. 


Incompressible Flow 


For this particular case the free-stream Mach Num- 
ber ./, is equal to 0, 2* = £, and %#, = 1. The solution 


has been indicated by Seban”’ and can be written 


where 
n 0.45, m= 5.12 0 < Ba? < 0.09 
n = ().45, m = 5.32, -0.06 < Bs;? < 0 
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Fic. 17. Dissipation and cooling effects upon the skin friction 
for flow over an arbitrary body with a = 1.0 
The function f”(0), g’(0), and Z;, which are necessary 


for the computation of the skin friction and heat 
transfer, are given in reference 23 for a Prandtl Num- 


ber of 0.7, typical of air. 


Nondissipative Flow with Variable Properties 


Here again, JJ, = 0 and s* = £&. 
thickness 8 is obtained from the simultaneous solution 


The momentum 


of the equations 


Bs? = (6 


d /@ 
1 = F(6z,") 
dx (. 


The functions F(@z,;?) and z; are plotted versus the 
form parameter 82,” in Fig. 14 for various values of the 
ratio 7,,/7, and for a Prandtl Number of 0.7. The 
gradients at the wall which determine the skin friction 
and heat transfer are given in Figs. 9 and 10. The 
linearity of the function F(8z,;*) is noted for specific 


D,) (dit, dx) 


cases only. These are: 


Telit 0.6 0.8 1.0 
0 < Ba? < 0.011 m=3.79 m = 4.28 m = 5.12 
0.04 < Bz,;? < O m= 4.45 m =5.10 m = 5.32 


with » = 0.45. The momentum thickness 6 is 


~— 


6° - m - m—Il oa 
= Ny My dx 
VD; 0 


Dissipative Flow with a Prandtl Number of 1 


When the wall is insulated g = 1 and 


— | — | 
3* = . M,? + (1 + L e Mit) § 
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and the problem reduces to the corresponding incom 
The functions f”(0) and z; are given in 


pressible one. 
thickness is obtained 


reference 23. The momentum 

from 

A 51+ [(y — 1)/2]M,7}"°"”""* 
=n = x 


V; i) 


v — n—2)/2 
eid Y ; 
/ i) ( a. ; AM, ) dx 
VU - 


where m and 7 are the constants specified in the section 
for incompressible flow. The form parameter is 


ail Z(: y-1 Vv ) S 
wile Dy 2 vo dx 


The derived equations are different from the solu- 
tions proposed by Rott and Crabtree” for the same type 


of flow—namely, 


0.45 b (1 ao 
2 


A = 
| 1 


When the wall is not insulated, the characteristics of 
the flow are determined by Eqs. (45) to (48). The 
functions F(8z;*) and z; are shown in Fig. 15. The 
gradients at the wall are obtained from Figs. 9 and 10. 
As evaluated from Fig. 15, the constants m and n ap- 



































pearing in Eq. (50) are, for 82,7 > 0 
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Fic. 18. Dissipation and cooling effects upon the heat transfer 
for flow over an arbitrary body with o = 1.0 








472 JOURNAL OF THE 









1.6 == T T 
f \ PRESENT METHOD | 


REFERENCE (19) 











° 
>» 0.8 —+ 
-_— 
2D)~ 
ho} ho) 
il 0.6 — 
04 
0.2;-—— 
“N 
ee 
N 
re) it 
0 0.4 0.8 1.2 1.6 2.0 2.4 
¥ 
Fic. 19. Comparison of results for dissipative flow over an 
insulated arbitrary body with 0 = 1.0 


High Mach Number Flow with a Prandtl Number of 0.7 


In the first part of this work, a limited number of 
results were presented for the case of dissipative flow 
with free-stream Mach Number values greater than 5. 
The gradients at the wall can be evaluated for any 
specific condition from Eqs. (34) and (35), and, in turn, 
these specify the velocity and temperature profiles." 
The functions F(§z,;") and z; can be computed then 
according to the method of reference 14. The solu- 
tions are defined here by the equations 


i — 1] li 
gs? = - (1 + 7 ue) 


V; dx 


_ d (& 
My ( ) = F(®8s;7) + 262;" 


dx Vy; 


since for large Mach Number 


is equal to 1. 


Applications of the Method of Solution 


The effects of dissipation and variation of properties 
upon the skin friction and heat-transfer coefficient can 
be demonstrated best by a numerical example. All 
the calculations presented here were carried out for the 
body specified in reference 19. The free-stream ve- 
locity and the conditions at the point a are taken as 
follows: 


MW = x(x + 5)/(¥ + 1)? 
304 m. per sec. = 997.4 ft. per sec. 
T, = 241.5°K. = 460°R. 


~ 
~ 
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The point a is chosen in the undisturbed flow up- 
stream of the body, and isentropic flow may be as- 
sumed to exist between the reference point a and the 


free stream. Thus 


Pa/ pr = (11/T,)?- 


The velocity and Mach Number distribution in the 
free stream are shown in Fig. 16. Calculations were 


carried out for the following cases: 


Incompressible flow, constant wall temperature, 
o = 1 

Nondissipative flow, 7;,/7, = 0.8, ¢ = 1 

Dissipative flow, 7),/7, = 0.8, ¢ = | 

Dissipative flow, insulated surface, o = | 


The results of the computation are shown in Figs. 
17, 18, and 19. Figs. 17 and 18 show the effects of 
dissipation and variation of properties. 

Fig. 19 is a comparison of the present results with 
those of references 19 and 27. Examination of Figs. 
17, 18, and 19 reveals the following facts: 

Addition of the dissipation term to the energy equa- 
tion decreases the value of the momentum thickness 
§. It increases the form parameter 2,2 for an accel- 
erating flow and decreases 82,” for a decelerating flow. 
These two effects combine to increase the skin friction 
for a positive gradient of the free-stream velocity and 
compensate for each other for a negative gradient of the 
free-stream velocity. The calculations show that the 
separation point moves closer to the forward stagnation 
point. 

Cooling reduces the skin friction and delays sepa- 
ration. The heat transfer is not affected. For an in- 
sulated surface, agreement within 10 per cent of the 
present method with those of references 19 and 27 is 
observed. 

The incompressible solutions obtained for a value of 
v, equal to | do not agree with the other results. Good 


agreement is obtained if one uses the relation 


Ke : 
6 ~ —m ~ m—I] ° 
— = ny My dx 
Vv; 0 


with 7, evaluated from the local free-stream conditions. 
This agreement may be fortuitous because of the flat- 
ness of the velocity curve and also because the amount 
of cooling considered in the numerical example was 
small. 
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APPENDIX (A)—DIFFERENTIAL ANALYZER SOLUTIONS 


Eq. (14) is a third-order differential equation, and 
Eq. (16) is of the second order. The simultaneous solu- 
tion to these equations involves a total of five integra- 
tions on the differential analyzer. The method of solu- 
tion consists of assuming the values f”(0) and g’(0) and 
proceeding to generate the functions f and g. The 
values of f”"(0) and g’(O) are found by trial and error. 
Different sets of values are assumed until one is found 
for which the boundary conditions at infinity are 
satisfied—namely, f’—> | andg—> 1. 

For each set of starting values, the functions /’, f, and 
g were recorded graphically, and inspection of these 
figures enabled the operator to decide whether or not 
the assumed values of f”(0) and g’(O) were satisfactory. 

The differential analyzer used was of the mechanical 
type with an estimated accuracy of 1 per cent per 
integration. The error of the solutions presented here 
is not expected to exceed 5 per cent for the five integra- 
tions leading to the results. In all cases the values 
of f”(0) and g’(O) were determined within three signifi- 
cant figures, and the error in these two quantities is 
not expected to exceed a fraction of 1 per cent. For 
incompressible flow the present values of f”(0) and g’(0) 
were within 0.5 per cent of the exact solutions.” ° 
Even though the gradients at the wall are defined ac- 
curately, the uncertainty in the velocity profile remains 
large as n > o if the ratio 7,,/7,, exceeds 1.5. For 
large values of 7,,/7,, the velocity profile is sensitive 
to the initial assumptions of f"(0) and g’(0), which 
makes the determination of the solution more difficult. 
Fig. 20 presents a typical case. A minor change in the 
value of g’(0) changes the velocity profile, and calcu- 
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Fic. 20. Velocity and temperature profile for 8 = 1.0, 0 = 1.0, and 
0M) = 2.0 
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lations beyond the cases A and B shown in the figure 
were not carried out because they would exceed the 


accuracy of the differential analyzer. Instead, asymp- 
totic solutions were used to compute the velocity and 
The solu- 
tions A and B coincide from 7 = 0 to 7 = 0.8, and the 


temperature profile at large values of 7. 


rest of the solutions can be determined from the asymp- 
totic equations for 7 > 0.8. 

Asymptotic solutions were utilized for 7, rs 
to 2. 


from above, which means that a reversal occurs in the 


equal 


For these cases, the function /’ approaches | 
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normal gradient of the velocity within the boundary 
layer. Such a solution was obtained by Hartree’ for the 
incompressible case with 6 < 0, but it was abandoned 
in favor of another possible solution for which the 
velocity gradient does not become negative. It is 
possible that there exists another solution of the system 
of Eqs. (14) and (16) for which the velocity gradient 
does not become negative for large values of g(0). No 
attempt was made to find whether this solution exists 
because it would mean that a discontinuity occurs in the 
values of f”(O) in terms of g(0) for a constant 6 value. 





Strength Under Combined Tension and Bending 
in the Plastic Range 


(Continued from page 453) 


It is of interest to record that the flexural rigidity, 
which reduces to /#/ under purely elastic coditions, is 
given by 

M 1 M 


(1/Ro) -— (1/R) « MW, 


CONCLUDING REMARKS 


The results that have been presented have been 
Ap- 
plications of the method to other compact symmetrical 


specialized to the case of the rectangular section. 
shapes can be made in a similar manner. For example, 
interaction charts have been calculated for heavy-wall 
round tubing. They differ little from the curves for an 
equivalent rectangle of the same area and _ section 
modulus. It would appear that the charts for the 
rectangle afford a useful approximation to the behavior 
of other symmetrical sections of compact type. 

The problem of fitting design has become more im- 
portant, as well as more difficult, with the present trend 
of air-frame design. Studies of the type reported here 


have proved useful in this field. More analytical work 


and experimental verification, however, are needed to 
amplify the somewhat specialized results presented in 
this paper. 
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Analysis of Flutter in Compressible Flow 
of a Panel on Many Supports 


JOHN M. HEDGEPETH,* BERNARD BUDIANSKY,? ann ROBERT W. LEONARD: 
Langley Aeronautical Laboratory, NACA 


ABSTRACT 


Two-dimensional linearized compressible flow theory is used, 
together with elementary beam theory, to analyze the dynamic 
stability of an infinitely long panel on equally spaced supports, 
having an air stream of arbitrary speed on one side and dead 
ir on the other. It is assumed that the flutter mode has a 
spatial periodicity of two bays. Longitudinal stresses in the 
panel are taken into account 

An exact solution for the condition of neutral stability yields 
two independent equations in the form of rapidly converging 
infinite series. The results of some numerical calculations made 
on the basis of the first few terms of one of the series are shown, 
ind states of stability and instability on either side of a flutter 
The implications of the 


boundary are carefully distinguished 


numerical results are discussed 


INTRODUCTION 
7 NHE POSSIBILITY OF FLUTTER of thin-sheet panels 
on the surface of supersonic aircraft or missiles 
The provision of adequate 


flutter 
particularly for 


is currently of concern. 
indeed be a 
controlling lightly 
loaded parts of the structure where extremely thin 


stiffness to prevent such may 


design factor, 
skins satisfy static strength requirements. 

Theoretical investigations, with somewhat puzzling 
results, have recently been performed in an attempt 
to provide rational criteria for the prevention of panel 
flutter. Miles! 
linearized supersonic flow theory to analyze the dynamic 


used two-dimensional, quasi-steady, 
stability of a two-dimensional, simply supported plate 
(1.e., a beam); his work, however, leads to the unhappy 
conclusion that all plates will flutter at supersonic 
Mach M = v2. This anomaly 
was only partially removed by Shen,*? who analyzed 
the same structure on the basis of unsteady linearized 
Shen discovered the 
existence of stable configurations at \J = v2, his 
results show stability for only a finite range of plate 
stiffnesses, with flutter occurring for stiffnesses above, 


Numbers below 


supersonic flow theory; while 


as well as below, this range. 

In the present paper, the problem of panel flutter 
is re-examined on the same theoretical basis used by 
Shen, but for a different structure—namely, an in- 
finitely long panel (beam) on equally spaced inter- 

Presented at the Aeroelasticity Session, Twenty-Second 
Annual Meeting, IAS, New York, January 25-29, 1954. 

* Head, Dynamics and Aeroelasticity Section 

+t Head, Structural Mechanics Branch 

t Aeronautical Research Scientist 


mediate supports, having an air stream of arbitrary 
speed on one side and dead air on the other (Fig. | 

The motivation for the present investigation is twofold 
studied bears closer geometrical 


first, the model 


similarity to an array of panels on actual aircraft 
surfaces than does a single isolated panel; and, second, 
the analysis leads to a relatively simple exact solution, 
the numerical predictions of which for a wide range of 
Mach Numbers may afford insight into the peculiarities 
that seem to be inherent in the problem of panel flutter 


SYMBOLS 


D = plate stiffness, Et°/12(1 — yp? 
E = Young’s modulus 
Fy, F = reaction forces at the supports 
WV = Mach Number 
V = Mach Number for stationary wavy wall 
P = buckling load per unit width, 72D /a? 
R = ratio of tension load to buckling load 
1 = middle plane tension per unit width 
l = velocity of airflow over the upper surface of the 
panel 
a = distance between adjacent supports 
ay = complex amplitude of the mth term of the Fourier 
expansion for the deflection 
( = velocity of sound in air 
f(t, n = function contained in the air force coefficients 
fo, f = complex amplitudes of the reaction forces 
g = damping coefficient, ratio of viscous bending 
moment to elastic bending moment at w 
= VvV-l1 
= reduced flutter frequency, wa / m 
k, = reduced first natural frequency, w,a/ 7 
Kx, = lift force per unit length of the panel 
m,n = integers 
Apu, Api = pressure perturbations on upper and _ lower 
surfaces, respectively, of the panel 
t = panel thickness 
w = panel deflection 
a = rectangular coordinates (see Fig. | 
x = longitudinal coordinate for stationary wavy wall 
= complex variable 
8, = quantity defined after Eq. (10 
¥ = density of plate material 
5(x = Dirac delta function; 6(x) = 0 for x # 0, 
tg (x) dx = 1 
€ = mass ratio parameter, pa/ yt 
N = V|M?-1 
wu = Poisson’s ratio 
p = air density 
T = time 
¢g(x1, ¥) = perturbation velocity potential 
Xn = air force coefficients 
w = circular flutter frequency 
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w) = first natural circular frequency of the panel, 


(x?/a?)V/ D/yt 


STATEMENT OF THE PROBLEM 


The configuration under consideration consists of 
an infinitely long, uniform plate, simply supported at 
equal intervals, with air flowing over the top and dead 
air below (Fig. 1). The plate is subjected to a middle- 
plane force in the longitudinal direction which is 
assumed to remain independent of the deflection. 
The problem is considered to be two dimensional, all 
quantities being constant in the direction perpendicular 
to the x-y plane, and the fluid is assumed to extend 
to infinity both above and below the plate. When 
the plate is stationary and undeflected, the pressure 
of the air below the plate is taken equal to the free- 
stream static pressure of the air above, so that there 
is no pressure differential across the plate. In the 
ensuing analysis, elementary beam theory is assumed 
to apply to the plate, and the behavior of the air both 
above and below the plate is described on the basis 
of unsteady linearized compressible flow theory. 

The term ‘‘flutter,’’ when used in a mathematical 
sense, generally refers to a steady, self-excited harmonic 
oscillation that can exist only for certain critical combi- 
the pertinent physical and geometrical 

In the next section, an equation defining 
these “flutter 
equation” It is 
to be emphasized, however, that the identification 
of the flutter states does not settle the question of 
whether or not a given plate in a given air stream is 
unstable. First of all, the a steady 
harmonic mode of oscillation does not by any means 


nations of 
parameters. 
critical combinations—the — so-called 


is derived for the present problem. 


existence of 


preclude the simultaneous existence of a different 
unstable mode of oscillation. Second, even if it can 
be shown that such an unstable mode does not coexist 
with a given flutter mode, it is still necessary to know 
whether this neutrally stable flutter mode becomes 
stable or unstable when the parameters of the problem 
are varied away from their critical values. These 
questions will be examined in detail after the flutter 
equation is derived and some of its numerical conse- 
quences are explored. 


DERIVATION OF THE FLUTTER EQUATION 


An oscillation of neutral stability (flutter) is specified 
by assuming a steady harmonic vertical motion for 
each point of the plate, with the various points not 
necessarily in phase with one another. In addition, 
it will be assumed that the motion has a spatial perio- 
dicity in the x-direction of two bay lengths. Then, the 
deflection w(x, 7) of the plate can be written with 
complete generality as 


w(x, 7) = Re ( . wr “Je (1) 
n=— « 
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Infinitely long, two-dimensional panel with flow above 
and dead air below 


where the Fourier coefficients a, are complex numbers 
and w is the (real) flutter frequency. 

Since the individual terms of the Fourier series do 
not satisfy the condition of zero deflection at the 
supports, it is necessary that the Fourier coefficients 
satisfy the constraining relationships 


» an = s- 


pe , ‘e (—1)"a, = O (2) 


These relationships will be introduced in the process 
of satisfying the condition of equilibrium of the plate. 


Derivation of Aerodynamic Forces 


It is convenient at 
magnitude of the pressure differential acting on the 


this stage to determine the 


plate while it executes the motion expressed by Eq. 

(1), with a flow to the right of velocity l’ (Mach 

Number .J/) above the plate and dead air below. 

Since linearized flow theory will be used, it is possible 

to consider separately the pressures arising from the 

individual terms of Eq. (1), and then superpose them. 
A general term of Eq. (1) is 


Re[a,e'"™*“e'*"] 


which, for n + O, can be written as 


al x+(wa/nr)r)) 


Re(a,e'"” 


term traveling wave of 


moving in the negative x-direction with 


Clearly, this 
amplitude | a, 
the velocity wa nz. 

If we now introduce the substitution 
(wa/nm)r, the problem of finding the pressures acting 


represents a 


m= xt 


on the traveling wave is reduced to the determination 
of the pressures acting on a stationary wavy wall of 
shape Re[a,e’""*’“] having a flow of velocity U + 
(wa/nm) above and a flow of velocity wa/nz below 
(see Fig. 2). 

The Ackeret solution*® to the wavy wall problem is 
followed in Appendix (A) in the derivation of results 
for the pressure perturbation at the upper surface of a 
stationary wavy wall due to a flow of Mach Number 
\/J,. It is found that: 
the right or left, 


for subsonic flow (17, < 1) to 


T pc” 


M,? 
Ap, = Re fh ee ll ; (3a) 
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FLUTTER 
for supersonic flow (17, > 1) to the right, 


™ pc* : M,? 7 
= Re| —in ae” | (3b) 


a Vv M,?- l 
and for supersonic flow (./,; > 1) to the /eft, 


— rpc i M;* eae 
Rel in an” | (3c) 
a V M,?--1 


For the stationary wavy wall of Fig. 2b, the pressures 


Ap. = 


Ap, = 


on the upper surface are found from Eqs. (5a)—(5c) 
with 17, = |M + (wa/nmc) ; the flow is to the right 
for [M + (wa/nxrc)] > O and to the left for [17 + 
(wa/nac)] < 0. 

The pressures below a wavy wall due to any flow are 
clearly the negative of those caused by the same flow 
above the wall. Hence, Ap,, the pressure perturbation 
below the equivalent wavy wall of Fig. 2b, is given by 
Eqs. (8a)—(3c) with the signs changed and with 1/; = 
wa/nmc ; the flow is to the right or left for wa nc 
positive or negative, respectively. 

The pressure differential across the wall acting in 
the upward direction is given by 


l, = Ap, — Ap, 


Summarizing these facts and replacing x, by x + 
(wa/nr)r, we find that the aerodynamic lift due to 
the motion w = Re[a,e’""*’ “e’*’] is given by 


J ™ pc” tuaK wr 
L,(x, T) = Re[x,(.M, R)ane Cae (4) 
a 
where, for n ~ 0, 
k [k 
xn(M,k) = ni f({( M+ ni+fli—m (5) 
n n 
with 
k = wa/ne 
and 
fin) = -i®/VB—1  (€>1 
. th i 
= e/V1—é (\ég I 
n 
= 1?/Ve — 1 (g¢< —1 
Finally, the result for the case » = 0 is found by a 


limiting process to be given by Eq. (4) with 
xX, = —2zk 


The total lift due to the assumed motion given by Eq. 
(1) then becomes, by superposition, 


T pc" = a . 
I(x, 7) = Re > wee | (6) 
I 


a 


Solution of the Equation of Equilibrium 


By virtue of the assumed spatial periodicity of the 
plate deflection, it suffices to satisfy equilibrium in 
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the interval (—0, 2a — 0). In this range, the equation 
of equilibrium may be written 


DwWrrrr — Tre + ylw,, = I(x, 7) + Fob(x) + 


Fié(x — a) (7) 


where 6(x) is the Dirac delta function and Fy and F; 


are the reactions at the supports at x = 0 and x = a, 
respectively. 
Letting 
Fy = Rel[foe'*"| 
F, = Re| fie’*” | 
and substituting Eqs. (1) and (6) into Eq. (7) gives 
ni3\* _f ae \* 1 pc* 
Ye | ) + 1("") — yw? — x 
n a a a 
fyo(x hun f,6(x — (1) (S) 


xn (M, A) | a,e""* " = 


Multiplying Eq. (S) through by e'""* “ and integrating 


from —0 to 2a — 0 gives, for all n, 


ni\* _ (nn \* 1 pc” 
| >( ) aa r\ ) — ylw? — xo Me) fa, = 
a a a 


fo + (—1)"h 


whence 


/ + 


(10) 


a, = : 
(2a)(yt)(re/a)*B, 


for B, + O, where 





By, = (nt + n?R)ki? — k® — exn( M, k) 
with 
y 
= os 


Yin. 
a, ad 5 4 
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(b) EQUIVALENT STATIONARY WAVY WALL 


Fic. 2. A traveling wave and its equivalent stationary wavy 
wall 
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Fic. 3. Modes involved in three-term solution 

R = T/P., (ratio of tension load to buckling load) 
k = wa/nc (reduced flutter frequency) 
ki = wa nc (reduced first natural frequency) 
€ = pa/ryt (mass ratio parameter) 
M = U/c (Mach Number) 


Now the restraint condition on the a,’s, Eqs. (2), 
which may be rewritten 


> a, = 0 (11) 
(Hv even) 
> a, =0 (12) 


n x 
(n odd) 


must be satisfied. Substituting Eq. (10) into these 
conditions gives 


(fo+ fi) > (1/8) = 0 


n —= 
(nm even) 


and 


(fo—-fi) > (/B,) = 0 
n x 
(n odd) 
It can now be seen that the simultaneous satisfaction 
of the equilibrium equation, Eq. (8), and the constraint 
conditions, Eqs. (11) and (12), can be achieved in a 


nontrivial fashion in several independent ways. We 
can have either 
Sfo —= fi 
la, = 0 for n odd 
> 1/6, = 0 (13) 
n —« 
(neven) 
or 
ffo = hi 
la, = 0 for n even 
> 1/4, = 0 (14) 
a= - 
(n odd) 


or, conceivably, Eqs. (13) and (14) could hold simul- 


taneously. 
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Eqs. (13) and (14) are conditions for two independent 
types of flutter. the 


terms in the Fourier expression for w 


Examination of nature of the 


reveals that 
Eq. (13) corresponds to a flutter mode that is identical] 
in each bay, whereas Eq. (14) corresponds to a motion 
that is identical from bay to bay, but with alternating 
direction. 

In the derivation of the above criteria it was assumed 
that 6, ~ 0; the implications of vanishing £8,.'s should 
be explored. It can be seen from Eq. (9) that (for 
n even, say) the vanishing of at least one 8, requires 
that fo = 
from zero are those for which the corresponding 8,'s 


—f;; then the only a,’s that can be different 
vanish. Consequently, the constraint relation, Eq. 
(11), as well as Eq. (9), can be satisfied by nonzero 
a,’s if two or more 8,'s vanish simultaneously. Similar 
considerations apply to the satisfaction of Eqs. (9 
and (12) by the vanishing of two or more 8,’s for n 
odd. 

To summarize, the criterion for flutter is Eq. (13), 
or Eq. (14), or the simultaneous vanishing of two or 
more £,'s of the same type (i.e., odd or even). 


Discussion of Flutter Equations 


The flutter criteria, Eqs. (15) and (14), are complex 
equations defining critical combinations of the five 
nondimensional parameters .1/, k, ki, «, and R. Conse- 
quently, one may expect that, when any three of these 
parameters are specified, the remaining two may be 
determined. 
the parameters k;, «, and R are known beforehand, 


For a given panel in a given air stream, 


and, in principle, it would therefore appear possible 
flutter Mach Number V/ and flutter 
In practice, however, it is much more 


to solve for a 
frequency k. 
convenient to assume values of 7, k, and R and solve 
for k; and ¢, since .V/ and & are included in the equations 
in a complicated manner as a consequence of the way 
they appear in the expressions for the aerodynamic 
forces. Indeed, it appears desirable to survey the 
entire range of critical combinations of the parameters 
of the problem, since, as has been mentioned, the 
establishment of a flutter speed for a given configura- 
tion does not constitute a conclusive determination of 
stability or instability. 

An indication of the results of such a survey is found 
in the next section, wherein numerical calculations 
are carried out for the modes of flutter corresponding 
to Eq. (13), with R = 0 and with the retention of only 


three terms of the infinite series. 


RESULTS FOR THREE-TERM SOLUTION 


If only the terms for n = —2, 0, and 2 are retained, 


Eq. (13) becomes 
(1/Bo) + (1/82) + (1/B-2) = 0 (15) 


Examination of the Fourier terms corresponding to 


a, Qs, and a—» and consideration of the constraint 
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condition, Eq. (11), show that the three-term solution 
involves two degrees of freedom corresponding, effec- 


9 


tively, to the two modes of deflection shown in Fig. 3. 


Flutter Loci at Various Mach Numbers 


Eq. (15) can be solved, with R = 0, by finding combi- 
yations of e€ and k, which satisfy the equation for 
various combinations of J and k. At the same time, 
consideration should be given to the possibility of the 
existence of flutter modes corresponding to the simul- 
taneous vanishing of two (or all three) of the 8,’s. 
The solutions will be discussed separately for subsonic 
and supersonic Mach Numbers, since they are of differ- 
ent character for these two régimes. 

Subsonic Speeds (M < 1).—It can be seen from the 
definition of x,(./, k) that for « > 0, k ¥ O, and for 
all \J < 1, the imaginary part of 6) cannot vanish and, 
furthermore, must have the same sign as the imaginary 
parts of 6: and 8 Consequently, Eq. (15) cannot 
be satisfied in the subsonic régime for any k # 0; 


for k = O, 8) vanishes, and hence Eq. (15) is no longer 
applicable. For this case of k = O, however, all three 


of the 6,’s vanish simultaneously whenever 


k,? = «M?/8V 1 — M* (16) 
For k ¥ 0, it is impossible for any two 8,,’s to vanish 
simultaneously. Hence, at all subsonic speeds, flutter 
occurs only for those combinations of ¢, ki, and A/ 
which satisfy Eq. (16). Since this flutter occurs at 
zero frequency, it actually corresponds to a condition 
of neutral sfatic stability. It may be noted that, 
since the solution of Eq. (16) is of the type that stems 
from the simultaneous vanishing of two or more 8,,s, 
it remains valid even for the case where all of the terms 
are considered; the 


of the original Fourier series 


additional terms, however, would provide additional 


solutions. 


contrast to the 


Supersonic Speeds (M > 1).—In 
situation at subsonic speeds, it can be shown that no 
two of the 8,,s can ever vanish simultaneously at super- 
sonic speeds, so that flutter modes can exist only by 
virtue of the satisfaction of Eq. (15). It turns out 
to be possible, for the presently considered three-term 
solution, to obtain from Eq. (15) explicit formulas for 
eand &, in terms of \/ and k [see Appendix (B)]. For 
each supersonic Mach Number, these formulas provide 
a parametric relationship between ¢ and k, in terms of 
k, so that curves of « versus k; can be constructed at 
each Mach Number. 

The result of such a construction, for many different 
Mach Numbers between 1 and 2.5, is the rather 
complicated and bewildering array of curves in Fig. 4, 
some of the details of which are discussed in Appendix 
(B). The situation is somewhat clarified by Fig. 5, 
which shows schematically a pattern of flutter loci 
for each of several Mach Number ranges within which 


the pattern is typical. 
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In Figs. 4 and 5, the arrows on each flutter locus 


show the direction of increasing & values. In the 
subsonic range, there is, of course, only one flutter 
locus (Fig. 5a), given by Eq. (16), along which k = 0 


For supersonic Mach Numbers, there may be as many 
as four separate curves, which have been designated 
as types A, B, C, and D. Each type results from a 
distinct range of variation of reduced frequency k, 


order corresponds to 


alphabetical 
(The numerals in each of the 


and increasing 
increasing levels of k. 


diagrams in Fig. 5 will be explained later. 

It is interesting to trace the metamorphosis of various 
types of flutter loci as Mach Number is increased. The 
type B locus, for example, starts out just above 1/ = 1, 
as shown in Fig. 5b; pulls away from the origin at 
M = WV/4/3 (Fig. 5c); 
1.24493 (Fig. 5d); “pulls through’ at some (unde 
termined) Mach Number below 2 (Fig. 5e); starts 
again at some Mach Number slightly less 
pulls through again and extends 
and, finally, at 1J = 2.00410, 
lies down The type A 
growing out from the origin gradually at IJ = y2 
(Fig. 5f), and in the vicinity of 17 = 2.05817 it distorts 
radically and whips out to an infinite length (Figs. 
5k, 51). The type C curve, remarkably, appears full 


begins to form a loop at JJ = 


looping 
than 2 (Fig. 5g); 
to infinity (Fig. 5h); 


the &)-axis. locus starts 


on 


grown at \J = 2 (Fig. 5i) and undergoes little change 
thereafter. The type D curve, on the other hand, 
starts to sprout at 1/ = 2 and grows gradually (Figs. 
5j, 5k, and 51). 

The significance of the various flutter loci must 


now be assessed. At any given Mach Number, every 
point in the k,-e plane defines a particular plate-air 
configuration, and we wish to know whether this 
configuration can have an unstable mode of oscillation. 
Further, it is to be realized that the configuration 
may possibly have more than one mode of unstable 
motion. We shall consider the number of such un 
stable modes as the degree of instability of the con 
figuration. Now, certainly, a configuration corre 
sponding to a point on one of the flutter loci possesses a 
neutrally stable mode of vibration, and it is plausible 
that, irrespective of the behavior of other possible 
modes, this mode of neutral stability becomes unstable 
on one side of the flutter locus and stable on the other 
side. Thus, at each Mach Number, the flutter loci in 
the ki-e plane separate regions of various degrees of 
These loci might therefore be more ap 
“stability boundaries.’ 


instability. 
propriately called 

It is clearly of interest to determine the degree of 
instability in each of the regions separated by the 
“safe” 

This 
determination has been made, and the numerals in 
the various regions of Fig. 5 show the degree of in- 
stability; the method of calculating the degree of 
instability is outlined in the next section. 


stability boundaries in order to discover the 
regions where the degree of instability is zero. 
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Fic. 4. Flutter loci for three-term solution with R = 0. 
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. (continued ) 
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Identification of Regions of Instability 


An unstable (divergent) mode of oscillation of the 
plate is represented by Eq. (1), with w a complex 
number having a negative imaginary part. Conse- 
quently, if the aerodynamic loading function x,(.l/, k) 
could be properly defined for k in the lower half of the 
complex k-plane, the flutter equations, Eqs. (14) and 
(14), could be interpreted as frequency equations for 
unstable modes of oscillation, as well as for neutrally 
stable modes. The required extension of the domain 
of applicability of x,(.1/, k) is indeed possible; it can 
be verified that the (unique) analytic continuation of 
xn(., k) into the lower half k-plane is defined as 


k k sa 
x(M,k) = 0 (ur + ) + e( ) (14) 
n i 


{for Im(k) < 0, n ¥ O| 


follows: 


where 
(s) = —i27/V 2? — 1 


iu 
Ss 


and the single-valued branch of the multiple-valued 


function g(z) that is chosen is the one for which the 


z-plane has a branch line on the real axis from —1 to 


l,and g(\/2) = —22. For Im(k) — 0, this definition 
coincides with the one previously given [see Eq. (5)]. 


The analytic continuation of X» for complex k is simply 


—2ik. It can be shown, rigorously, by means of 
Fourier transform 
theory, that the formula for lift, Eq. (6), remains valid, 
with the use of Eq. (17) for x,, when w has a negative 


[It is of interest to note that there 


analysis and complex variable 


imaginary part. 
does not exist an analytic continuation of the pre- 
viously defined x, into the upper half k-plane; conse- 
quently, it appears that Eqs. (13) and (14) cannot be 
used to calculate complex frequencies of stable (de- 
caying) modes of oscillation. | 

With the use of the definition (17) for x,, a solution 
of the flutter equation that gives a & in the lower half 
plane means that the configuration represented by the 
selected combination of values of 1/7, R, k;, and € is 
(It can be shown, incidentally, that the 
the 


unstable. 


possibility of an unstable mode arising from 
simultaneous vanishing of two $’s for a complex k can 
be ruled out.) Thus, in the presently considered 
three-term solution, one could, in principle, determine 
the stability of a point in the k-e plane for any Mach 
number by solving Eq. (15) for k and seeing whether 
there are any “unstable” roots (i.e., roots in lower 
half k-plane). the 
roots can be determined without their exhibition by 
use of the ‘Nyquist diagram’ concept. This concept, 


which is familiar to electrical engineers and has already 


Fortunately, existence of such 


been used in connection with wing flutter by Dugundji,* 
is explained below. 
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Consider, in general, a function F(z) that is analytic 
except for poles, for Im(z) < 0. Let N be the number 
of zeros and P the number of poles of F(z) having 
negative (Each zero of order a@ is 


imaginary parts. 


counted a@ times, and similarly for multiple poles.) 


Consider further the contour C shown in Fig. 6, con- 
sisting of an infinite straight line in the lower half plane, 
parallel to and infinitesimally close to the real axis, 


together with an infinite semicircle in the lower half 























plane. Let z traverse C once in the clockwise direction, 
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Fic. 6. Contour used in constructing the Nyquist diagram 


which 
in the F- 
Then, according to complex variable theory, 
the number of clockwise 


and trace out the corresponding path of F(z) 
we shall call here the “ 
plane. 


Nyquist diagram”’ 


encirclements of the origin 
made by the Nyquist diagram is equal to NV — P. 


is taken as (1/69) + 

1/82) + (1/B-2), and k is the independent complex 
variable. Clearly, for each Mach Number the number 
of unstable zeros of F 


In the present problem, F 


is constant within each of the 
regions boundaries in the 


The poles of F occur only at 


separated by stability 
ki-e plane (see Fig. 5). 
the zeros of the 6,'s; since none of the 8,,’s themselves 
have poles for Im(k) < 0, P can readily be found by 
constructing an auxiliary Nyquist diagram for each 
of the N, and 
hence the degree of instability in each region of the 


6,'s. It is thus possible to evaluate 
construction of 
Actually, it suffices to 
construct a complete Nyquist diagram for F in only 


ki-e plane, simply by the Nyquist 


diagrams in each region. 
one of the regions and simply evaluate the change in 
the number of clockwise encirclements when a stability 
boundary is crossed. (The Nyquist diagram for each 
configuration represented by a point on a stability 


boundary must, of course, pass through the origin. 


The results of such a study by means of Nyquist 
diagrams are summarized in Fig. 5. The integer 
within each region is the degree of instability as given 
by the number of unstable roots of Eq. (15) determined 
by the Nyquist-diagram analysis. That the degree 
of instability in the unstable region is always even is 
caused by the fact that complex k-roots always occur 
in pairs of the form +a — 1b, corresponding to pairs 
of unstable really 
can be seen from Fig. 5 that, as one proceeds along the 


stability boundaries in the direction of increasing k 


modes which are identical. It 


the side of the curve for increased stability is always 
to the right for the type A and C boundaries and to 
the left for the B and D boundaries. 
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Other Criterta for Instability-—A word of caution 
might be inserted concerning the applicability of two 
other techniques that are in current use for determining 
the stability of a given configuration. One method 
of widespread popularity in wing flutter analysis is 
to calculate the amount of artificial structural damping 
required to permit the existence of a neutrally stable 
mode; if the damping required is positive, the original 
configuration is considered to be unstable. However, 
to the authors’ knowledge, it has not been proved 
that a decrease in structural damping always leads to 
decreased stability. Indeed, an example will be given 
later where this is not necessarily the case for the 
present problem. 

Another technique that has been applied to wing 
flutter? and was used by Shen for panel flutter’ is to 
calculate the energy per cycle required to vibrate the 
structure harmonically by means of an external force; 
the structure is considered unstable if this energy is 
negative (i.e., if the external agency absorbs energy). 
Unfortunately, this concept is not always valid, as 
1939, and, according to 
Frazer 


was shown by Frazer® in 
Frazer, before him by Southwell in 1927; 
showed that the total work per cycle done by an excit 
ing force may well be negative even when the structure 
is stable. 
Results for a Particular Altitude and Material 

The results embodied in Figs. 4 and 5 may be more 
clearly visualized by considering their implication for 
a particular plate material and a particular flight 


altitude. The parameters « and k; may be written as 


e = (1/2)(p/7)(a/t) 


mt E 
€ @V12(1 — u*)y 
Thus, the product 
p E 
eX 12(1 — w)y 


ek, = 


depends only on p and c (determined by altitude) and 
on £ and y (determined by the plate material). Conse- 
quently, for altitude and material fixed, the points of 
interest in the k-e plane lie along a hyperbola. Fig. 
7 shows, in the k-e plane, the region of practical 
interest; the stability boundaries for 1f = 1.4, and 2 
are shown for comparison purposes. The bounding 
hyperbola of the cross-hatched practical region corre- 
sponds to magnesium at sea level; as altitude and ma- 
terial density are increased, hyperbolas further inside 
this region are obtained. Once altitude and material 
are specified, both « and k; depend only on the thickness 
ratio t/a; stable and unstable regions can then be 
shown in a chart of ¢/a versus Mach Number by cross- 
plotting the data of Fig. 4 and using the information 
in Fig. 5 to specify the degree of instability. This 
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has been done for the case of aluminum at sea level, 
with the results shown in Fig. 8. 

The shaded regions in Fig. 8 correspond to unstable 
conditions, and the degree of instability has been noted 
in the various regions. (There is actually an extremely 
narrow region of degree 6 in the upper shaded region 
near \f = 2 which cannot be shown with the scale 
used.) The implications of the results shown, as 
well as the significance of the dotted line and the 
sprinkling of symbols, are discussed in the next section. 


DISCUSSION 


At first sight, Fig. 8 presents a discouraging picture 
indeed and one that is even less optimistic in many 
respects than the results of Miles and Shen. No 
stability at all appears possible, for any thickness, at 
supersonic Mach Numbers below approximately 1/ = 
1.25, and it appears that, even above this Mach 
Number, a sufficient inctease in thickness can always 
render a stable panel unstable. The situation may 
be even worse than shown, for only three terms have 
been considered in the ‘“‘even’’ flutter equation, Eq. 
(13). Furthermore, the flutter analysis of the panel 
is not complete without consideration of the “odd” 
solution, Eq. (14). In the subsonic range, the highest 
static stability boundary 
solution is actually given by the dotted curve in Fig. 8 


resulting from the odd 


and corresponds to greater thicknesses than the even 
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solution. No numerical results have been obtained 
for the predictions of the odd solution at supersonic 
speeds. 

Even if we disregard this last consideration, the 
results of Fig. 8 obviously require modification or 
reinterpretation before they can be considered useful 
for design. It is easy to discount the subsonic results 
close to J = 1. At sonic velocity, linearized aero- 
dynamic theory predicts infinite pressures in steady 
flow, and, since the points on the subsonic stability 
boundary correspond to zero flutter frequency, their 
validity breaks down close to 1J = 1. A reasonable 
empirical modification of the results in the subsonic 
range would be to assume a horizontal cutoff on the 
stability boundary at about JJ = 0.95. 

Consider next the supersonic regions; the expression 
for the aerodynamic lift also becomes infinite at certain 
combinations of k and J/. These singularities in the 
air forces actually appear only because the plate is 
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infinitely long, and they correspond to a sort of resonance 
between the plate and the infinite field of flow; such 
singularities do not occur when the pressures on finite 
panels are calculated by linearized theory. It is 
significant that the flutter frequencies associated with 
most of the points on the higher stability boundaries 
are extremely close to these “‘resonant’’ frequencies; 
the infinities consequently may be partially responsible 
for the peculiarities of the solution. However, a 
simple rational procedure for removing the effects of 
these infinities is by no means as evident as it was in 
the subsonic region. Indeed, the infinities must not 
be entirely responsible for the unreasonable results, 
since they were not involved in Shen's results for a 
finite panel at J = 2. How, then, can the present 
results at supersonic speeds be 
attempts to do so, one unsuccessful, the other promising, 


rationalized? Two 


follow. 


Effect of Structural Damping 


The idea presented itself that perhaps the inclusion 
of a small amount of structural damping in the analysis 
might wipe out vast regions of instability in Fig. § 
and render the results more reasonable. Accordingly, 
the term (Dg/w;)W,::,, representing viscous (Sezawa 
damping was included in the differential equation 
of equilibrium, Eq. (7). [The addition of this type 
of damping term was chosen in preference to the more 
common use of a complex stiffness D(1 + ig), because 
the latter is strictly applicable only to steady harmonic 
motion and leads to difficulties when complex fre- 
quencies are considered; the damping forces arising 
from both procedures are identical when w = «).| 


Fig. 9 shows the stability boundaries for .1J = 1.8 
obtained with values of g equal to 0.01 and 0.03, as 
well as the original stability boundaries for g = 0. 
As can be seen, not only are the stable regions inside 
the loops not greatly increased in extent by the addition 
of damping but the stable region in one of the loops is 
actually decreased in size. 


An incidental outgrowth of this study is shown in 
Fig. 10, which, again for 1/J = 1.8 and for e = 0.004, 
shows calculated values of g required to produce 
steady plotted against k;. The curves 
pass through zero at the stability boundaries; but the 
remarkable fact to be noted is that, well within the 
second stable region, the addition of damping leads to 
a steady oscillation. Thus, the concept that removal 
of damping tends to destabilize is evidently not valid 


oscillation 


here, and, as was mentioned previously, doubt is cast 
on the use of the “artificial damping’ concept for 


determining whether a given system is stable. 


Rates of Growth of Unstable Oscillations 


A more fruitful attempt to rationalize the results of 
Fig. 8 was made by actually calculating (approximately 
the magnitude of the unstable roots at various points 
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TABLE 1 


Ranges of k Corresponding to Flutter Loci 
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Type Type 
VW A B 
l 0,0.2 
2 0.2698 , 0.4 


3 0.536 0.6 
0.7602 ,0.8 


' 


; 0, 0.59236 0.9738 , 1.0 

6 0, 0.91080 L.t7e8 12 

a 0, 1.17648 1.3846 , 1.4 

S 0, 1.4194 1.5876 , 1.6 

Q 0, 1.65000 1.79012, 1.8 

2.0 0, 1.87285 1.99184, 2.0 
2:1 0, 2.00000 
2.2 0, 2.00000 
3.3 0, 2.00000 
2.4 0, 2.00000 
a 0, 2.00000 

When |k, > 2(.M — 1), no solution exists because 


of the fact that all three 8,’s always have negative 
imaginary parts (for e > 0). All of the real solutions 
to Eq. (15) are therefore given by Eq. (B1). 


Examination of the x,’s reveals that yz and yo;x — 2; 
are even functions of k and that yo; + x_», is an odd 
function of k; hence, the values of « and k; for k < 0 
are the same as those for k > 0, and our attention can 
be confined to the range 0 < k < 2(\J — 1). Results 
for many values of \/ obtained by varying & in this 
range are shown in the form of plots of ¢ against &; in 
Fig. 4. At any one Mach Number, there are seen to 
be as many as four types of flutter loci; each locus 
arises from a different range of k, and the distinct 
ranges result from the fact that the quantity inside 
the square-root symbol in the expression for e, Eq. 
(B1), changes sign as many as three times as k is 
increased from zero to 2(.\J — 1). The ranges of k 
for which this expression is positive (and, hence, for 
which real values of k; and ¢ result) are given in Table 
1 for Mach Numbers between 1.1 and 2.5, in increments 
of 0.1. 
corresponding k)-e curves in Fig. 4 are designated as 


At each Mach Number, each range and the 


types A, B, C, or D, in order of increasing k. 


It should be noted that many of the k-ranges in 
Table | are extremely small and could be easily over- 
looked if care is not taken. In addition, the values of 
k in the ranges often are extremely close to those values 


and M. No 


Type ype 
Cc D 
2.00000, 2.00000 
2.14712, 2.18670 2.19670, 2.2 
2.32684, 2.3836 2.39664, 2.4 
2.52189, 2.5832 2.59697, 2.6 
2.72396, 2.78365 2.79726, 2.8 
2 9288 , 2.9854 2 99746, 3.0 


of k where one or another of the y,’s becomes infinite: 
these critical values are k = 2 and 2(.\/ — 1), the 
latter value coinciding with the upper limit of possible 
k’s. 

In connection with Fig. 4, it should be pointed out 
that, at some Mach Numbers not listed in Table 1, 
some curves undergoing rapid transition with change 
in Mach Number have been shown. In one instance, 
JJ = 2.058171, only a short segment of the type A 
curve is shown to illustrate the start of the rapid 
transition that occurs just below JJ = V24+ v5 = 
2.058171027. 
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On the Approximate Solution of the Laminar 
Boundary-Layer Equations 


ITIRO TANI?# 
Cornell University 


SUMMARY 


A simple method is developed for solving approximately the 
equations of the laminar boundary layer. The velocity profile is 
assumed as a member of a one-parameter family of curves, the 
parameter being different from the usual Pohlhausen parameter. 
The Pohlhausen parameter is calculated by a simple quadrature 
formula, which is derived from the energy-integral relation. The 
relation between the profile parameter and Pohlhausen parameter 
js determined from the combination of both momentum-integral 
ind energy-integral relations. The accuracy of the method is 
examined by comparing the results with those of exact solutions 
for the flow of incompressible fluid. Satisfactory agreement is 
obtained. The method of solution is then extended to the flow 
of compressible fluid along a heat-insulated wall, assuming the 
Prandtl Number slightly different from unity and the Sutherland 
formula for the variation of viscosity with temperature. As an 
example of the application, the effect of compressibility on separa- 


tion point is discussed. 
List OF SYMBOLS 
r = distance measured along the wall 


y = distance measured normal to the wall 
velocity component in the direction x 


ll 


velocity component in the direction y 


ll 


1 = absolute temperature 
= speed of sound 

p = density 

m7 = viscosity 

v = kinematic viscosity 

k = thermal conductivity 


= specific heat at constant pressure 
ratio of specific heats 
c = Prandtl number, Eq. (42) 
= 1 — Va, Eq. (46) 
boundary layer thickness 


» 


ll 


displacement thickness, Eq. (5) 


6 = momentum thickness, Eq. (6) 


energy thickness, Eq. (7) 


\ = Pohlhausen parameter, Eq. (24) 

1 = velocity profile parameter, Eqs. (11), (64) 

: oe 4 = coordinate in the transformed flow, Eq. (54) 

U,V = velocity component in the transformed flow 

A = boundary layer thickness in the transformed flow 
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On leave from the 


A* = displacement thickness in the transformed flow, 
Eq. (61) 

+s = momentum thickness in the transformed flow, Eq. 
(62) 

o” = energy thickness in the transformed flow, Eq. (63) 

A = Pohlhausen parameter in the transformed flow, 
Eq. (77) 

A = a constant, Eq. (23) 

B = a function of x, Eq. (56) 

c = a function of 7, Eq. (48) 

D, E, F = functions of a, Eqs. (17), (18), (19) 

G, H = functions of a, Eq. (13) 

P,@ = functions of a, Eqs. (20), (21) 

g J = functions of a, Eqs. (73), (74) 

S = a constant, Eq. (48) 

? = a function of u/u;, Eq. (50) 

Re = Reynolds Number, 4;x/», 

M = Mach Number, u/c; 

Cy = coefficient of skin friction, Eq. (80) 

Subscripts 

1 = condition at the outer edge of the boundary layer 

0 = reference condition, Eq. (69) 

m = mean value 

Superscript 
0 = condition at x = 0 


INTRODUCTION 


ia HAS BEEN SHOWN by Stewartson! that the momen- 
tum equation for a compressible laminar boundary 
layer is transformed into that for an incompressible 
laminar boundary layer, assuming that the wall is 
thermally insulated, that the Prandtl Number is unity, 
and that the viscosity is proportional to the absolute 
temperature. If any one of the assumptions is dropped, 
the transformed equation is no longer of the same form 
as that for an incompressible flow, but it is still easy to 
solve because of the disappearance of variable density. 
Therefore it seems adequate to consider first the lam- 
inar boundary layer in incompressible fluids. 

The most frequently used method of solving the 
boundary layer equations approximately is the momen- 
tum-integral method, in which the momentum equation 
is integrated over the boundary layer thickness, and is 
hence satisfied only in average. By assuming a definite 
form for the velocity profile as a function of the normal 
distance, an ordinary differential equation is obtained, 
with the distance along the wall as the independent 
variable. In the original method due to von Karman 
and Pohlhausen,’ a fourth-degree velocity profile is 
assumed, whose coefficients are determined by the 
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boundary conditions at the wall as well as at the outer 
edge of the boundary layer. Consequently, the velocity 
profile is expressed as a member of a one-parameter 
family of curves, with 


A = (67/v)(du,/dx) 


as the parameter, which makes its appearance through 
the condition that the momentum equation should be 
satisfied at the wall. 

This method gives a reasonably accurate solution in a 
region of accelerated flow, but its adequacy in a region 
of retarded flow has been questioned. Separation of 
flow may actually occur where the method fails to give 
it. Various attempts have been made to increase the 
accuracy of the method. For example, Howarth’ used 
the velocity profile obtained from the exact solution for 
the case when the velocity outside the boundary layer 
decreases linearly with the distance; Schlichting and 
Ulrich‘ assumed a sixth-degree polynomial for the ve- 
locity profile, so that it could satisfy additional condi- 
tions at the wall as well as at the outer edge of the 
boundary layer. These refinements generally yield 
fairly good results when compared with the original 
method due to Pohlhausen. Strictly speaking, however, 
none of the methods is sufficiently accurate because the 
velocity profile is expressed as a member of the family of 
curves with the Pohlhausen parameter A, whereas the 
results of certain exact solutions suggest that \ cannot 
fix the velocity profile. In particular, \ should have a 
certain definite value at the separation point according 
to the assumption, but this is not the case according to 
the exact solutions, the value of \ at separation de- 
pending considerably on the velocity distribution out- 
side the boundary layer.*® 

It is reasonable to suppose that the difficulty may 
be reduced by using a velocity profile with two parame- 
ters, one of which is A and the other new. This has 
actually been done by Wieghardt,® who assumed an 
eleventh-degree velocity profile and determined the 
parameters by using the momentum-integral as well as 
the energy-integral relations, the latter of which is ob- 
tained by multiplying the momentum equation by 
velocity and integrating across the boundary layer 
thickness. Satisfactory results can be obtained, but 
the method appears quite tedious for practical applica- 
tion. An interesting modification has been proposed by 
Walz,’ who used a one-parameter family of curves again 
but gave up the fulfillment of the condition that the 
momentum equation should hold at the wall. The pro- 
cedure is thus considerably simplified, but results of 
comparatively high accuracy can still be obtained. 
Truckenbrodt* has recently extended the method so 
that the laminar and turbulent boundary layers may be 
treated along similar lines. 

It is to be pointed out that the improvement of ac- 
curacy is attained by dropping one of the boundary 
conditions at the wall and satisfying the energy-integral 


relation instead. Since the Pohlhausen parameter 
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enters only through the boundary condition, this js 
equivalent to saying that better results may be obtained 
by abandoning A for use as the profile parameter and by 
adopting another profile parameter that should be de. 
termined such that the momentum-integral and energy 
integral relations are both satisfied. The conclusion 
seems to be quite natural; because \ cannot fix the ve 
locity profile, as mentioned above. 

This method of approach is adopted in the present 
investigation. A simple fourth-degree polynomial js 
used as the velocity profile, quite the same as the 
original Pohlhausen scheme, but dropping one of the 
boundary conditions at the wall. Consequently, one of 
the coefficients of the polynomial remains undeter- 
mined, and it is used as the parameter for the velocity 
profile. Both the momentum-integral and energy-inte 
gral relations are considered, from which the relation 
between the profile parameter and the distance along 
the wall is determined. The analysis is given in the first 
part of the paper, and it is shown that satisfactory re- 
sults are obtained when compared with the exact solu 
tions. 

In the second part of the paper, the analysis is ex- 
tended to the laminar boundary layer in compressible 
fluids. The wall is assumed to be thermally insulated, 
but the other restrictions imposed by Stewartson are re- 
moved. The Prandtl Number may be different from 
unity, but the difference from unity is assumed small. 
The so-called Sutherland formula is adopted for the 
variation of viscosity with temperature, because it is 
considered desirable to use as accurate a formula as 
possible at high Mach Numbers where the variations in 
temperature are considerable. In order to simplify the 
analysis, however, a simple relation between temperature 
and velocity is assumed, the effect of Prandtl Number 
being taken into account approximately. The momen- 
tum equation is then transformed by the method pro- 
posed by Rott,’ which is a useful extension of the original 
Stewartson transformation. The resulting equation is 
no longer of the same form as that for the incompressible 
flow, but it is still amenable to solution by a similar 
technique as mentioned above. The accuracy of the 
solution is checked in the case of a flat plate by com- 
paring with Crocco’s results,’ and it is found to be 
satisfactory. The method of solution is then applied for 
estimating the separation point when the velocity out- 
side the boundary layer decreases linearly with the dis- 


tance. 


(1) LAMINAR BOUNDARY LAYER IN INCOMPRESSIBLE 
FLUIDS 


Basic Equations 


Consider first the steady two-dimensional laminar 


boundary layer flow of an incompressible fluid. The 
continuity equation is 
Ou/Ox + Ov/Oy = 0 (1) 
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APPROXIMATE SOLUTION OF 


and the momentum equation is 


Ou Ou du, O°u (9 
e —_ o } 


: oun | oy oo dx Oy" 
where x and y are the distances measured along and nor- 
mal to the wall respectively, « and v are the velocity 
components in the directions x and y, respectively, m; is 
the velocity outside the boundary layer, and » is the 
kinematic viscosity. Integrating Eq. (2) with respect 
toy from vy = 0 to y = 6, where 6 is the thickness of the 
boundary layer, and taking account of Eq. (1), we have 


the momentum-integral relation: 


6? du, 6* 26/Ou ‘o 
+ 92 2 4. = ( (3) 
Vv dx 6 it) Oy y=0 


Multiplying Eq. (2) by uw and integrating with respect 


lt} dé- 


vy dx 


toy irom y = 0 to y = 4, we have the energy-integral 


relation: 


u, de* * 6%? du, 40* f° (du\" 
6 = = dy (4 
vy dx vy dx uw," Jo \Ov, 
Here, 5*, 6, and 6* are the displacement thickness, 


thickness, and energy thickness, respec- 


Uu 
1 — ) ay (9) 

1} 
(1 = “) ay (6) 

uy 

u* " 
(: _ :) dy (7) 
u\~ P 


Following Pohlhausen,? we assume the velocity profile 


momentum 


tively, 
Ps 
° = 
0 
6 
u 
Oo Uy 
6 
a f u 
0 uy 
Velocity Profile 


in the form 


2 c yi 


Rf y _¢ ¢ 
=ata~ +d great (3) 


uu 

lt; 
and determine the coefficients ao, a, 6, c, and d, which 
are functions of x, from the conditions 


y = 0: u = 0 t (9) 
y= 5: u um, Uu/Oy = 0, Ou/dy?>=O0f ~ 


The usual condition 


0°u/Oy? = — (m/v) (du,/dx) (10) 


which states that Eq. (2) should be satisfied at the wall, 
is now dropped, so that one of the coefficients, say a, re- 
mains undetermined. This coefficient a is simply 
adopted as the parameter for the velocity profile, be- 
cause it possesses a clear physical meaning in that it is 
proportional to the shearing stress at the wall. We thus 


obtain 


ae (« sx 43 ) at (: a} (11 
= p= e 1 eS ) 
62 5 iil 5 
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Introducing the nondimensional quantities 
6*/5 = D, 6/6 = E, *~6=F 


H = §*/@ = D/E, G = 6* 6= F/E (13) 


26 / Ou 40* {-* (du\? 
( = P. dy =Q (14) 
u, \Oy/ y =0 u,- JO Oy 
we can write Eqs. (3) and (4) in the form 
uy de" 6° duy ” 
+ 2(2+ H) = Pp (19) 
y dx vy dx 
u, dG°*6" G6" du, 
+ 6 =@Q (16) 
vy dx vy dx 
where 
D = 2/5 — a/20 (17) 
E = 4/35 + @/105 — a*/252 (18) 
; S76 73 23 =; a 
F= — + —2e- -t# - (19) 
9005 9005 5460 2860 
P = Zak (20) 
O = [(4 35) F](48 — 4a + 3a’) (21) 


The quantities D, E, F, G, HW, P, and Q are functions of a 
only, and their numerical values are tabulated in Table 
- 


Quadrature Formula for Momentum Thickness 


It is readily shown that Eq. (15) may be integrated 
in the quadrature formula 


. 


#° : 
= Ag, tT! u," dx 
Vv 0 


provided that the distance x is measured from the point 
where 6 = 0, and that P — (3 + 2//)\ can be replaced 
by a linear function of \ such that 


(22) 


» 


P — (3+ 2H)\ = A — ma (23) 


where A and are constants and 


A = (6°/v)(du,/dx) (24) 


is the so-called Pohlhausen parameter. This method of 
approximate integration was put forward independently 
by Hudimoto" and the present writer'* early in 1941. 
The values A = 0.44 and n = 5 proposed by the present 
writer were obtained by evaluating the left side of Eq. 
(23) from Howarth’s solution for the case when the 
velocity outside the boundary layer decreases linearly 
with the distance.* Formulas of a similar type have 
also been found by Young and Winterbottom,'* Walz,"* 
and Thwaites,'® independently. 

It is not possible, however, to proceed similarly in 
the present calculation, because \ is no longer used as 
the profile parameter. Instead, Q is a slowly changing 
function of a, so that it may be replaced by a constant 


value Q,,. Eq. (16) is then integrated in the form 
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Fic. 1. Both sides of Eq. (35) are plotted as functions of a. 


Both curves have a local maximum at a = 4 but do not exactly 
coincide with each other. 


6? 
G? = Onu,—* f uy dx 
v 0 


Replacing, moreover, the function G by a constant value 


9 x 
6 8 6 _ 
= — u,° dx 
v Ga* 0 


As will be seen later, Q = 1.079 and G = 1.567 for the 
Using these values for Q,, and 


25) 


Gm, we have 


(26) 


flow along a flat plate. 


2 ee ” 
= 0.44 u,-* I u,> dx (27) 
v 0 


which is exactly the same as previously obtained from 
the momentum-integral. It was Truckenbrodt* who 
first found that the quadrature formula for the momen- 
tum thickness may also be deduced from the energy- 
integral, although he used the velocity profile obtained 
from the exact solution for the case when the velocity 
outside the boundary layer is proportional to the power 


Gm, we get 


of distance. 

We now consider Eq. (27) as the first approximate 
solution for 6*/v. We then calculate \ by Eq. (24). 
Therefore, we determine \ as the function of x. But, 
since it is a, and not A, that is used as the parameter for 
the velocity profile, it is necessary to obtain a relation 
that relates \ and a. 


Determination of the Parameter 


Eliminating d6?/dx from Eqs. (15) and (16), we have 


H(r- 2) 
2 G 


from which we can determine the relation between \ and 


din G 


(28) 
dln uy 


MH — 1) = 
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a. First, we neglect the second term on the right side, 
and we obtain the first approximation to the relatioy 
between \ and a. Then we estimate the neglected 
term, and proceed to the second approximation, and so 
on. 

After determining the relation between \ and a, we 
are in a position to check the accuracy of the quadrature 
formula (27). The function Q may now be plotted 
against \ and replaced by a linear function 
(29) 


Q0=Q, + rr 


where Q,, and r are constants, so that Eq. (16) is again 
integrated in the form 


9 B 
a, 
as ae ea uy" dx 
Vv ~ 0 


Thus we get the 


(30) 


where G is treated as dependent on X. 
second approximate solution for 6?/v. In most cases, 
however, the effects due to the variability of Q and G 
cancel each other, and the result obtained from Eq. (30) 
remains almost the same as that from Eq. (27). It 
turns out therefore that the simple quadrature formula 
(27) is a sufficient approximation. 

It is to be noticed that Truckenbrodt* derived from 
Eqs. (15) and (16) a differential equation for the quan- 
tity L = f§(H — 1)~'d In G, and solved it by an ap- 
proximate method. LZ is also a function of a, and the ad- 
vantage of using such a quantity may be readily seen 
from the Eq. (28) which, as a matter of fact, is es- 
sentially the same as the differential equation of Truck- 
enbrodt. However, there is a difference in the treat- 
ment of this equation: it is used by Truckenbrodt to 
determine LZ as a function of x, while in the present 
analysis it is used to obtain the relation between d and a. 
Another difference is that in Truckenbrodt’s solution 
the velocity profile is adopted which is obtained from 
the exact solution for the special case when the velocity 
outside the boundary layer is proportional to the power 
of the distance, while in the present analysis the velocity 
profile is assumed as a simple polynomial, which is more 
convenient to be extended to the case of compressible 
fluids. 


Particular Values of the Parameter 


As clearly indicated by Eq. (28), the relation between 
\ and a is not universal, but depends on the velocity 
distribution outside the boundary layer. However, 
there are a couple of particular values of a correspond- 
ing to certain definite conditions. For example, the 
separation point is characterized by the particular value 
a = 0. Other particular cases are the flow along a flat 
plate and the flow at the stagnation point. 

For the flow along a flat plate, \ = 0 and G = con- 


stant.t Therefore we have the condition 


t At the minimum pressure point of a flow with pressure 
gradient \ vanishes but G is not necessarily stationary. 
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APPROXIMATE SOLUTION OF 
QO 
P = — (31) 
G 
fom which we obtain 
a = 1.857 (32) 


as the particular value of the parameter for the flat 


plate. We have also 


Q= 


while the Blasius solution for the flat plate gives 


P = 0.439, 1.079, G= 1.567, H=2.60 (33) 


P=0.441, Q=1.090, G=1.572, H=2.59 (34) 


It appears therefore that the present approximate 
method predicts the flat plate characteristics with 
satisfactory accuracy. 

In order to find the condition at the stagnation point, 
we put # = O in Eqs. (15) and (16) and eliminate 
du,/dx from them. We then have 

P/(2+ H) = 0/3G? (35) 
as the equation for determining the value of a. Un- 
fortunately, however, it is found that there is no root of 
the equation in the physically significant range of a. 
Similar difficulties have already been experienced in 
connection with the refinement of the Pohlhausen solu- 
tion. It will be found, however, that the curves of 
P/(2 + HH) and Q/3G? plotted against a come very 
close to each other in the neighborhood of the point 
where both curves have a local maximum, as shown 1n 


Fig. 1. We tentatively take the value of a at this point 

a = 4.00 (36) 

as the value corresponding to the stagnation point. We 
then have 

P= i711, A = 0.084 (37) 


On the other hand, the exact solution for the flow u; = x 
gives!® 
(3S) 


P = 0.721, A = 0.085 


It is seen from this comparison that the approximate de- 
termination (36) is sufficiently accurate for practical 
purposes. 

Comparison with the Exact Solutions 


In order to check the accuracy of the solution as a 
whole, calculations are made for the case when the ve- 


locity outside the boundary layer is given by “, = ™,° X 
I x”), where 1," is the velocity at x = 0. Exact 
solutions have been obtained by Howarth* for the case 
n = 1, and by the present writer® for the cases” = 2, 4, 


and S. Results are shown in Fig. 2, where the Pohl- 
hausen parameter \ is plotted against the nondimen- 
sional shearing stress P, instead of against the profile 
parameter a. This method of representation is prefer- 


able, because there is no analogue to a in the exact solu- 
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It will 
be seen that the accuracy is satisfactory except for the 
case n = 8. Although it is difficult to extend the calcu 
lation right up to the separation point because of the 
(d In 1), it is 


tion, and Pa isa slowly changing function of a. 


considerable change in the value of (d In G) 
still possible to arrive at a point so near to the separa- 
tion point that the latter point may be extrapolated 
with sufficient accuracy. The result seems to be en- 
couraging, since the method can predict fairly ac- 
curately how the relation between \ and P (hence \ and 
a) changes in accordance with the velocity distribution 
outside the boundary layer. Such a change was demon- 
strated by the present writer,’ but until now there has 
been proposed no satisfactory solution which takes this 
effect into account. 

BOUNDARY LAYER IN COMPRESSIBLE 
FLUIDS 


(II) LAMINAR 


Basic Equations 


For the steady two-dimensional laminar boundary 
layer in compressible fluids, we have the continuity 


equation 
(0/Ox)(pu) + (0/dy)(pv) = O (39) 
the momentum equation 
Ou Ou du, ra) Ou . 
pu “a + pv = = ply rs + | =a (40) 


and the energy equation 
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oT 2 dity O oT 

plc, tT pvc, = piltlt; y k 7 
Ov Ov dx Oy Oy 


Ou \- 
M ( ) ( }] ) 

Oy 
where p is the density, yu is the viscosity, & is the ther 
mal conductivity, c, is the specific heat at constant 
pressure, 7’ is the absolute temperature, suffix 1 denotes 
the condition at the outer edge of the boundary layer, 
and other symbols have the same meaning as defined in 
Part I. wand & are considered as functions of 7°, while 
c, and the Prandtl Number 


o = uc,/k (42) 


are assumed constant. Since the pressure remains con- 


stant across the boundary layer, we have 
pm = 1,7 (43) 


Outside the boundary layer we have the isentropic re 


(") ‘i aaa 
po T» 


where y is the ratio of specific heats, and the suffix 0 de 


lationship 


notes a certain reference condition. 
It is assumed that there is no heat transfer through 
the wall. Therefore we have the condition 


y = 0: Ol ov = 0 (45) 


The Prandtl Number o is assumed to be slightly dif 
ferent from unity such that the square of the quantity 


| - Vo th) 


é 
may be neglected compared to unity. Finally, the 
viscosity is assumed to be given by the Sutherland 
formula 

mM poC (T/T) (47 
where 

C = (T/T) “[((To + S)/(T + S)] (4S) 
and Sis a constant which for air has the value 216°R. 

C is not assumed constant, as usually done, but 
treated as a function of the temperature, because it 
seems highly desirable to use as accurate a formula as 
possible for the variation of viscosity with temperature, 


Temperature Profile 


For a Prandtl Number of unity and zero heat transfer 
at the wall, the analysis can be considerably simplified 
because the energy equation (41) is replaced by a 
simple relation between temperature and velocity, viz., 


ri y—1u;? ( “) 
— =1+4 -(1-—— (49) 
T; 2 " uy," 
where ¢; is the local sound speed at the outer edge of 
the boundary layer. Assuming the Prandtl Number 
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slightly different from unity, we write 


- = 1+ ie g 
T; 2 C) 
where ® is an undetermined function of % 1, such that 


&® = 0 when uw 1, 0. This approximate expression is 


exact when o = 1, and, moreover, it gives the tempera- 


ture at the wall 


which seems to be supported by the results of calcula- 
tion due to Tifford and Chu" even in the presence of a 
pressure gradient. 

Assuming the temperature boundary layer has the 
same thickness as the velocity boundary layer,? and 
integrating Eq. (41) with respect to y from y 0 to 
y = 6, we have 


- . 1 fe 
/ puc,(T — 7)) dy = ; pile ( iy u-)dy (52 
J/0 Jo 


Substituting Eq. (50) we obtain 


. ; ‘ “7 \ 
/ pum dy / pit ( l } ay O83 
0 0 \ lt) 7/ 


which yields an integral condition to be satisfied by the 


function ®. 


Stewartson-Rott Trensformation 


It was shown by Stewartson! that the momentum 
equation (40) may be transformed into that for an in 


compressible fh 


ud, provided that there is no heat trans 
fer, the Prandtl Number is unity, and the viseositv is 
proportional to the absolute temperature. The trans 
formation was modified by Rott® such that the effect of 
Prandtl Number is taken into account, by assuming the 
temperature profile in a form similar to Eq. (50) where 
the term &® is omitted. The transformation is given by 


l 


T\ayve 
y=(2) / P ay | 
To /S0 pr ij 


where x and y are the coordinates of the original flow, 
and X and JY are the coordinates of the transformed 
flow. 

Applying the transformation (54) to the momentum 
equation (40), and making use of the temperature pro- 
file (50) we obtain 
U(OU/OX) + VOU /OY) = 

(1 + EB) [(U,(dU1/dX)] + (0/0V)[C(OU/OY)] 


(55) 


t See also the subsequent footnote. 
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APPROXIMATE SOLUTION OF 


where (’ and I are the velocity components in the direc 
tions \ and J, 1s the velocity outside the boundary 
laver in the transformed flow, and 
a l ) 
a } - ~ 
B Mm 
V al(4 1) /2] (a7 ‘¢ 
There exists a relation between (; and mu, such that 
l 7 T) sie al D7 
and, moreover 
U/U, = u/u (SS 


Eq. (55) is no longer of the saine forin as Eq. (2) for 
the incompressible flow, because the restrictions im- 
posed on the Prandtl Number and the viscosity-tem 
perature relation have been removed. It is still easy to 
solve, however, since the variable density has dis- 
appeared. 

Integrating Eq. (55) with respect to VY from V = 0 
to } = Awhere Ais the thickness of the boundary layer 
in the transformed flow, we have 


02 dl, a>. ww ; 
9 ( 9 + |. | b dj = 
0 0 0 


dX 
~ (3) 7 
(09) 
U;, \OY/y ~¢ 


Multiplying Eq. (55) by U’ and integrating with respect 


U; d8? 
— _ 
yw aX Vy 


to } from Y = Oto Y = A, we have 
2 ¢B a U , 
T 3 6 j "7" ) = 
40* [° /au\? ._ 
= C ( -_} dY (60) 
l 1” 0 Oo} 


These equations are the analogues of the momentum- 
integral relation (3) and energy-integral relation (4) 
respectively, and A*, 0, and 0* are the displacement 
thickness, momentum thickness, and energy thickness 
in the transformed flow, respectively, 


U, do* 2 4 6 e*? dU, (: 
» ax te an 


BOT 


NDARY-LAYER EQUA 


Velocity Profile 
velocity profile in the form 
U y? ( } y y ) 

= 6—S + 3 + a (1 64 
U; A- A A A A 


where a is the profile parameter. 


We now assume the 


This is exactly the 
same as used in incompressible flow, Eq. (11 

l, by 
virtue of Eq. (5S), they may be considered as functions 
of Y/A. 


Since C and ® are functions of u ™, or of l 


We assume therefore 


C = ao + ay( V/A) + aro( Y?/A?) + a3( V%/A (O65 
®@ = B + B( V/A) + BY? A?) + 
B3( ¥3/A*) + 6,(¥4/ A‘) (66 


where the coefficients, a's and #’s, are determined by 


the conditionst 


Y=0: C=C, oOC/OY = @, 
@=0, 0¢/d0Y =0 
(67 
Yu A: C=C, 0C/eY =% 
@=0, 0¢/0Y =0 


and 

+ These boundary conditions are based on the assumption that 
the temperature boundary layer has the same thickness as the 
velocity boundary layer. This is not correct except when the 
Prandt!] Number is unity 
tion is considered small, because the function C is multiplied by 
(0U/O¥)? which is small near the outer edge of the boundary 
layer, and also the function ® is multiplied by £ which is assumed 


However, the error due to this assump- 


small in the present analysis. 





TABLE | 

a D E F G In G H P @) Q/G? J K 

0 0.4000 0.1143 0.1759 1.531 0.4262 3.500 0 0.960 0.4094 0.333 2.382 
0.25 0.3875 0.1164 0.1784 1.532 0.4268 3.328 0.0582 0.962 0.4098 0.317 2.330 
0.50 0.3750 0.1181 0.1812 1.835 0.4285 3.176 0.1181 0.968 0.4109 0.298 2 282 
0.75 0.3625 0.1192 0.1834 1.539 0.4311 3.041 0.1788 0.979 0.41338 0.278 2.238 
1.00 0.3500 0.1198 0.1859 1.544 0.4345 2.921 0.2397 0.994 0.4169 0.257 2.198 
1.25 0.3375 0.1200 0.1860 1.550 0.4383 2.813 0.3000 1.014 0.4219 0.236 2 161 
1.59 0.3250 0.1196 0.1862 1.557 0.4425 2 716 0.3589 1.038 0.4282 0.214 2.125 
1.75 0.3125 0.1188 0.1858 1.564 0.4471 2 630 0.4158 1.066 0.4357 0.198 2 092 
2 00 0.3000 0.1175 0.1845 1.571 0.4518 9 554 0.4698 1.097 0.4448 0.173 2 060 
9 95 0.2875 0.1156 0.1825 1.579 0.4566 2.486 0.5203 1.130 0.4536 0.154 2 029 
2.50 0.2750 0.1133 0.1797 1.586 0.4613 2.427 0.56965 1.165 0.4632 0.137 1.999 
3.00 0.2500 0.1071 0.1714 1.600 0.4700 2 333 0.6429 1.234 0.4821 0.107 1.940 
3.59 0. 2259 0.0990 0.1595 1.611 0.4767 3. Sis 0.6931 1.290 0.4970 0. O84 1 87 
100 0.2000 0.0889 0.1436 1.615 0.4796 2.250 0.7111 1.313 0.5931 0.067 1.812 
$59 0.1750 0.0768 0.1235 1.608 0.4752 2.279 0.6911 1.281 0.4952 0.054 1.731 
5.00 0.1500 0.0627 0.0989 1.578 0.4551 9 392 0.6270 1.165 0.4678 0.046 1.622 
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1-4 T l : 1 / 584 146 19763 ” 
| wad ne ~ 57915 aa”) — 
Crocco LE \S140 VtILO Olid 
3 ~~ Chapman-Rubesin | J = (i Q)} 16/35 — (13/105)a + a? 84} 74 
“—_ 7 Present method K and J are also dependent only on a, and their numeri- 
CyVRe ha 2 cal values are given in Table l. 
1-2 nike Eqs. (71) and (72) are almost identical with Eqs. 
a a (15) and (16) for incompressible flow, the only dif- 
T,= 648 R ference being in the term proportional to €B and the term 
i proportional to (C; — 1). These additional terms repre 
sent the effect of Prandtl Number different from unity. 
and the effect of viscosity-temperature relation dif- 
| ferent from a simple proportionality, respectively. 
1-0 
O | 2 M Approximate Method of Solution 
Fic. 3. Skin-friction coefficient of a flat plate plotted against Eqs. (71) and (72) may be solved in exactly the same 


the Mach Number. Free-stream temperature 648°R., Prandtl 


Number 0.72. 



































1-5 | | 
1-4 — 
Cs Re 
¢ Crocco 
i: en eee Chapman-Rubesin 
° Present method 
1“ | | | 
0 | 2 3 4 5 
M 
Fic. 4. Skin-friction coefficient of a flat plate plotted against 
the Mach Number. Free-stream temperature 72°R., Prandtl 


Number 0.72. 


oo 
| @ — dY = 0* (68) 

0 U; 
the last condition being obtained from Eqs. (53) and 
(63). 

The reference condition 0 is now chosen such that the 
temperature 7) is given by the relation 

I) = T;}1 + V ally — 1) 2] (uy? c7) (69) 


In other words, 7) is equal to the wall temperature 
[Eq. (51)]. With this definition we have 


Co = 1, Cr = (11/70) “(To + S)/(T1 + S)] (70) 
Eqs. (59) and (60) are then written in the form 
U; de? 2 6? dU, 
2(2 + H + tBK) =P (71) 
Vo dX r sd r Vo dX 
U, dG°e? 2 ) G0? dU, 
61 EB = 
vw dX i ( t yw aX 
Q{1 + J(Gi — 1)] (72) 


respectively, where G, H, P, and Q are given by Eqs. 
(13), (17), (18), (19), (20), and (21), and 


way as in incompressible flow. The quadrature formula 
for the momentum thickness in the transformed flow is 


rr. &.. fs”... ; ' 
Se l — uy [1 + J,(C 1— 1) } dX (75 
: 0 


Vo G n- 


which is obtained from Eq. (72) by neglecting the term 
(2/3)éB and replacing Q, G, and J by constant values 


corresponding to the flow along a flat plate.j Then, 


eliminating d0?/dX from Eqs. (71) and (72), we obtain 


A[H — 1+ éB(K — 2)] = (1/2)P — (Q/2G") 
[1 + J(C, — 1)] + Al[(d In G)/(d In U,)] (76) 

which yields the relation between the profile parameter a 

and the Pohlhausen parameter in the transformed flow 


A = (0?7/)(dU;/dX) (77) 


Eqs. (75), (76), and (77) correspond to Eqs. (26), (28), 


and (24), respectively. 


Skin Friction on a Flat Plate 


In order to check the accuracy of the method, the skin 
friction on a flat plate is calculated and compared with 
that given by other solutions. First, the parameter a is 
determined from the condition 


P = (Q/G*)[1 + J(G, — 1)] 78 


which is obtained by putting \ = 0 and G = constant in 
Eq. (76). Next, the momentum thickness 0 in the 
transformed flow is determined by integrating Eq. (71 


in the form 
6? = P(X /U;3) 


(79) 


Then the coefficient of skin friction based on the free- 
stream characteristics is given by 


C, = 4 ( ay 7 
1 pum? \" dy a 


} For most practical purposes, values corresponding to incom- 
pressible flow (a = 1.857, Om/Gm? = 0.44, Jm = 0.184) may be 
used with sufficient accuracy. 


ov P 
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APPROXIMATE SOLUTION OF 


where Re = “x, » is the Reynolds Number based on the 
distance from the leading edge. 

Values of C,V Re are calculated for various values of 
Mach Number MV = u/c; assuming y = 1.4 and o = 
0.72 (€ = 0.15). Results for the free-stream tempera- 
ture 7, = 648°R. and 7; = 72°R. are represented in 
Figs. 3 and 4, respectively. Included for comparison 
are results calculated by Crocco" and by Chapman and 
Rubesin.'* Crocco’s results are considered to be most 
reliable because he used the Sutherland formula (47), 
while Chapman and Rubesin replaced C by a constant 
value. It is seen from this comparison that the present 
method gives results which are sufficiently accurate for 
practical purposes. It is also to be noted that the 
effect of Prandtl Number is rather small compared with 
the effect of viscosity-temperature relation different 
from a simple proportionality. 


Compressibility Effect on Separation 


Transforming back to the original flow, we obtain 
from Eqs. (75) and (77) 
2, - 1 
yO =» 2 bd 
Qm (To\ 1-1 ~°V% du r 
abet &— u,~* + 
Gat Ts dx 0 


re 7~5  7Ve . 
IuiGs —_ 1)} (=) u,° dx (81) 


where 7) 7, is given by Eq. (69). This formula enables 
us to calculate A directly from the velocity distribution 
in the original flow. The relation between A and a may 
be determined from Eq. (76), where In Ll should 
preferably be replaced by In m4 — (Yo/2) In (7), T)). 
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Fic. 5. Relation between the Pohlhausen parameter A in the 


transformed flow and the nondimensional shearing stress P 
for the velocity distribution u,; = ™°(1 — x) 
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Fic. 6. Effect of compressibility on separation for the velocity 
distribution u, = ™°(1 — x). The velocity ratio at the separa- 


tion point is plotted against the Mach Number. 


As an example of the application of the method, the 
calculation was made for the boundary layer when the 
velocity outside the boundary layer is given by 


UW, = u,°(1 — x) (82) 


where 7," is the value of u, atx = 0. It is assumed that 
y = 1.4, o = 0.72 (€ = 0.15), and the free-stream 
temperature at x = Ois 7)" = 648°R. and 7)" = 72°R., 
respectively. Results are shown in Fig. 5, where the 
Pohlhausen parameter A is plotted against the nondi- 
mensional shearing stress P for various values of the 
Mach Number M° = u,°/c,°._ As in the skin friction on 
a flat plate, the effect of compressibility is greater for 
the case 7° = 648°R. than for the case 7,° = 72°R. 
The separation point is determined by locating the 
value of x for which A attains the value corresponding to 
P = 0. Results of this calculation are shown in Fig. 6, 
where the velocity ratio “;/u,° at the separation point 
is plotted against the Mach Number 17°. Included also 
is the result calculated by Stewartson' by assuming that 
the Prandtl] Number is unity, the viscosity is propor- 
tional to temperature, and the value of A at separation is 
the same as for incompressible flow. In spite of these 
assumptions, Stewartson’s results appear to be com- 
paratively good; this is because of the fortunate cir- 
cumstances that the separation is delayed by the effect 
of viscosity-temperature relation different from a simple 
proportionality, while it is hastened by the effect of dis- 
tortion of velocity distribution in the transformation 


(54) and (57). 
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A Note About the Effects of Product of Inertia 
on Lateral Stability 


C.N. Tsu 
Design Speciaiist, Northrop Aircraft, Inc., Hawthorne, Calif 
February 18, 1954 
- REFERENCES 1-3, the importance of the products of inertia 
on lateral stability of aircraft has been shown. This note 
points out a step in evaluating the effects of the product of inertia 
on aircraft dynamic lateral stability which may simplify to some 
extent this evaluation. 
The lateral equations of motion from reference 1 are repeated 
below for convenience. 


(J,D? —- Ci;D )o — (/,,D* + Ci;D)y¥ _ Cig 8 =QVQ (1 

(1D? — CuyDW — UIrD? + CigD)6 — Crg8 = 0 (2) 

[((2bu/V) (Dy + DB)|] — (CL + Cy 3D )o — (Cy, tan y + 
CyyD)¥ — Cyg8 = 0 (3) 


It is observed from the above equations that the product of inertia 


terms J,,D*y and J,,D*o are only in Eqs. (1) and (2). By the 
manipulations of [/.(1) + J:(2)] and [J,(2) + J,(1)] Eqs 
(1) and (2) now may be rewritten as (1’) and (2’). 
(1,D? -— Ci;,’D — €,,;'Dy - Ci,'B = (0 er’ 
(1,D? — Cny D)y —- Cng Do - Cng’B = 0 (2 
where 
Cig’ = (Cry + (Ue2/Te)Cng\/{l — (ee/Tz) e2/Iz)]} (4) 


Ta)Cngl/{1 — (Tse/Te) (Iee/Te)}} (5) 


Cig’ = [Cig + (Tee/Te)Cugl/{1 — [(Te2/Iz) Ie2/I2)]} (6) 
Cry’ ~ [Cry + (Le Tz) Cig] bli [(Es/ls) C/E C7 
Cug’ = (Cng + Ue2/Tz) Cig] it = (../f.) Us/li)It (8 
Cug’ = (Cug + (Uze/Tz)Cig\/ — [(Ie2/Te) (ee/I2))}{ (9) 


Thus the form of Eqs. (1’) and (2’) become identical to that of 
Eqs. (1) and (2) with /,., the product of inertia, made zero. Con- 
sequently, the effects of product of inertia on dynamic lateral 
stability may now be evaluated both qualitatively and quanti- 
tatively through the use of equivalent lateral stability deriva- 
tives as defined in relationships (4) through (9). Known effects 
of stability derivatives on dynamic stability may now be utilized 
to anticipate the effects of product of inertia. Stability boundary 
plots such as contained in the references may now be plotted in 
terms of these equivalent rolling moment and yawing moment 
derivatives with side force derivatives remaining unaltered. 
When rudder and aileron controls are involved, it can be shown 
that the resulting equivalent yawing and rolling moment control 


$96 


derivatives are related to the natural derivatives in the same way 


as the stability derivatives 
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Concluding Remarks on Detached Shock Waves 


Robert Wesley Truitt 

Professor of Aeronautical Engineering, Virginia Polytechnic 
Institute, Blacksburg, Va. 

February 18, 1954 

that 


Laitone’s statement 


ig REFERENCE 1, | took issue with 
4 was not 


Mocckel’s result? for shock 
valid for thin wedges near Mach Number one. It 


detachment distance 


was shown 
quite clearly in Fig. 1,' that Moeckel’s values for semiangle @ = 
$.5° at Mach Numbers near unity are in excellent agreement 
with experiment, whereas Laitone’s formula showed extremely 
poor agreement. 
Laitone is correct when later,’ in explaining his poor agreement 
with experiment, he concludes that his own formula [Eq. (3)*]: 
. rightly may be restricted to Mach Numbers less than 1.01! 
For Mach the results 


Laitone’s Eq. (3)*] do not apply 


Numbers greater than this—i.e., when 


there cannot be any important 


effect of nose shape, ... 
It should be pointed out, however, that earlier in Fig. 7,° Lai- 
with Moeckel’s result* 


< 15 


tone does apply and compare his Eq. (3), 
in the free-stream Mach Number range (1.1 < lJ) and 
Fig. 7* proves there is a considerable effect for 


10°, 15°) near Mach Number of unity.” 


concludes: 
slender wedges (6 = 5°, 
It is clear that, at that time, Laitone considered Mach Numbers, 
(1.1 < Wp < 1.5), sufficiently near Mach Number one to com- 
pare his results with Moeckel’s. Laitone also shows /is experi- 
mental results in Fig. 7° to be in accord with his Eq. (3)! for the 
same Mach Number range. 

In view of reliable experimental results shown in Fig. (1),' 
it must be concluded that Laitone’s theoretical and experimental 
results shown in Fig. 7,5 are not only incorrect but misleading 


REFERENCES 
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Plane Bodies, Readers’ Forum, Journal of the Aeronautical Sciences, Vol. 20, 


No. 9, p. 651, September, 1953 
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Factors Affecting Transition Reynolds Numbers 
on Models in Supersonic Wind Tunnels 


ohn Laufer 

enior Reseerch Engineer, Jet Propulsion Laboratory 
tute of Technology, Pasadena, Calif 

February 19, 1954 


California 


the transition 


the 


that 


lavers on 


PUBLICATIONS have indicated 


R' CENT 
Revnolds > 


me in different wind tunnels varied greatly in no obviously pre 


umber of boundary same 10 


ictable fashion. The great scatter of the experimental data was 
believed to be caused by disturbances of various degrees present 
in the tunnels 

These investigations indicate clearly that if a quantitative 
study of the mechanism of the transition phenomenon is to be 
undertaken in a wind tunnel, the elimination or understanding of 
perturbation fields of the flow in wind tunnels is essential 

In general, these disturbances found in supersonic wind tunnel 
nozzle flows may be of three different types*: (a) vorticity fluc- 
tuations (or transverse waves); (b) sound, or Jongitudinal waves 
including weak compression waves emanating from the imperfect 
nozzle walls; and (c) temperature fluctuations. Although these 
disturbances are coupled especially at higher Mach Numbers, one 
can make some statement as to their possible origin in wind tun- 
nels. For instance, a considerable amount of vorticity fluctu- 
ations is generated by blowers, turning vanes, honey combs, etc 
Longitudinal wave type disturbances can originate from poorly 
fabricated nozzle walls or possibly from turbulent shear layers.* 
Finally, temperature fluctuations might be caused by inadequate 
mixing of air in coolers or by large local turbulent dissipation 


Unfortunately, at the present time no adequate experimental 




















method is available that would permit a direct study of these 
effects 
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FORUM 


The method of approach to this problem adopted at the Jet 


Propulsion Laboratory of the California Institute of Technology 


is an indirect one. It is proposed to attempt a successive elimi 


d three contributions to the pertur 


nation of the above mentione 


bation field in the supersonic tunnel Strictly speaking, of 


course, this is not possible because of the interaction existing 


However, since one has some physical picture of 
might 


between them 


the origin of the large st part of these disturbances, one 


artificially change the conditions under which they are generated 


Specifically, the variation, if any, of the boundary laver transi 


Reynolds Number 


model under two conditions; 


tion could be observed on a wind tunnel 


first, when disturbances of tvpe (a 


settling 


are reduced by installing damping screens in the tunnel 
) 


chamber; second, when they are greatly increased by the in 


stallation of a grid before the contraction. Similar experiment 
can be carried out involving disturbances of types (b) and (« 
At present only the first set of experiments, as described above 
this L 
Phe 


two conditions were 


has been carried out at iboratory in the IS in. by 20 in 


supersonic wind tunnel turbulence levels in the settling 
and 7 per 


Mach 


tunnel 


chamber under the 0.6 per cent 


These 


Number rang¢ 


cent values were approximately constant for the 


1.4 to 4.0, and nearly independent of the 


stagnation pressure rhe temperature fluctuations were negli 
gibly small 

A 5° included angle cone, 100 em. in length, was used as the 
wind tunnel model The method of detecting transition was 
similar to that reported in reference 1 rhermocouples were 


placed 2 cm. apart and surface temperatures were recorded. The 
sudden rise in temperature indicated transition 
Some typical results are given in Figs. 1, 2, and 38 The values 


of the ordinate are proportional to the difference of the model sur 
face temperature, 7, and the temperature at the outer edge of the 
laver, 7 The Number 
of free stream conditions and wetted length on the conc 
the factor A, 


for the laminar portion of the 


boundary Revnolds is defined in terms 


surface 
The results were normalized by to vield identical 
values of temperature difference 


boundary layer 


The results may be summarized as follows: In the low Mach 
Number flows (.1/ 2.5) the turbulence level of the settling 
chamber had a strong effect on the boundary laver transition 


Revnolds Number However, no such effect could be detected 
for flows at MW > 2.5 was further substantiated by 
similar experiments conducted in the JPL 12- by 


working section of this tunnel has 


This result 
12-1n. super 


sonic tunnel. Thc a much 


condition and a contraction ratio only one 


the 18- by 2‘)-in 


different entrance 
third 
on the transition Reynolds Number of a flat plate 


as large as tunnel. Some experiments 


boundary layer 
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conducted in the 18- by 20-in. tunnel were also consistent 
with the above conclusion. 

It is interesting to note that in the Mach Number range 1.4 - 
M <4 the transition Reynolds Number decreased monotonically 
with increasing Mach Number for constant Reynolds Number 


per inch, a trend also reported in reference 5. 
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Transonic Flow Past Wedge Profiles by 
Hydraulic Analogy 


E. V. Laitone* and Helmer Nielsent 

University of California, Berkeley, and Ames Aeronautical 
Laboratory, NACA, Moffett Field, Calif., Respectively 

February 19, 1954 


B’ UTILIZING the analogy existing between the gravity waves 
propagated in shallow water and the two-dimensional flow 
of a perfect gas with a specific heat ratio of 2, an experimental 
study was made of the transonic flow about symmetrical slender 
wedge profiles. 

In the hydraulic analogy the local water depth corresponds to 
both the density and the temperature of the transonic gas flow. 
This permits direct evaluation of the equivalent gas flow velocity 
distribution by simply measuring the variation of the water depth 
However, since the tests were made by towing the model through 
'/;-in. deep water (the best conditions for the analogy'), it was 
practically impossible to measure the surrounding water depths 
by pressure measurements or simple depth gages. Consequently, 
the water heights were evaluated by measuring the variation in 
the electrical resistance of 0.004-in. diameter piano wire probes 
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These water height measurements were used to compute the 
velocity distribution along the stagnation streamline from the 


detached shock wave to the nose of the wedge. The experimenta| 
data shown in Fig. 1 were thus obtained by Nielsen,? and ar 


presented in terms of the standard pressure coefficient defined a 


Cp = (pb — pi)/[C1/2) a1 Vi?) = (2/yM?) [(p/p1) - 1 (] 
or corresponding to the hydraulic analogy with a specific heat 
ratio of y = 2as 


Cp = (1/M,?) [(A/m)? — 1] = (ghi/ Vi?) [(A/m)? — 1] (2 


The data presented in Fig. 1 were those obtained with a wedge 
base width of */, in. (Y, = One-half base width = */, in.) and 
with semivertex angles (6) of 5° and 20 The distance Y was 
measured along the stagnation streamline ahead of the wedge as 
shown in Fig. 1. 

Fig. 2 summarizes the experimental data of reference 2 in terms 
of the detachment distance (Y = ZL) measured to the foremost 
position (normal shock portion) of the detached bow shock wave 

The principal aim of these tests was to check the validity of the 
simplifying assumptions of the small perturbation theory pre- 
sented in reference 3 for determining the location of the de- 
tached bow shock wave as the Mach Number approached unity 
The theory developed in reference 3 was based on the observation 
that in the limiting case as the Mach Number approached unity, 
the detached bow shock wave and the sonic line became nearly 
parallel lines, far apart, and both extending an unlimited distance 
from the body. Then the interaction between the sonic line and 
the detached shock would be negligible, so the velocity variation 
along the stagnation streamline should be practically the same 


as that produced in a subsonic flow having \/7, M, = (2 - 
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VW Then it was shown that the shock wave detachment dis 


tance would be determined as the distance ahead on the stagna- 
tion streamline where the velocity drop through the normal shock 
portion (with 4, ahead and \/, behind) produced the required 
velocity ratio corresponding to the subsonic flow velocity vari 
ation along the stagnation streamline. This vields the expression 
for the detachment distance L from the shoulder of the wedge as 


> te 


(4V 2 w/v + 10) (M4 - 


In [1 — ( @L)|~-! = 


which agrees with the transonic similarity relation. The experi 


rest wedge 


mental data on the detachment distance fcr the sk 
tested in reference 2 compared favorably with this theory, as 
shown in Fig. 2, but, as would be expected, only for Mach Num 
bers approaching unity. 

Fig. 1 compares the velocity or pressure coefficient variation 
used in the small perturbation theory of reference 3 with the ex 
the measured local water 


perimental values calculated from 


depths. It is seen that there is some justification in assuming that 
the flow variation behind the detached shock is similar to that 
in subsonic flow, but in agreement with the basic assumption, 
only for the slenderest wedge at Mach Numbers near unity. At 
the higher Mach Numbers or larger nose angles the velocity gra 
dient became radically different from the subsonic variation as- 
sumed in reference 3, as would be expected. 

It is interesting to note that the experimental work being con- 
ducted at the present time indicates that as the channel walls are 
moved nearer to the model, the detached shock moves farther 
ahead of the nose, especially for Mach Numbers near unity. The 
measured distance, just before choking or a strong Mach-reflec- 
tion occurs, tends to approach Moeckel’s* result, which predicts 
a detachment distance independent of the wedge angle. 
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Errata to ‘‘The Calculation of Wall Shearing 
Stress from Heat-Transfer Measurements in 
Compressible Flows’’! 


Nick S. Diaconis 


Lewis Flight Propulsion Laboratory 
May 6, 1954 


EK* (9), (12), AND (13) are incorrect 
of Eq. (6) results in the following 
1)/2)M 


oe f, ' + ae! | f e~** dis 


temperature distribution has no effect 


NACA, Cleveland hic 


The solution 


correct 


2 Kn? 


This correction in the 


on the expression for the Nusselt Number, Eq. (16), or the con- 
clusions 
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Aerodynamic Coefficients of an Oscillating 
Airfoil in Two-Dimensional Subsonic Flow 


van dz Vooren, and J. H. Greidanus 
titute, Amsterdam- 


R. Timman, A. | 


National Aeronaut cal Research Ir W Holland 


February 7, 1954 


|. ernie VALUES for the aerodynamic coeflicients of an 
oscillating wing in two-dimensional subsonic flow, calcu 
lated by aid of an analytical solution in terms of Mathieu func 
tions, have been published by the authors in THis JOURNAI 

Considerable question regarding the accuracy of these values was 
put forward by Fettis* and Jordan‘ with the argument that the 
reciprocity relation, which was not known at the time that the 
authors published their data, was not satisfied. Consequently it 
was investigated whether the formulas giving the final solution 
were in agreement with the reciprocity relation. This 
to be the case as has been shown by Timman,’ who also simplified 


appeared 
the formula used for calculating the strength of the pressure 
singularity at the leading edge. The reason that the numerical 
values did not satisfy the reciprocity relation was that in spite 
of numerous checks one small error had been made in the numeri 
cal evaluation. Meanwhile, this error has been eliminated by 
the Mathematical Centre and revised values satisfying the rec- 
iprocity relation up to a few units in the fifth decimal figure, are 
now available. The values for the wing alone are given sub- 
sequently, but values for the wing with control surface and for 
the same parameter values as in reference 2 may be applied for 
at the above mentioned address 
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Al = amplitude of translation in mid-chord point, posi- 
K = aerodynamic force, positive downward tive downward 
M = aerodynamic moment about mid-chord point, posi- B = amplitude of rotation, positive if trailing edge down- nd 
tive tail-heavy ward 
py = air density (in undisturbed state ) ka’, Ra” = real and imaginary part of ka, etc 
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Characteristic Solution for Axially Symmetric 
Transonic Flow 


rik Mollo-Christensen 
Research Engineer, Department of Aeronautical Engineering, 
Massachusetts Institute of Technology, Cambridge, Mass 


February 20, 1954 


INTRODUCTION 


J D. CoLe has given the expression for simple waves in slightly 
e perturbed two-dimensional transonic flow.' It can be shown 


that the potential corresponding to Cole’s solution is of the form 


ox, y) = X(x) V(y) 


where V(x) and }(y) are functions of x only and y only, respec 
tively 

Supposing that the potential equation for axially symmetric 
flow has a similar separated solution, and proceeding to seek it, 
me finds that such a solution exists, and is a characteristic solu- 
ton 

Both 
ound a bend, and the flows are explicitly given in the physical 


the characteristic solutions correspond to expansion 


plane 


(1) Two-DIMENSIONAL TRANSONIC CHARACTERISTICS 


The equation for the perturbation potential ¢ is [2 


Set 


g(x, y) = X(x) V(y) (2 
where V(x) and }(y) are functions of, respectively, x and y alone. 


Substituting in Eq. (1), one obtains: 


yt 1 X(x)X"(x) V"(y) 
; = . (3) 
a* X (x) [¥(y)]? 
Solving the two equations 
a* (x + ¢)8 
(x,y) = - (4) 
Ky + 1) (y + ce)? 


where c; and c. are arbitrary constants. To this solution can be 


idded the obvious solution of Eq. (1): 


!, being a constant. 
The components of perturbation velocity are: 


a* x+a\? 
u=>¢, = 
+i\y +a 


ind 
Sa" x + “) 
v=y=- \4 
3 7 + l ¥ + Ce 


ind 7 can be expressed as a function of u alone, viz., 


+ 1 
u (8S) 


which is the expression given by Cole for the transonic character 
istics. We note that the potential given by Eq. (4) vanishes at 
t= —c, as do its first derivatives, which correspond to uniform 
sonic flow. It will also be noted from Eq. (8), that when u is 
negative, v is imaginary, so the solution is only valid for super- 
sonic local Mach Numbers. 
(2) AxIALLY SYMMETRIC FLow 

Let ry be the radial distance. The potential equation for axially 

symmetric flow over a slender body of rotation is 
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FORUM 
y+ 1 l 
rrr = ~ Or + ¢ (9) 
a* r 
Trying a solution of the form 
g(x, r) = X(x)R(r) (10 
yields, upon substitution, 
y+ 1 X(x)X"(x) _ r 
. = - R(r) + R"(r)/[R(r)]? (11) 
a* X(x) r 
Solving the two ordinary differential equations yields 
2 a® (xX + G1)* 
Ax, rT) = (12) 
Q « +. | r 
The velocity components are 
Ss 2 x+oca\? 
“=, = 
37+ 1 r 
(13) 
: 2» x+Ca ' 
ial hal Qg Y T l r 
And they are interrelated as follows: 
\/2/3 ,¥+1, , 
v= t 2/3 u (14 
\ a” 


The solution found is therefore a characteristic solution 

The solution is not valid when r vanishes since then there can 
be no further opportunity for the flow to expand. The solution 
is also restricted to supersonic flow, because of the square root in 
Eq. (14), which makes v imaginary for negative u 
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On the Interaction of a Plane Shock Wave with 
a Wedge 
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(1) INTRODUCTION 


ET A PLANE SHOCK WAVE of any strength traveling horizontally 
| F into a still air, encounter a corner formed by two plane solid 
It is known by experiments that the reflection process 
regular reflec- 


walls. 
takes either of two qualitatively different forms: 
tion or Mach reflection. Unsteady flow takes place within the 
region bounded by reflected shock (and Mach shock) and the 
walls. 

When the viscosity and heat conductivity of the air are neg- 
lected there is no intrinsic length in the problem; the flow pat- 
tern remains similar to itself (i.e., is pseudo-stationary) and may 


Here, R 


and @ are the plane polar coordinates with the corner as origin 


be described in terms of the variables R/t and a@ alone 


and ¢ is the time measured from an instant the incident shock 
struck the corner. 

Extensive studies of the problem have been carried out experi- 
mentally.! However, theoretical investigations of the field of 
flow have so far been confined to cases when the vertex angle 26 
of the wedge is small or when it does not differ much from z. 
This paper provides a preliminary account of an attempt to solve 
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approximately the pseudo-stationary flow field for the cases of 
arbitrary values of 6 and shock strength 

Solutions of the linearized equations for velocity components 
and pressure are obtained which enable us to solve the boundary 
value problem. It is shown that a singularity of rapid expansion 
takes place on the line bounding the region of influence of the 
corner on the flow field. This effect is clearly verified+ by the 
experimental results of White.” 


(2) EQUATIONS FOR PSEUDO-STATIONARY FLOW 


Let the velocity components in RK and a@ directions, pressure, 
density, and velocity of sound at any point be U, V, P, Q, and C, 
respectively, and let the quantities at a certain standard point 
be denoted with suffix 2 which is moving uniformly in such a way 
that the state is constant in time. If now the transformations 
R/Cot = 1, U/C, = M+ 4, V/C2 

U2/C, = M, V2/C, = N, c 


: P/(%2C?) = p, Q/Q: = p 


be made, and if the fact that u, v, c, p, p depend on r and a@ be 
used, the equations of conservation of mass, momentum, and en 


tropy become as follows: 


Op N+v70p 
(M+u-—r) -+ 
or r Oa 
E M + u 1 ov 
p + + = (0 (1) 
Or r r Oa 
ou N+v0u (N + v)? 1 Op 
(M+u-—r) + _ + =() (2) 
Or r Oa r p Or 
ey N+ Oo 
(i +u-—r) t + 
or r Oa 
(M+ u)(N + 7) 1 Op 
t = (3 
r pr Oa@ 
Op V+u70p 
(lJ +u—r) / + 
or r Oa 


Oo N+v70 
cyT(M+u-—r) : + — | er (4) 
or r Oa 


Eliminating p and p between these equations one obtains 


Ou (MV +u—r)(N + 7) Ou 
(M+ u — r)? — c?] ' ~ 
or r Oa 
~_M+u4 : Ov 
c? +(M+u—r)(N +02) t 
r or 
(N + v)? — c? Ov : 
+(N+v27)? =0 (5) 
r Oa 
It can be shown that the system of Eqs. (1)-(4) with c? = yp/p 
is hyperbolic for (17 + u — r)? + (N + v)? > c? and elliptic for 


(M+u-—r)?+(N 4+ 0)? < c?. 

Considering that the whole field suffers a uniform expansion in 
time about the corner, one may regard r to represent the ‘‘phase 
velocity” of a point having fixed property. Then it is convenient 
to think of the “‘local field velocity’’ having components 7 + 
u—rand VN +v. It is the flow velocity relative to an observer 
traveling with the point in question. The flow field changes its 
character when the local field velocity exceeds the local speed of 
sound. 

From now on we consider the case when N = O and 4, v are 
small as compared with 14. Then Eqs. (1)-(4) can be approxi 


mated as 


+ The author is indebted to Professor Walker Bleakney of Princeton Uni 


versity for valuable discussions and for sending him several papers 
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(M — r)(Ou/Or + (Op/dor) = 0 f 
}Ir(.M — r)] (Ov/or)} + Mi Op/da 0 7 

}[r(. — r)] (Op/or)} r(Ou/dor) + M+ u 
(Ov/Oa@ xX 


Elimination of p between Eqs. (6) and (7) results in 


O*u O% 
(W--r — 2(M —1r) —rrM-r 
Oroa Or Or 
For r # AM, this equation can be integrated once with respect to 
to give 


Ou/Oa — rov/or — v = fla g 


Here f(q@) is an arbitrary function of a, and it is to be noted that 


the vorticity is given by —f(a)/? Elimination of p and wu fron 
Eqs. (6), (8), and (9) gives a second-order equation in v: 
O*; 0 
1 — (M — r)?|r? + [38 — 201 — 1r)?]r 
or? or 
O77 
fla 0 (10 
Oa’ 
If we put v = v* + g(a@) and let g"(a) + g(a) + fla 0, ther 
the equation for v* becomes 
1—-(M—r)y]r + (3 — 20M — r)*]r — 
or? Or 
0* 
0 (11 
Oa- 


Integrating Eqs. (9) and (7) with respect to @ it will be see 


after some considerations that the solutions are given in terms of 


* 


and an arbitrary function g(@) as follows: 


v*|da — g(a) — M 12 


V/ ‘ a da = 
US ea)da (14 


in this approximation the equation becomes elliptic in the circle 


r< M+ 1for WV < 1orinthering region: M-1l<r<M4 
for J > 1. It may be noted that all the characteristics touch 
the circles, = lorr= M-— 1 


(3) SOLUTIONS FOR THE CASE OF REGULAR REFLECTION 


To solve Eq. (11), we assume v* in the form 
v* = Z(r) sin (kra/a 15 
and make use of a substitution 
g = (M? — 1)/2r — (M — 1)/2 it 
Then we get the following equation 
aZ ( l | ) dZ 
s(1 — 3) t — _ : mm 
dz? 2M +2 Vw? — 1 dz 
R22 /ay? — 1 a 
Zz () 
WV? — 1 


find that Z is given in terms of hypergeometri 


Thus, we 
tions as follows: 


Z= B, Flay, bh . tas 3 + Bf — F(ay — cy 4 


= B.(—=- kK Flay, ap — ¢ 
P(—s)~” 


where FP; — &;, are constants and 


Note 
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P,/P; = 0.2. The line BD must coincide with the circle r 

M+ 1. If one takes the line OA as a = 0, the constant ap should 
7 be taken as r — 6soas to make v = 0on OAandon OC. Then, 
‘ —1/2M +2 | for any integer k, a, > 0 and & < 0 for Wf < 1. Since v must 
vanish at the origin, the second term in Eq. (19) must be dis 


To settle the problem, we consider hereafter the case of regular 
carded; g(a) = 0 by the same reason 


oe reflection as shown in Fig. 1, in which the reflected shock is di : bai : 
ided into two parts: a straight part AB after which a region Consequently, by superposition we find that 
{BD of uniform flow follows, and a curved part BC which bounds a. > we _ oS l 
the region BCOD of pseudo-stationary flow If we choose the oe Bs Ax(—2) FY ai, a Ck + 1, de ee ) : 
gint D as the standard point, the various quantities with suffix oon 
wre equal to those immediately behind the reflected shock at sin : 21 
rect fe { and can be obtained by the simple theory of regular reflection Ss 
Then Vis less than unity except for strong incident shocks with and the expression valid near ry = J + Lis scen to be 
qQ 
fed that 
‘ fron 
4a." Ml —c)F (a — 0 if ee 
r>—6LI(1 — drMa—ct]l 
Mc — 1)M(a — 0 l 
(—< Fia-—-ct+1, b-—-c+1,2—6¢;2 (22 
Mia)r(c — 0 
fy () 
Note that in this and in the following equations the suffixes & of a;,, by, and c are omitted for convenience. By use of Eqs. (12) and 


14) wand pare obtained as follows: 


Mc —1)M(a —b +1 i 1- w\fra — ora — b + 1) abd 
—s)' Fia — « T lb—¢ T 1l2- Cc; gs t Dog A 


12 M(a)l(c — b 2 M1 — bMa—-ct+ ri 
13 . Ii(c — 1)I(a — b+ 1 : 
Fla —- l, ¢ | g (c -— 1 — Fia-—-c+1,6—-c+1,2—¢ 
ria)l(e — b 
Mc —1)Ma-b+1)(a-—-ce+1)(6—c4+1 : | 
l4 ° 3 Filia —c+2,6-—-c4+2,38 -—¢ 24 
I(a)l(c — 0 2 — ¢ ( 
circle ¢ 
~ rT =—*@ kra j . ] 
+1 coust SoA cos WU(-: 'Fla,a—-c+la-b+1; ) = 
_— kr r— 6 ( - 
touch 
i-_ sai : ] ala ( ] 
VW aq —2 ™ F a,a-ctla-b ] — * 
9 ss 
2 \ Z aay 
> 1\ jt 
(— -F (0 +la-c+2,a-60+4+2 ) { 25 
9) a = kra Ml — c)l(a — b + l 
>A cos Sy F(a, b, ¢; 5) 4 
— kn e— sd M1 — b)Ma—ct+l 
Mic — 1)l(a — b l , l V 
If V (—2Z Fla — «¢ lLlb—« 2 - Z — \ x 
l(a I'(< — ~D 2 
r(1 — c)l(a — b + 1) ab Mc — 1)r(a — 6 l 
lbob+i,c+ 1 x 
Mil — b)ra—cetl) « Ma)t(c — b 
Mic — 1l)Ma—b + 1 a—-¢ l(b —<¢ l 
- - F(a — 1lb—« 12-—c; 2) + x 
M(ajl(e — b 2-— 
(—:)! F(a—-c+2,b6-—-c+2,3—C¢; z (26 
Phe constants .1, are to be determined so as to make u,v, and — of experimental results when available. Considering the fact that 


assume Vi »s on BC which are given by the shock rela : ; , . : 
ete “ the alues ¢ ‘ 2 & if 4 entropy is constant along a line which is tangential at any point 
tions. Considerable difficulties arise, however, inasmuch as the 
(19 form of the reflected shock is not known beforehand. One will 


have to rely upon some methods of trial and error, or make us¢ tions of velocity and pressure 


to the local field velocity, we can find density from the distribu 
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Fic. 2. Velocity distribution on the wedge. 





The singularity of rapid expansion, though it is unfamiliar in 

the stationary flow, has its verification in the curves of pressure 

Fic. 1. Configuration and notation of regular reflection. distribution on the wedge which were obtained by White? making 
use of wedges with 6 = 30° and 38°. 

(4) Stncucarity oF Raprp EXPANSION If one reduces the vertex angle of the corner for fixed intensity 
of the incident shock the triangular region ABD of uniform flow 
reduces and finally toa point. Then, the singularity of rapid ex- 
pansion occurs immediately behind the compression by reflected 
shock at point A. It may well be possible that the singularity 
gives influence on the ordinary relations of regular reflection. It 


A few qualitative descriptions of the features of the flow field 
were obtained from preliminary computations making use of three 
terms in the expressions (21)—(26) for a case of P)/P; = 0.8 and 
6 = 45°. The velocity distribution on the wedge is shown sche- 
matically in Fig.2. The tangent of the curve is vertical at D. It 


(2M+2) in the expression (24), might further be responsible, at least partly, for the well-known 


is due to terms (—z)~“ = (—:z)! 
the exponent of which is independent of k. 
The corresponding singularity appears in the pressure dis- Lastly, the writer wishes to express his sincere thanks to 


discrepancies between theory and experiment for Mach reflection. 


tribution on the wedge which has a similar character as the curve Professor S. Tomotika for his kind inspection of the paper. 
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